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Resume: Ce travail est consacre a l’etude de la commande d’une machine asynchrone alimentee par un 
onduleur triphase commande par la MLI a elimination selective d’harmoniques en temps reel. Cette MLI 
a pour but d’eliminer les harmoniques ciblees et les inconnus sont les angles de perturbation. Pour calculer ces 
angles en utilisant l’algorithme d’optimisation par les essaims particulaires. La methode des tableaux l’une des 
plus connue pour une implementation temps reel de la MLI a elimination selective d’harmoniques. 
La nouvelle MLI a elimination selective d’harmoniques temps reel optimisee, e’est-a-dire la generation des 
instants de commutation est basee sur l’intersection d’un signal sinusoidal de reference et un signal dit 
porteuse sinusoidal modifiee. Ce principe est similaire au principe de la MLI triangulo-sinusoldale. 
Pour assurer une generation temps reel du signal de la porteuse modifiee avec une variation continu de 
l’indice de modulation on propose l’exploitation de l’interpolation polynomial de Lagrange. La commande 
vectorielle indirecte a flux rotorique oriente de la machine asynchrone avec des regulateurs PI flous presente 
une meilleure application pour notre nouvelle MLI. La validation experimental des differents resultats 
obtenus en simulation a ete effectuee sur un banc d’essai a base de dSPACE 1103. 

Mots cles : MLI a elimination selective d’harmoniques temps reel, angles de perturbation, algorithme des 
essaims particulaires, porteuse sinusoidal modifiee, interpolation polynomial de Lagrange, machine 
asynchrone, commande vectorielle, regulateurs PI flous. 


Summary: The present research work is dedicated to the control of an induction AC machine fed by a three 
phase inverter, which is controlled using a selective harmonic elimination PWM real time. This PWM signal 
aims to eliminate targeted harmonics, the unknown parameters consist of the perturbation angles. While the 
particle swarm optimisation (PSO) algorithm is adopted for determining the numerical values of these angles, 
the well-established look up table method is used for a real time implementation of the selective harmonic 
elimination PWM. For this novel optimisation real time selective harmonic elimination PWM, the switching 
instants are generated by comparing a sinusoidal reference signal with a modified sinusoidal carrier signal. 
This approach is similar in principle to the sine-triangle PWM. To address real time requirements and such to 
ensure a real time generation of the modified sinusoidal carrier signal with a varying modulation index, a 
Lagrange polynomial interpolation has been successfully used. The indirect field oriented control technique 
using fuzzy controllers for controlling the AC induction machine revealed to be a better choice for our new 
PWM. Different simulation results have been validated using a dSPACE 1103 platform. 

Keywords: Selective harmonic elimination PWM real time, perturbation angles, particle swarm optimisation 
algorithm, modified sinusoidal carrier signal, Lagrange polynomial interpolation, AC induction machine, 
field oriented control, PI fuzzy controllers. 
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Introduction generale 


1 


Aujourd’hui, tous les secteurs industriels, font de plus en plus l’appel a des sources 
d’energie statiques, plus fiables, moins complexes et moins polluantes, pour assurer 
le fonctionnement des differents equipements. Le developpement des nouvelles structures 
des convertisseurs statiques grace au perfectionnement des composants de puissance 
commandables a l’ouverture et a la fermeture tels que GTO et IGBT, comme les onduleurs 
multi-niveau, les convertisseurs matriciels et les onduleurs en Z (Z-source Inverter) 
[Bag-99], [Cas-05], [Bou-08], [Lal-09], [Lou-10], [Ish-12], Tous les progres dans le 
domaine de la microinformatique se traduisent par la creation d’une generation des 
calculateurs numeriques plus puissant et rapide, ce qui permet la facilite de 
T implementation des differents algorithmes de commande, il s’agit de la carte dSPACE 
(digital Signal Processing and Control Engineering) et la carte FPGA (Field Programmble 
Gate Array). Tous ces progres ont permis d’envisager des applications a vitesse variable de 
maniere efficace ou les associations de machines electriques et de convertisseurs statiques 
trouvent de plus en plus d’ applications en moyenne et forte puissance [Sun-94], [Ali-08], 
[Woe-08], [Rus-10], 

L'amelioration de la forme de la tension de sortie des convertisseurs est un axe de 
recherche tres riche. L’objectif d’ameliorer la tension de sortie de l’onduleur peut etre 
atteint soit par la modification de la structure de l’onduleur, soit par le choix de la strategic 
de sa commande. II existe plusieurs strategies de commande par la modulation de large urs 
d’impulsions (MLI), les plus connues sont la modulation triangulo-sinuso'rdale, 
la modulation par hysteresis et la modulation vectorielle qui sont caracterisees par une 
frequence de commutation tres elevee [Roc-96], [Bow-06]. Pour les applications de 
moyenne tension avec un onduleur a deux niveaux, les composants de puissance sont 
limites en frequence et en tension inverse. Pour remedier cet inconvenient, la technique de 
MLI par elimination selective d’harmoniques (MLI-ESH) est la plus adaptee par l’avantage 
de faible frequence de commutation done de faible pertes par commutation, ce qui permet 
d’augmenter la duree de fonctionnement des composants de puissance [Jas-05], [Hui-06], 
[Fan-09], [§te-12], 

Le concept de MLI-ESH a ete propose dans les annees 1960 par Trumbull, et puis 
developpe par Patel et Hoft en 1973 [Pat-73], [Pat-74]. Son principe est de choisir un profil 
de la tension de sortie dans le but d’eliminer un certain nombre d’harmoniques et d’imposer 
la valeur desiree du fondamental [Moh-95]. La MLI-ESH possede deux degres de liberte 
pour la forme soit bipolaire ou unipolaire, de plus on a trois degres de liberte pour la 
symetrie e’est-a-dire la symetrie par rapport au quart et a la demi- periode, par rapport a la 
demi-periode et pas de symetrie done on a six degres de liberte. La formulation 
mathematique de cette technique consiste a etablir un systeme d’equations, dont ses termes 
sont obtenus a partir de la decomposition en serie de Fourier du profil de la tension 
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souhaitee ou les inconnus sont les angles de commutation ou bien les instants de 
commutation, la nature de ce systeme d’ equations est non lineaire a fonctions 
trigonometriques [Joh-04], [Dah-07], [Fei-10]. De nombreuses methodes mathematiques 
sont proposees dans la litterature pour la resolution de ce systeme d’equations, par exemple 
Newton-Raphson, Simplex, la theorie resultante et les polynomes symetriques, les 
algorithmes genetiques (GA), colonies de fourmis (CF) et les essaims particulaires (PSO : 
Particle Swarm Optimisation) [Sun-92], [Joh-05], [Kin-07], [Bar-08], [Dah-10], [Sar-10]. 

La methode classique ou standard pour l’implementation de la MLI-ESH c’est 
la methode des tableaux (Look up table), dans la pratique les tableaux sont des memoires 
electroniques, sa capacite depend du nombre d’harmoniques a eliminer et le pas 
d’echantillonnage de l’intervalle de l’indice de modulation. Recemment, plusieurs travaux 
essaient de proposer des versions de la MLI-ESH sans passer par l’etape de stockage des 
angles de commutation dans des memoires [Phi-04], [Jas-07], [Hir-11]. Le principe de la 
methode proposee dans [Hir-12.1], [Hir-12.2] est inspire de la MLI triangulo-sinuso'rdale, 
au lieu d’utiliser une porteuse triangulaire, on utilise une porteuse a forme sinuso'rdale 
modifiee, cette derniere est comparee avec un signal de reference, par consequence la 
generation des signaux de commande. Les resultats presentes dans ces travaux sont limites 
du point de vue controle de la valeur du fondamental et du nombre d’harmoniques a 
eliminer. D’autre methodes de generation des tensions desirees par la MLI-ESH en temps 
reel ; exigent le passage par l’etape soit de simplification ou d’ approximation des 
trajectoires des angles de commutation calculees par la MLI-ESH. La strategic MLI-ESH a 
echantillonnage regulier (MLI-ESH-ER) [Bow-99], [Bow-01], la MLI-ESH linearisante 
(MLI-ESH-L) [Azl-06], la MLI-ESH simplifiee (MLI-ESH-S) basee sur l’algorithme du 
centre de gravite (Centroid PWM : CPWM) [Vil-0.5] qui est l’amelioration de la MLI-ESH 
avec l’algorithme des surfaces egales (Equal Area Algorithm) [Che-97]. Ces techniques 
assurent un fonctionnement sur toute la plage de variation de l’indice de modulation. 

Dans cette optique, nous nous interessons, au cours de cette these, a 1’ etude de la 
commande des entramements electriques a base de la machine asynchrone triphasee (MAS). 
Les systemes d’ entramements electriques sont beaucoup plus occupes par les MAS qui 
presentent l’avantage d’etre robustes, de construction simple et peu couteuses, mais leurs 
commandes sont plus complexes que celle d’une machine a courant continu. 
Pour la MAS, il est difficile d’obtenir le decouplage effectif entre les deux grandeurs 
de commande qui sont le flux magnetique et le couple electromagnetique. 
L’idee etant d’utiliser la strategic de la commande vectorielle qui permet de rendre 
son comportement identique a celui de la machine a courant continu [Gab-01], [Bou-08], 
[Lar-13], [Lui-13], 
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Le regulateur classique le plus utilise c’est le proportionnel integral (PI) caracterise 
par sa rapidite et sa simplicite de mise en oeuvre. L’inconvenient de ce type de regulateur 
est lie au calcul de ses parametres qui dependent aux valeurs des parametres de la MAS. 
L’intelligence artificielle est employee soit pour optimiser les valeurs de regulateur PI 
[Tis-10], [Far-11] soit pour remplacer ce dernier, parmi les methodes avancees 
on cite la logique floue qui est la plus appliquee [Moa-01], [Rad-04], [Nor-04], 
[Dum-08], [Ara-12]. 

Le contenu de cette these est structure en quatre chapitres presentes comme suit : 

Dans le premier chapitre, nous presentons un etat de Part sur les differentes formes 
et symetries de la MLI-ESH. La resolution des systemes d’equations non lineaires pour le 
calcul des angles de commutation de cette technique est obtenue par l’application d’un 
algorithme d’ optimisation dit les essaims particulaires. 

Le deuxieme chapitre sera consacre au developpement d’une nouvelle technique de 
MLI-ESH en temps reel (MLI-ESH-TR) pour les onduleurs monophases et triphases. Cette 
technique presente notre premiere contribution originale dans cette these, on propose une 
generalisation avec une optimisation d’un signal dit la porteuse modifiee. Les resultats de 
simulation sont confronted aux resultats experimentaux. 

Le troisieme chapitre traite principalement P implementation en temps reel de 
differentes variantes de la MLI-ESH presentees dans la litterature. Nous detaillons les 
algorithmes d’implementation des cinq techniques qui sera la deuxieme contribution 
originale dans cette these car si on revient aux resultats presentes dans ces travaux on ne 
trouve pas des resultats qui demontrent vraiment un changement continu de la tension de 
sortie de l’onduleur. De meme, les resultats de P experimentation sont similaires aux 
resultats de simulation. 

La commande vectorielle a flux rotorique oriente d’une MAS associee a un onduleur 
pilote par la MLI-ESH-TR, fait l’objet du quatrieme chapitre. La focalisation a cette 
application a pour but de montrer la robustesse de notre schema propose de la commande 
vectorielle avec le bloc de generation des impulsions de commande. Ce bloc regoit comme 
entree deux variables independantes qui sont la pulsation et la valeur de la tension desiree 
(indice de modulation). 

Nous terminons par une conclusion generale qui evaluera l'ensemble des resultats 
obtenus et nous proposons des perspectives a ce travail. 




Cha pitre I 

s 

Etat dart de la strategic 
MLI par elimination 
selective d'harmoniques 
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Introduction 

On rencontre plusieurs strategies de commande des onduleurs, telles que la MLI 
triangulo-sinuso'rdale, la MLI par hysteresis, la MLI vectorielle et la MLI par elimination 
selective d’harmoniques. La MLI-ESH (en anglais : selective harmonic elimination pulse 
with modulation SHE-PWM) est reconnue comme etant la plus performante et la plus 
utilisee pour les applications a faible frequence de commutation des interrupteurs de 
l’onduleur ce qui permet de surmonter le probleme des pertes par commutation, ainsi que le 
probleme d’interference electromagnetique ce qui a comme consequence de 1’ amelioration 
du rendement du point de vue la duree de vie de convertisseur [Had-07], [Bar-08], [Bou-08], 

La strategic MLI-ESH a ete proposee pour la premiere fois pendant les annees 1960 
par Trumbull, et puis developpee par Patel et Hoft en 1973 [Pat-73], [Pat-74], elle presente 
une alternative a la strategie traditionnelle la triangulo-sinuso'rdale. Si on revient a la 
strategic de pleine onde, chaque semi-conducteur du meme bras est en etat de conduction 
pendant une demi-periode. Le principe de la MLI-ESH est basee sur l’introduction des 
encoches (coupures) dans la forme rectangulaire de la tension V am (figure 1.1 et 1.2), les 
instants de ces encoches dites les angles de commutation. Le nombre et la position de ces 
angles depend du nombre d’harmoniques a eliminer et l’amplitude du fondamental de la 
tension de sortie respectivement [M oh-95]. 

Pour calculer les angles de commutation il faut d’abord formuler un modele 
mathematique general de l’amplitude des harmoniques, base sur la decomposition en serie 
de Fourier du profil de la tension souhaitee. On obtient un systeme d’ equations non 
lineaires a des termes trigonometrique. Durant plusieurs annees la symetrie par rapport au 
quart et a la demi-periode presente la forme la plus appliquee pour determiner ces angles 
dans l’ordre de reduire la complexity des calculs. Recemment 1’ exigence a cette symetrie ne 
devient pas strictement necessaire du point de vue de la recherche des solutions multiples au 
meme probleme, en raison d’obtenir un faible taux de distorsion harmonique (TDH) 
[Vas-04], [Joh-04], [Jas-05], [Hui-06], [Vas-06], [Dah-10], [Fei-10], 

La strategie MLI-ESH est organisee en deux etapes, premierement la determination 
des angles de commutation par la resolution du systeme d’ equations non lineaires. Pas mal 
de methodes mathematiques sont proposees dans la litterature, en commengant par les 
methodes iteratives (Newton-Raphson, Simplex), d’autres methodes sont basees sur le 
principe de la theorie resultante et les polynomes symetriques [Joh-05], [Jas-07], [Dah-10]. 
Autrement les methodes basees sur des algorithmes metaheuristiques sont de plus en plus 
appliquees aux problemes fortement non-lineaires. II s’agit des algorithmes genetiques 
(GA), colonies de fourmis (CF) et les essaims particulaires (PSO : Particle Swarm 
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Optimisation) [Kin-07], [Sar-10]. Deuxiemement ces angles sont stockes dans des 
memoires pour generer en temps reel [Jas-07]. 

Dans ce chapitre, on etudie les differentes variantes de la strategic MLI-ESH. Ainsi, 
on commencera par elaborer son modele mathematique, sans a priori sur la methode de 
resolution, en utilisant l’algorithme PSO. Ensuite, on fait une etude comparative. 

1.1- Principe de la strategic 

L’onduleur triphase a deux niveaux (Figure 1.1) est constitue de trois bras, chaque bras 
comporte deux interrupteurs en serie. Chaque interrupteur est compose d’un transistor et 
d’une diode montes en tete beche. Pour l’onduleur monophase on utilise seulement deux 
bras. 



triphase 



Figure I. lb Structure d’un onduleur 
monophase 


Figure 1.1 Structure des onduleurs utilises 

Comme dite dans l’introduction, pour une commande en pleine onde et pour une commande 
en modulation de largeurs d’impulsions par elimination selective d’harmoniques, la forme 
de la tension V am de la phase a par rapport au point milieu m de la source est presentee a la 
figure 1. 2a et la figure I.2b respectivement. 




la strategic pleine onde la strategie MLI-ESH 


Figure 1.2 Forme des tensions de sortie 
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La strategie MLI-ESH presente trois types de symetrie : la symetrie par rapport au quart et a 
la demi-periode, par rapport a la demi-periode et enfin pas de symetrie (dissymetrie), en 
plus de la symetrie il est a noter qu’on peut distinguer deux formes : MLI-ESH bipolaire et 
MLI-ESH unipolaire [Joh-04], [Dah-07] comme montre les figures 1.3 et 1.4 : 



Figure 1. 3a La symetrie par rapport au 
quart et a la demi-periode 
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Figure 1. 3b La symetrie par rapport 
a la demi-periode 



Figure 1.3 Differentes symetries de la strategie MLI-ESF1 bipolaire 



Figure I.4a La symetrie par rapport au 
quart et a la demi-periode 
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Figure I.4b La symetrie par rapport 
a la demi-periode 
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Figure I.4c Pas de symetrie 


Figure 1.4 Differentes symetries de la strategie MLI-ESF1 unipolaire 
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1.2- Serie de Fourier 

Le developpement en serie de Fourier est un outil mathematique qui permet de presenter 
une fonction periodique par une somme infinie des fonctions trigonometriques (sinus ou 
cosinus), ces termes trigonometriques ont des frequences multiples de la frequence de la 
fonction. Le but d’utiliser la transformee de Fourier est de determiner l’expression 
mathematique des harmoniques en fonction des angles de commutation selon le nombre 
d’harmoniques a eliminer et l’amplitude du fondamental desiree [Mck-04]. L’expression 
mathematique de la transformee de Fourier d’une fonction periodique quelconque f(t) est la 
suivante: 


f(t)=<.„+2X ,cos(n.2.7r. f 0 .t ) + b n .sin(n.2.7r. f 0 .t ) 

n = 1 

Avec : 

n : nombre entier dans l’intervalle [1, 2, 3, 4,. . oo [ 

fo : la frequence de la fonction f(t) 

a 0 , a n et b n : presentent les coefficients de Fourier 

1 t 0 +T 

a o =- J f(t)dt 

f o 

a n=~ \ f(t).cos(n.2.7T.f 0 .t)dt 
T i 

K =^~ f f (t).sin(n.2.7r.f 0 .t)dt 

T t 

{ fo 


( 1 . 1 ) 


(1.2) 


1.3- Mise en equations 

La mise en oeuvre de cette strategie necessite d’abord a imposer un profil de la tension 
desiree qui possede une certaine forme et symetrie, ensuite on cherche l’expression de 
l’amplitude du fondamental et des harmoniques en fonction des angles de commutation ceci 
est obtenu par le developpement en serie de Fourier de la tension desiree. Ceci conduit a 
etudier les huit cas suivants : 
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1.3.1- La modulation de largeurs d’impulsions par elimination selective 
d’harmoniques bipolaire (MLI-ESH-B) 

Pour cette strategic, la tension de sortie entre le point milieu d’un bras de l’onduleur et le 
point milieu de la source bascule entre la valeur E et -E ou -E et E, done la valeur 
crete-a-crete de cette tension est egale a 2.E. Un inconvenient lie a cette strategic est le 
probleme des interferences magnetiques engendrees par le passage brusque de la tension de 
l’onduleur de la valeur E (-E) a la valeur -E (E). 

Les symetries etudiees par la suite sont les suivantes : 

1.3.1. 1- La symetrie par rapport au quart et a la demi-periode 

Pour annuler N harmoniques et controler le fondamental, on a besoin de (N+ 1) angles a,-, 
avec cette forme, la decomposition en serie de Fourier de la tension ne contient que les 
harmoniques impairs en sinus. L’amplitude de l’harmonique de rang n est donnee par 
1’ expression suivante : 


h 


n 



N + 1 

.cosin.a,) 

i= 1 


(1.3) 


Avec la contrainte sur les angles : 

0 < a 1 < a 2 < a 3 < < a N+1 < nj2 (1.4) 

Comme exemple, nous considerons que N = 4 il est possible d’eliminer quatre harmoniques 
d’ordre impaire. On distingue deux cas : 

Cas 1 : pour un onduleur monophase le rang d’harmonique n prend les valeurs suivantes : 
3, 5, 7, 9. On aboutit au systeme d’ equations non lineaires suivant : 

^ = 1+ 2.(-cos(a 1 ) + cos(« 2 ) - cos(a 3 ) + cos(or 4 ) - cos(a 5 )) - M = 0 
1 2 =1+ 2. (-cos^.cq) + cos(3.« 2 ) - cos(3.a 3 ) + cos(3.cr 4 ) - cos(3.a 5 )) = 0 
< J 3 = 1 + 2. (-cos^.cq) + cos(5.« 2 ) - cos(5.a 3 ) + cos(5.« 4 ) - cos(5.a 5 )) = 0 
J 4 = 1 + 2.(- 005(7.^) + cos(7.ar 2 ) - cos(7.a 3 ) + cos(7.a 4 ) - cos(7.« 5 )) = 0 
I 5 = 1 + 2.(- 005(9.^) + cos(9.« 2 ) - cos(9.a 3 ) + cos(9.a 4 ) - cos(9.a 5 )) = 0 


(1.5) 
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Cas 2 : pour un onduleur triphase (avec un recepteur a neutre isole) les tensions simples et 
composees ne contiennent pas les harmoniques multiples de trois, done le rang 
d’harmonique n prend les valeurs suivantes : 5, 7, 11, 13. On aboutit a un systeme 
d’ equations non lineaires suivant : 

Jj = 1 + 2.(- cos(cq) + cos(a 2 ) - cos(« 3 ) + cos(« 4 ) - cos(a 5 )) - M =0 

I 2 = 1+ 2. (-cos(5.cq) + cos(5.or 2 ) - cos(5.a 3 ) + cos(5.« 4 ) - cos(5.a 5 )) = 0 
< I 3 = 1+2. (-cos(7.cq) + cos(7.« 2 ) - cos(7.a 3 ) + cos(7.or 4 ) - cos(7.a 5 )) = 0 (1.6) 

1 4 =1+ 2.(-cos(ll.a 1 ) + cos(ll.a 2 ) -cos(ll.a 3 ) + cos(ll.a 4 ) -cos(ll.a 5 )) = 0 

1 5 = 1 + 2.(-cos(13.a 1 ) + cos(13.a 2 ) - cos(13.a 3 ) + cos(13.a 4 ) - cos(13.or 5 )) = 0 


Avec : 

M =x.m,/4,E=l (1.7) 

m, : l’indice de modulation (0< m i <1) 


1.3.1. 2- La symetrie par rapport a la demi-periode 

Avec cette symetrie, la decomposition en serie de Fourier de la tension ne contient que les 
harmoniques impairs en sinus et en cosinus. 

II faut noter que chaque harmonique requiert deux angles, e'est-a-dire pour eliminer la 
partie reelle et la partie imaginaire. Pour annuler N harmoniques et controler le fondamental 
on a besoin 2(JV+1) angles. 

L’amplitude de l’harmonique de rang n est donnee par l’expression suivante : 


a = ■ 


b = 


n.n 

4 

n.n 


2.(iV+l) 

^(- 1 ) i+1 .sm(n.a i ) 

Z=1 


2.(1V+1) 

1+ ^(-l)'.cos(n.a / ) 

1=1 


h, =^af+bf 


( 1 . 8 ) 


Avec la contrainte sur les angles : 


0 < a 1 < a 2 < a 3 < < a N+1 < n 


(1.9) 


Pour le fondamental (n = 1) il faut noter aussi qu’il y aura un nombre infini des paires 
(a^bj) et que l’amplitude (fq) et la phase (</>{) du fondamental sont calculees comme suit : 
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'h l =Ja[+b[ 

fa =atag(b l /a ] ) 


( 1 . 10 ) 


Lorsque la phase (dephasage) ne change pas la forme du spectre d’harmoniques ; on peut 
imposer bi = 0 par consequence l’amplitude du fondamental depend uniquement de cq 
[Fei-10], Nous considerons le cas pour N = 4, on a besoin de dix angles. 

Cas 1 : pour un onduleur monophase le rang d’harmonique n prend les valeurs suivantes : 
3, 5, 7, 9. On aboutit au systeme d’equations non lineaires suivant : 

I h = 1- cos((q) + cos(a 2 ) - cos(« 3 ) + cos(a 4 ) - - cos(a 9 ) + cos(a 10 ) - M = 0 

1^ = sin(«, ) - sin(« 2 ) + sin(« 3 ) - sin(« 4 ) + + sin(a 9 ) - sin(a 10 ) = 0 

I 2 = 1 - cosp.cq) + cos(3.a 2 ) - cos(3.« 3 ) + cos(3.« 4 ) - - cos(3.« 9 ) + cos(3.« 10 ) = 0 

I 2 = sinp.cq) - sin(3.« 2 ) + sin(3.« 3 ) - sin(3.« 4 ) + + sin(3.« g ) - sin(3.« 10 ) = 0 

< 2 (IH) 


I 5 =1-005(9.^) + cos(9.« 2 ) - cos(9.« 3 ) + cos(9.« 4 ) - - cos(9.« g ) + cos(9.« 10 ) = 0 

I 5 = sin^cq) - sm(9.a 2 ) + sin(9.« 3 ) - sin(9.« 4 ) + + sin(9.« g ) - sin(9.a 10 ) = 0 

Cas 2 : pour un onduleur triphase (avec un recepteur a neutre isole) les tensions simples et 
composees ne contiennent pas les harmoniques multiples de trois, done le rang 
d’harmonique n prend les valeurs suivantes : 5, 7, 11, 13. On aboutit au systeme 
d’equations non lineaires suivant : 

1, = 1- cos(aj) + cos(a 2 ) - cos(a 3 ) + cos(a 4 ) - - cos(a 9 ) + cos(a 10 ) - M = 0 

I u = sin(«, ) - sin(a 2 ) + sin(a 3 ) - sin(a 4 ) + + sin(a 9 ) - sin(a 10 ) = 0 

f 2| = 1- cos^xq) + cos(5.a 2 ) - cos(5.a 3 ) + cos(5.a 4 ) - - cos(5.a g ) + cos(5.eq 0 ) = 0 

I 2 = sin(5xq) - sin(5.a 2 ) + sin(5.a 3 ) - sin(5.a 4 ) + + sin(5.a g ) - sin(5xq 0 ) = 0 

( 1 . 12 ) 


I 5 = 1- cos(13xq) + cos(13.a 2 ) - cos(13.a 3 ) + cos(13.a 4 ) - - cos(13.a g ) + cos(13.a 10 ) = 0 

I 5 = sin(13xq) - sin(13.a 2 ) + sin(13.a 3 ) - sin(13.a 4 ) + + sin(13.a g ) - sin(13.cq 0 ) = 0 
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1.3.2- La modulation de largeurs d’impulsions par elimination selective 
d’harmoniques unipolaire (MLI-ESH-U) 

Au contraire de la strategie precedente, la tension de sortie entre le point milieu d’un bras 
de l’onduleur et le point milieu de la source bascule entre les valeurs (0 et E ou E et 0) et 
(0 et -E ou -E et 0), ce qui donne une valeur crete-a-crete de cette tension egale a E. 
Un avantage de cette strategie par rapport a la MLI-ESH-B est qu’elle produit moins de 
perturbations electromagnetiques a cause du passage de la tension de sortie de la valeur E 
(0) a la valeur 0 (-E). Les symetries etudiees par la suite sont les suivantes : 

1.3.2. 1- La symetrie par rapport au quart et a la demi-periode 

Pour annuler N harmoniques et controler le fondamental on a besoin de (N+ 1) angles. 
La decomposition en serie de Fourier de la tension ne contient que les harmoniques impairs 
en sinus. L’amplitude de l’harmonique de rang n est donnee par l’expression suivante : 


h = 


n.n 


N+l 


J](- 1 ) ,+1 . cos (n.a,.) 


i=l 


(1.13) 


Avec la contrainte sur les angles : 

0 <a 1 <a 2 <a 3 < < a N+1 < n/2 (1.14) 


On distingue deux cas : 

Cas 1 : elimination des harmoniques du 3 eme jusqu’au 9 eme 

/ : = cos(a 1 ) - cos(a 2 ) + cos(« 3 ) - cos(a 4 ) + cos(a 5 ) - M = 0 
I 2 = cos^.cq) - cos(3.a 2 ) + cos(3.a 3 ) - cos(3.« 4 ) + cos(3.« 5 ) = 0 
< I 3 = cos(5 xq) - cos(5.« 2 ) + cos(5.cr 3 ) - cos(5.a 4 ) + cos(5.« 5 ) = 0 
I 4 = cos^.cq) - cos(7.a 2 ) + cos(7.cr 3 ) - cos(7.a 4 ) + cos(7.a 5 ) = 0 
J 5 = cos(9 .a 4 ) - cos(9.a 2 ) + cos(9.« 3 ) - cos(9.« 4 ) + cos(9.a 5 ) = 0 


(1.15) 
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Cas 2 : elimination des harmoniques du 5 eme jusqu’au 13 eme 

1 3 = cos(a 1 ) - cos(a 2 ) + cos(a 3 ) - cos(cr 4 ) + cos(a 5 ) - M = 0 
I 2 = cos(5.cq) - cos(5.or 2 ) + cos(5.« 3 ) - cos(5.or 4 ) + cos(5) = 0 

< J 3 = cos(7.cq) - cos(7.a 2 ) + cos(7 .a 3 ) - cos(7 .a 4 ) + cos(7.or 5 ) = 0 (1.16) 

1 4 = cos(ll.cq) - cos(ll.a 2 ) + cos(ll.a 3 ) - cos(ll.« 4 ) + cos(ll.a 5 ) = 0 

1 5 = cos^S.aJ - cos(13.a 2 ) + cos(13.a 3 ) - cos(13.a 4 ) + cos(13.a 5 ) = 0 


I.3.2.2- La symetrie par rapport a la demi-periode 

Avec cette symetrie la decomposition en serie de Fourier de la tension ne contient que les 
harmoniques impaires en sinus et en cosinus. 

II faut noter que chaque harmonique requiert deux angles, c'est-a-dire pour eliminer la 
partie reelle et la partie imaginaire. Pour annuler N harmoniques et controler le fondamental 
on a besoin de 2(iV+l) angles 

L’amplitude de l’harmonique de rang n est donnee par l’expression suivante : 


a„ = 


n.n 


2.(1V+1) 

^(-ly.sin (no,) 


i = 1 


ry 2.(N+1) 

<b n = Yj (“ 1 ) <+1 - cos(n.a ; ) 

n.7rl 1 1 

k =V a n +b n 


(1.17) 


Avec la contrainte sur les angles : 

0<a 1 <a 2 <a 3 < < a N+1 < n (1. 18) 

Pour le fondamental (n = 1) il faut noter aussi qu’il y aura un nombre infini des paires 
(cq,£q) et que l’amplitude (fq) et la phase du fondamental sont calculees comme suit : 

K =^o[+b[ 

A = otag(bJa,) 


(1.19) 
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Lorsque la phase ne change pas la forme du spectre d’harmoniques on peut imposer bi = 0 
par consequence l’amplitude du fondamental depend uniquement de cq [Fei-10]. 
Nous considerons le cas pour N = 4, on a besoin de dix angles. 

Cas 1 : elimination des harmoniques du 3 eme jusqu’au 9 eme 

= coslcq) - cos(a 2 ) + cos(a 3 ) - cos(or 4 ) + + cos(a 9 ) - cos(or 10 ) - 2.M = 0 

Jj = -sin(cq) + sin(a 2 ) - sin(« 3 ) + sin(« 4 ) - - sin(« 9 ) + sin(« 10 ) = 0 

I 2 = cos(3xq) - cos(3.a 2 ) + cos(3.« 3 ) - cos(3.« 4 ) + + cos(3.a 9 ) - cos(3.or 10 ) = 0 

I 2 = -sin(3.cq) + sin(3.a 2 ) - sin(3.a 3 ) + sin(3.or 4 ) - - sin(3.a 9 ) + sin(3.a 10 ) = 0 

< I- 2 ( 1 . 20 ) 


J 5 = cos(9.cq) - cos(9.a 2 ) + cos(9.« 3 ) - cos(9.or 4 ) + + cos(9.a g ) - cos(9.a 10 ) = 0 

I 5 = -sin(9.cq) + sin(9.a 2 ) - sin(9.« 3 ) + sin(9.a 4 ) - - sin(9.a 9 ) + sin(9.or 10 ) = 0 

Cas 2 : elimination des harmoniques du 5 eme jusqu’au 13 eme 

= costcq) - cos(a 2 ) + cos(« 3 ) - cos(a 4 ) + + cos(a 9 ) - cos(a 10 ) - 2.M = 0 

Jj = -sh^cq) + sin(a 2 ) - sin(a 3 ) + sin(a 4 ) - - sin(a 9 ) + sin(cq 0 ) = 0 

I 2 = cos(5.cq) - cos(5.or 2 ) + cos(5.or 3 ) - cos(5.a 4 ) + + cos(5.a 9 ) - cos(5.a 10 ) = 0 

I 2 = -sin(5.cq) + sin(5.a 2 ) - sin(5.a 3 ) + sin(5.a 4 ) - - sin(5.a 9 ) + sin(5.a 10 ) = 0 

( 1 . 21 ) 


I- = cos(13.cq) - cos(13.a 2 ) + cos(13.« 3 ) - cos(13.a 4 ) + + cos(13.a 9 ) - cos(13.a 10 ) = 0 

L = -sintlS.cq) + sin(13.a 2 ) - sin(13.a 3 ) + sin(13.or 4 ) - - sin(13.a 9 ) + sin(13.a 10 ) = 0 

v 2 

1.4- Resolution d’un systeme d’equations non lineaires 

Apres la formulation mathematique de la strategie MLI-ESH, on cherche maintenant une 
methode pour trouver la solution de ces systemes d’equations non lineaires a fonctions 
trigonometriques. La difficulty principale de la strategie MLI-ESH est dans le calcul des 
angles de commutation, dont la complexity augmente avec le nombre d’harmoniques a 
eliminer. Plusieurs methodes de resolution des systemes d’equations non lineaires ont ete 
proposees dans la litterature, la methode traditionnelle de Newton-Raphson c’est la plus 
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connue des methodes iteratives dont la convergence a la solution depend fortement de la 
valeur initiale adoptee. Une autre approche est basee sur la transformation des equations a 
fonctions trigonometriques en polynomes. Cette derniere est basee sur la theorie des 
resultantes et les polynomes symetriques. Cette technique permet de determiner des 
solutions multiples sans aucune valeur initiale. D’autres methodes basees sur les methodes 
d’ optimisation intelligente classees comme la methaheuristique, il s’agit des algorithmes 
genetiques, les colonies de fourmis et les essaims particulaires [Mck-04], [Bar-08]. 

L’algorithme des essaims particulaires est une methode stochastique d'une analogie avec les 
comportements collectifs de deplacements d'animaux, introduite en 1995 par James 
Kennedy et Russel Eberhart, cet algorithme est inspire des interactions sociales des 
animaux, denommes particules vivant en essaim, comme les essaims d’abeilles, les vols 
d’oiseaux et les bancs de poissons [Flo-08], [Laz-09]. 

En effet, on peut observer chez ces animaux des dynamiques de deplacement relativement 
complexes, alors qu'individuellement chaque individu a une intelligence limitee et une 
connaissance seulement locale de sa situation dans l'essaim. Un individu de l'essaim n'a 
pour connaissance que la position et la vitesse de ses plus proches voisins. Chaque individu 
utilise done, non seulement sa propre memoire mais aussi l'information locale sur ses plus 
proches voisins pour decider de son propre deplacement. Des regies simples telles que ; 
aller a la meme vitesse que les autres, se deplacer dans la meme direction ou encore rester 
proche de ses voisins sont des exemples de comportements qui suffisent a maintenir la 
cohesion de l'essaim. L'intelligence globale de l'essaim est done la consequence directe des 
interactions locales entre les differentes particules de l'essaim [Sli-09], [Laz-09], [Lef-12]. 

Alors, L'algorithme PSO consiste a suivre les etapes suivantes : 

1/ Creer une population initiale avec une distribution aleatoire de particules et de leurs 
vitesse s ; 

2/ Calculer le vecteur de vitesse pour chaque particule, en utilisant l’experience de chaque 
particule et l’acquisition de vitesse de la population ; 

3/ Mise a jour de la position de chaque particule en utilisant le vecteur de vitesse et la 
position precedente ; 


4/ Evaluer la convergence. 
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Les equations formalisant le mouvement des particules sont [Boub-12] : 

vf +1 = wvf + c 1 r 1 (p i - xf ) + c 2 r 2 (g i - xf ) (1.22) 

xf +1 = xf+vf +1 (1.23) 

k : le numero d’iteration 

xf : la position de la particule z a l'iteration k, 

vf : la vitesse de la particule z a l'iteration k, 

c, , c 2 : les coefficients d’acceleration, 
q , r 2 : des nombres aleatoires dans l’intervalle [0,1], 
w : le coefficient d’inertie. 

Enfin p, est la meilleure position trouvee par la particule et g, est la meilleure position de la 
population. 

La Figure 1.5 presente un schema vectoriel des equations 1.22 et 1.23 ci-dessus [Ebb-12]. 



Figure 1.5 Representation vectoriel de l’algorithme PSO 

L’application de l’algorithme PSO dans notre probleme d’ehmination d’harmoniques 
consiste a chercher le minimum de la fonction objective exprimee par l’equation 1.24 et 1.25 
pour une symetrie par rapport au quart et a la demi-periode et par rapport a la demi-periode 
respectivement : 

F(a 1 ,a 2 ,a 3 ,a 4 ,a 5 ) = 0.5(1; + 1 2 + 1 3 2 + I 4 + 1 2 ) (1.24) 


E(cq, cc 2 , oc 3 , oc 4 ,oc 5 , cc 3 , ct 2 , ct 3 , otg, cqg) — 0.5(1^ + 1 1 + l 2l "t k 2 + 2 + 1^ + I 4 ^ + 1^ + 1 ^ ) (1.25) 
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Etape 1 

Introduire toutes les donnees de l’algorithme telles que : la taille de l’essaim, le nombre 
maximal d’iteration et les coefficients q, c 2 , q, r 2 et w. 

Etape 2 

- Initialisation de la position de chaque particule avec des angles aleatoires verifiant 
l’equation 1.4 (1.9) ; 

- Initialisation de la vitesse de chaque particule par un zero ; 

- Initialisation de la meilleure position p ; de la particule a sa position initiale ; 

- Evaluation de la meilleure position de l’essaim g { , par la recherche du minimum de la 
fonction objective en utilisant l’equation 1.24 (1.25). Ce minimum en realite presente la 
meilleure position p, de la particule. 

Etape 3 

Refaire jusqu'a le nombre maximal d’iteration 

- Mise a jour de la vitesse de la particule par l’equation (1.22) ; 

- Mise a jour de la position de la particule par l’equation (1.23) ; 

- Evaluation de la meilleure position de l’essaim g { , par la recherche du minimum de la 
meilleure position p ; de la particule. 

Dans la derniere iteration, g t presente la meilleure solution trouvee c’est-a-dire les 
meilleurs angles de commutation trouves qui respectent la condition de l’equation 1.4 (1.9). 

1.5- Resultats de simulation 

Dans tous ces cas etudies, on adopte les parametres de l’algorithme PSO pour determiner la 
solution des angles de commutation. La taille de l’essaim est fixee a 100 particules, 
le nombre maximal d’iterations est de 700, les coefficients q = c 2 = 2 et w = 0.9. 


Les angles de commutation sont determines de faqon a annuler les harmoniques de rangs 
trois, cinq, sept et neuf dans le cas d’un onduleur monophase et les harmoniques de rangs 
cinq, sept, onze et treize dans le cas d’un onduleur triphase. 




Chapitre I : Etat d’art de la strategic MLI par elimination selective d’harmoniques 


17 



0 0.1 0.2 0.3 04 05 06 0.7 0.8 09 1 

m, 

Figure 1. 6a Cas 1 
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Figure 1. 6b Cas 2 


Figure 1.6 Trajectoire des angles de commutation en fonction de m, pour une 
symetrie par rapport au quart et a la demi-periode (Bipolaire) 



Figure I.7a Cas 1 


Figure I.7b Cas 2 


Figure 1.7 Trajectoire des angles de commutation en fonction de m, pour une 
symetrie par rapport a la demi-periode (Bipolaire) 


Pour chacun des cas etudies, les solutions donnant les angles de commutation sont 
presentees aux figures 1.6, 1.7, 1.8 et 1. 9a. Elies montrent la variation de ces angles en 
fonction de m ; . L’ analyse de ces figures montre qu’on a des solutions le long de l’intervalle 
de variation de m, [0 , 1], on a aussi une variation quasi lineaire dans des intervalles bien 
determines. II apparait dans la figure 1. 9b une discontinuity dans l’intervalle ]0.6 , 0.7[ 
ce qui signifie pas de solution possible. 
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Figure 1. 8a Cas 1 
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Figure 1.8 Trajectoire des angles de commutation en fonction de m, pour une 
symetrie par rapport au quart et a la demi-periode (Unipolaire) 
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Figure 1. 9a Cas 1 



Figure 1. 9b Cas 2 


Figure 1.9 Trajectoire des angles de commutation en fonction de m, pour une 
symetrie par rapport a la demi-periode (Unipolaire) 


Afin de confirmer la precision des solutions des angles de commutation obtenus par 
l’algorithme PSO, les allures des differentes tensions et leurs spectres d’harmoniques sont 
illustrees aux figures 1. 10 et 1. 11 pour le cas de la modulation bipolaire et aussi aux figures 
1.12 et 1.13 pour le cas de la modulation unipolaire. Dans tous les cas les parametres adoptes 
sont : 

La tension continue a l’entree de l’onduleur est egale a 100 V, la frequence de la tension de 
sortie f= 50Hz et l’indice de modulation m ,• = 0.70. 


Tension V am (V) Tension V ab (V) 
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Figure 1.10a Cas 1 
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Figure 1. 10 Differentes tensions de sortie et leurs spectres d’harmoniques pour une 
symetrie par rapport au quart et a la demi-periode (Bipolaire) 







Tension V am ( V ) Tension V ab (V) 
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Figure 1. 11a Cas 1 









Tension V am (V) Tension V ab ( V ) 
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Figure I.12a Cas 1 
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Figure 1.12 Differentes tensions de sortie et leurs spectres d’harmoniques pour une 
symetrie par rapport au quart et a la demi-periode (Unipolaire) 
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Figure 1.13 Differentes tensions de sortie et leurs spectres d’harmoniques pour une 
symetrie par rapport a la demi-periode (Unipolaire) 
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Les figures 1.14 et 1.15 presentent la variation de l'amplitude du fondamental en fonction de 
l’indic de modulation pour les differents symetries de la MLI-ESH. On remarque que 
l'amplitude du fondamental est proportionnelle avec l’indice de modulation, done la 
MLI-ESH permit un reglage lineaire de l’amplitude du fondamental de la tension de sortie 
de l’onduleur dans l’intervalle de variation de m, [0 , 1], 

L’analyse spectral des differentes tensions montre clairement que les harmoniques d’ordre 
ciblees (trois, cinq, sept et neuf) et (cinq, sept, onze et treize) sont effectivement disparues 
tout en assurant que le fondamental prend sa valeur desiree. Ces resultats confirment 
l’efficacite de l’algorithme PSO du point de vue precision du calcul des valeurs des angles 
de commutation. 

Le Tableau 1.1 presente les valeurs du taux de distorsion harmonique (TDH) base sur la 
formule suivante : 


Cas monophase : TDH = 



Cas triphase : TDH = 



(1.26) 


(1.27) 


Une comparaison des huit cas etudies chacun un caracterise par leur valeur de TDH pour un 
indice de modulation egale 0.70 est reportee au tableau 1.1. 



Type de symetrie 

TDH % 

MLI-ESH 

Bipolaire 

MLI-ESH 

Unipolaire 

Monophase 

par rapport au quart et a la demi-periode 

150.57 

78.51 

par rapport a la demi-periode 

150.40 

78.46 

Triphase 

par rapport au quart et a la demi-periode 

92.34 

46.87 

par rapport a la demi-periode 

93.41 

27.23 


Tableau 1.1 Valeurs du TDH 

Les tensions qui ont une symetrie par rapport a la demi-periode ont un TDH presque 
le meme au cas d’une symetrie par rapport au quart de la periode pour la MLI-ESH-B. 
De plus la methode MLI-ESH-U nous a permis la diminution du TDH a environ la moitie 
de la methode MLI-ESH-B. Un autre point concernant la resolution des systemes 
d’equations non lineaires. Une symetrie par rapport a la demi-periode necessite un nombre 
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double d’inconnus ce qui provoque un temps de calcul plus que la symetrie par rapport au 
quart de la periode, pour le meme nombre d’harmoniques a eliminer. De plus on trouve que 
le changement de l’etat ferme a l’etat ouvert pour les interrupteurs du meme bras soit le 
double pour la strategic MLI-ESH-B par rapport a la strategic MLI-ESH-U. 

L’ implementation de la MLI-ESH-U est plus compliquee car la commande des 
interrupteurs du meme bras se fait d’une maniere independante au contraire de la strategic 
MLI-ESH-B, ou les interrupteurs sont commandes d’une maniere complementaire. 

Conclusion 

Dans ce chapitre, la commande des onduleurs monophases et triphases basee sur la MLI 
par elimination selective d’harmoniques est presentee. L’objectif est d’eliminer les 
harmoniques indesirables d’ordres specifiques et d’imposer le fondamental a une valeur 
desiree. 

La determination des instants (angles) de commutation des interrupteurs necessite 

la resolution des systemes d’ equations non lineaires a fonctions trigonometriques. 

La methode de resolution utilisee rend le probleme d’elimination selective d’harmoniques 
a un probleme d’ optimisation, il s’agit de 1’ optimisation par les essaims particulaires. 
Cette derniere a permis d’obtenir toutes les solutions des cas etudiees sur tout l’intervalle de 
variation de l’indice de modulation sauf pour un seul cas. Cette strategic permet aussi, non 
uniquement d’annuler les harmoniques non desirees mais d’assurer un bon controle de la 
valeur du fondamental. De plus, la tension produite a la sortie de l’onduleur (monophase ou 
triphase) par la MLI-ESH-U offre un taux de distorsion harmonique plus faible que celui de 
la MLI-ESH-B. 

L’ implementation en temps reel de la strategic d’elimination d’harmoniques exige 

le stockage des angles de commutations obtenus en utilisant la methode des tableaux 
(en anglais luck-up table). Pratiquement cette methode est basee sur l’usage des memoires 
electroniques, la capacite de cette derniere depend a la fois du nombre d’angles a utiliser 
(relative au nombre d’harmoniques a eliminer) et le pas d’echantillonnage de l’intervalle 
de l’indice de modulation, ce qui rend la realisation de cette strategic complexe et couteuse. 

Le chapitre suivant sera dedie au developpement d’une nouvelle approche de la strategic 
d’elimination selective d’harmoniques similaire a la technique triangulo-sinuso'rdale. 
Au lieu d’utiliser une porteuse triangulaire, on utilise une porteuse sinuso'rdale modifiee. 




C ha pitre II 

MLI par elimination 
selective d'harmoniques 
en temps reel 
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Introduction 

Depuis plusieurs annees, de nombreuses strategies de commande des onduleurs ont 
ete developpees. La MLI triangulo-sinuso'idale est l’une des techniques la plus connu. Son 
principe general consiste a u t iliser les intersections d’une onde de reference generalement 
sinuso'rdale, avec une onde de modulation (porteuse) sous forme triangulaire. La frequence 
de cette derniere est plus elevee que celle de la reference. La MLI-ESH detaillee dans le 
chapitre precedent presente une alternative de la MLI triangulo-sinuso'idale du point de vue 
diminution de la frequence de commutation [Jas-05]. La difference entre ces deux 
strategies est la maniere de generation des impulsions de commande des interrupteurs. 
Dans la MLI triangulo-sinuso'idale les instants de commutation sont generes a chaque 
instant et en temps reel par une comparaison entre la reference et la porteuse. Par contre la 
MLI-ESH consiste a determiner les angles de commutation par la resolution d’un systeme 
d’equations non lineaires, les valeurs numeriques obtenues des angles sont stockees dans 
des memoires. La lecture de ces angles se fait en temps reel [Phi-04], [Jas-07]. 

Recemment la MLI par elimination selective d’harmoniques en temps reel 
MLI-ESH-TR (en anglais real time harmonic elimination : RTHE), a ete proposee comme 
une alternative de la strategic MLI-ESH. Le principe de la MLI-ESH-TR est d’utiliser une 
porteuse a une forme differente de la forme triangulaire [Phi-04], [Jas-07], [Hir-11], 
le but essentiel de cette strategic est d’eviter l’utilisation des memoires. Ceci permet de 
faciliter la methode de generation des d’impulsions de commande, sans aucune contrainte 
en terme de precision c’est-a-dire on garde le meme spectre d’harmoniques et la valeur du 
fondamental desiree. 

Dans ce contexte, les travaux de [Hir-12.1], [Hir-12.2] proposent l’utilisation d’une 
porteuse a forme sinuso'rdale modifiee. Les resultats presentes dans ces travaux sont limites 
car le nombre maximum d’harmoniques a eliminer est trois sans aucun controle du 
fondamental qui devient un inconvenient de ces travaux. 

Par consequent, dans la premiere partie de ce chapitre, l’inconvenient discute 
precedemment est traite avec succes par la proposition d’un modele general de la strategic 
MLI-ESH-TR pour le cas d’un onduleur monophase et triphase. Pour evaluer l’efficacite de 
la methode proposee, la seconde partie sera consacree pour presenter les differents resultats 
de simulation et experimentale de la methode developpee. 
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II. 1- Principe de la strategic MLI-ESH-TR 


Afin de generer les signaux de commande des interrupteurs d’un onduleur par la technique 
MLI Sinus-Sinus, deux signaux sinuso'idaux sont utilises [Hir-12.1]. Un signal de reference 
v ref = sin( 2.n. f.t + 7r) et une porteuse v car = m i .sin(k.2.7r. f.t) comme illustre dans la 


figure II. la. Le signal resultant par la comparaison entre la reference et la porteuse est 
montre sur la figure II. lb, pour m, = 1, k = 9, la frequence f = 50Hz et t represente le temps. 
L’expression generale pour calculer les points d’ intersection tp, en radians (rad) est donnee 
par : 


®, = L71 . , z = [1 : - — — ],0 <q>, < — ,k = 5, 9, 13, 17, 21, ..., n, n+4 

1 (k + (-!)') 2 1 2 


(II. 1) 
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Figure II. la Forme des signaux utilises 
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Figure II. lb Resultats de comparaison 


Figure II. 1 Principe de la strategie MLI Sinus-Sinus 


Le tableau II. 1 montre un exemple pour k = 5, 9, 13, 17, 21. 
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* 

5ti/12 

5ti/16 

7l/4 

<P6 

* 

* 

3n/7 

n/3 

371/11 

<P7 

* 

* 

* 

7ti/16 

7ti/20 

<P8 

* 

* 

* 

4n/9 

471/11 

<P9 

* 

* 

* 

* 

9ti/20 

<PlO 

* 

* 

* 

* 

571/11 


Tableau II. 1 Valeurs des points d’intersection 
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II.2- Mise en equations 

Le signal presente dans la figure II. lb a une symetrie par rapport au quart et a la demi- 
periode (les harmoniques pairs sont nuls). La decomposition en serie de Fourier montre que 
l’amplitude de l’harmonique de rang n est donnee par l’expression suivante : 


h =■ 


rut 


1+2-X (- 1 )'- cos(n.ft.) 


(H-2) 


k et N sont liees par les relations suivantes : 

k = 2(N + 1) + 1 (II. 3) 


Dans le but d’eliminer N harmoniques et de controler le fondamental, on a besoin de (N+ 1) 
angles de perturbation par quart de periode [Mak-14.1]. L’equation (II. 2) devient : 


h 


n 



IV +1 

2 -Z( _1 )'- cos ( n -(^ ~A)) 

i=l 


(H-4) 


Comme exemple, nous considerons que k = 9 done N = 3, il est possible d’eliminer trois 
harmoniques d’ordre impair et de controler le fondamental. On distingue deux cas : 

Cas 1 : pour un onduleur monophase, le rang d’harmonique n prend les valeurs suivantes : 
3, 5, 7. On aboutit au systeme d’ equations non lineaires suivant : 


Jj = 1 + 2.(-cos(^ - /3 X ) + coprp 2 - /? 2 ) - cos(^ 3 - /? 3 ) + cos(^? 4 -/? 4 ) -M =0 
1 2 = 1+ 2.(-cos(3.(ft -A)) + cos(3.(^ 2 -p 2 ))~ cos(3.(^ 3 -#,)) + cos(3 .{rp A - p A )) = 0 ^ ^ 

I 3 = 1+ 2.(-cos(5.(^ - /?,)) + cos(5.(^ 2 -p 2 ))~ cos(5.(^ 3 - pp) + cos(5 ,{cp A - PP) = 0 
h = 1+ 2.(-cos(7 .{(p x - PP) + cos(7 .{rp 2 -p 2 ))~ cos(7 .{rp 3 - pp) + cos(7 .{rp A - PP) = 0 

Cas 2 : pour un onduleur triphase (avec un recepteur a neutre isole), les tensions simples et 

composees ne contiennent pas les harmoniques multiples de trois, done le rang 
d’harmonique n prend les valeurs suivantes : 5, 7, 11. On aboutit au systeme d’equations 
non lineaires suivant : 
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1+ 2.(-cos(^ 1 -/?,) + cos(^ 2 - P 2 ) - cos(^ 3 - A) + cos(^> 4 - A) - M =0 
1 + 2.(-cos(5.(^ - A)) + cos(5.0 2 - A)) - cos(5.(^ 3 - A)) + cos(5.(p 4 - A)) = 0 ^ 

1 + 2.(-cos(7.(ft - A)) + cos(7.(<p, -P 2 ))~ cos(7.(^ 3 - A)) + cos(7.(p 4 ~PJ) = 0 
1+ 2.(-cos(ll.(^ - A)) + cos(11.0 2 -A))- cos(ll.(^ 3 - A)) + cos(ll.(p 4 - A)) ^ 0 

Avec : 

M = K.mjA (II. 7) 

m, : l’indice de modulation (0< m, <l) 

Afin de determiner les valeurs desirees des angles de perturbation pi pour n'importe quelle 
valeur de m, et k, le systeme non lineaire a fonctions trigonometriques decrit en (II. 5) 
et (II. 6) doit etre resolu. 

On va adopter pour resoudre ce probleme l’algorithme d’ optimisation par essaim 
particulaire, en minimisant la fonction objective suivante : 

F (A>A>A> ,Av + i) = 0.5 (A + i 2 2 + Il+ + A + i) (n.8) 

Avec une contrainte sur les angles : 

0 ^ ft ~ A ^ <Pi ~ A ^ - A ^ ^ <Pn+i - Pnm ^ V 2 ( n - 9 ) 

Pour generer le signal desire de commande des interrupteurs obtenu a partir de 
l’intersection entre le signal de reference et de la porteuse, on doit determiner l’expression 
du signal de la porteuse modifiee utilisant les angles de perturbation calculees par 
l’equation (II. 8). Ce signal a une symetrie par rapport au quart et a la demi-periode d’une 
part, et d’ autre part on doit eliminer les harmoniques indesirables et controler le 
fondamental. 

L’expression du signal de la porteuse modifiee proposee par [Hir-12.1], [Hir-12.2] : 

f(@) = A L .sin(k 1 .6) + A 2 .sin(k 2 .9) + (II. 10) 

Avec Ai et k, sont les inconnus. 

La forme proposee de f(9) par [Mak-14.1], [Mak-14.2], [Mak-14.3] : 

N + 1 

f(^) = Z C cA > C, = sin(q,.0) 



(11-11) 
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Avec Ci c’est un constant, 6?. = rp i - [i i ct q, prennent les valeurs impaires dans l’intervalle 

[3, k\. L’avantage de notre formule proposee de la porteuse modifiee est qu’elle a une forme 
lineaire ce qui facilite sa resolution par n’importe quelle methode de resolution des 
systemes d’equations lineaire [Mak-14.1]. 

Pour plus de details on va prendre le cas de k = 9, le systeme d’equations qui permet de 
determiner l’inconnu A, est decrit comme suit : 

f (A) = A sin(qA) + A sin(q 2 A) + A 3 sin(q 3 ^) + A 4 sin(q 4 ^) 
f ( A) = A sin(q,# 2 ) + A 2 sin (q 2 0 2 ) + A s i n (t? 3 ^ 2 ) + A sin(q 4 0 2 ) 

( 11 . 12 ) 

f (A) = A sintqA) + A sin(q 2 A) + A sin(q 3 A) + A sin(q 4 A) 
f (A) = A sin (QiA) + A sin (q 2 A) + A sin (q 3 A) + A sin (q 4 A) 


Avec: 

f (A ) = sin(6>. + n ) , = q>. - j3 i ,i = 1, 2, 3, 4 et qj=3, q 2 = 5, q 3 =7, q 4 =9. 

La strategie proposee peut se resumer dans les etapes suivantes : 

Etape 1 : Choisir la valeur de l’indice de modulation et fixer le nombre d’harmoniques 
a eliminer en relation avec la valeur de k ; 

Etape 2 : Determiner les valeurs des angles de perturbation par la resolution du systeme 
d’equations (II. 5) pour le cas 1 et (II.6) pour le cas 2 en utilisant l’algorithme PSO ; 

Etape 3 : Construire un systeme d’equations utilisant l’equation (II. 11) et les angles de 
perturbation determines a l’etape 2 pour calculer l’inconnu A, ; 

Etape 4 : Introduce les valeurs de A, et q, dans l’equation (II. 11) ; 

Etape 5 : Generer le signal de commande par la comparaison entre le signal de reference et 
le signal de la porteuse modifiee determinee a l’etape 4. 

II.3- Resultats de simulation 

Afin de confirmer la validite de notre strategie proposee, on va presenter des resultats dans 
le cas ou tous les harmoniques sont elimines jusqu’a l’harmonique 7, 11, 15 et 19 pour un 
onduleur monophase, et jusqu’a l’harmonique 11, 17, 23 et 29 pour un onduleur triphase. 
On adopte les parametres de l’algorithme PSO pour determiner les valeurs des angles de 
perturbation. La taille de l’essaim est fixee a 100 particules, le nombre maximal d’iterations 
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est de 700, les coefficients c± - C 2 = 2 et w = 0.9. Le pas de l’indice de modulation 
est fixe a 0.1. 

La tension continue a l’entree de l’onduleur est egale a 100 V, la frequence de la tension de 
sortie f= 50Hz et l’indice de modulation m; = 0.70 

Les valeurs des angles de perturbation pour un onduleur monophase et triphase avec 
k = 9, 13, 17 ,21 sont presentees dans les tableaux II. 2, 3, 4, 5 et II. 6, 7, 8, 9 
respectivement. Dans les resultats obtenus, nous avons constate qu’on a des solutions pour 
chaque pas de l’indice de modulation. Ces solutions prennent des signes negatifs ou positifs 
done on aura un deplacement des points d’intersection originaux soit vers la gauche ou vers 
la droite. 


m ,- 

Pi 

P2 

Ps 

fP 

0.1 

0.0380 

- 0.0581 

0.1162 

- 0.1222 

0.2 

0.0332 

- 0.0453 

0.1025 

- 0.1043 

0.3 

0.0292 

- 0.0315 

0.0898 

- 0.0857 

0.4 

0.0263 

- 0.0167 

0.0784 

- 0.0662 

0.5 

0.0246 

- 0.0005 

0.0687 

- 0.0451 

0.6 

0.0244 

0.0170 

0.0614 

- 0.0217 

0.7 

0.0261 

0.0364 

0.0582 

0.0057 

0.8 

0.0306 

0.0584 

0.0634 

0.0411 

0.9 

0.0398 

0.0851 

0.0897 

0.0970 

1.0 

0.0681 

0.1340 

0.2144 

0.2452 


Tableau II. 2 Valeurs numeriques des angles de perturbation pour k = 9 


mi 

Pi 

P2 

Ps 

P * 

Ps 

Pe 

0.1 

0.0174 

- 0.0286 

0.0528 

- 0.0587 

0.0897 

- 0.0915 

0.2 

0.0151 

- 0.0220 

0.0461 

- 0.0476 

0.0791 

- 0.0791 

0.3 

0.0133 

- 0.0148 

0.0403 

- 0.0356 

0.0692 

- 0.0664 

0.4 

0.0120 

- 0.0070 

0.0356 

- 0.0227 

0.0599 

- 0.0531 

0.5 

0.0112 

0.0016 

0.0323 

- 0.0084 

0.0514 

- 0.0391 

0.6 

0.0112 

0.0110 

0.0309 

0.0075 

0.0442 

- 0.0238 

0.7 

0.0121 

0.0213 

0.0321 

0.0259 

0.0392 

- 0.0063 

0.8 

0.0141 

0.0330 

0.0375 

0.0480 

0.0388 

0.0157 

0.9 

0.0183 

0.0470 

0.0504 

0.0770 

0.0520 

0.0513 

1.0 

0.0422 

0.0884 

0.1325 

0.1743 

0.2457 

0.2624 


Tableau II. 3 Valeurs numeriques des angles de perturbation pour k- 13 
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mi 

Pi 

Pi 

Ps 

p 4 

Ps 

Ps 

P7 

Ps 

0.1 

0.0100 

- 0.0170 

0.0301 

- 0.0345 

0.0507 

- 0.0531 

0.0721 

- 0.0729 

0.2 

0.0086 

- 0.0130 

0.0261 

- 0.0273 

0.0443 

- 0.0440 

0.0637 

- 0.0634 

0.3 

0.0076 

- 0.0086 

0.0228 

- 0.0195 

0.0386 

- 0.0343 

0.0557 

- 0.0538 

0.4 

0.0068 

- 0.0038 

0.0203 

- 0.0109 

0.0337 

- 0.0238 

0.0480 

- 0.0438 

0.5 

0.0064 

0.0014 

0.0187 

- 0.0014 

0.0299 

- 0.0123 

0.0409 

- 0.0334 

0.6 

0.0064 

0.0072 

0.0183 

0.0091 

0.0276 

0.0005 

0.0345 

- 0.0222 

0.7 

0.0069 

0.0136 

0.0194 

0.0211 

0.0276 

0.0155 

0.0294 

- 0.0097 

0.8 

0.0081 

0.0209 

0.0228 

0.0352 

0.0315 

0.0342 

0.0270 

0.0057 

0.9 

0.0104 

0.0294 

0.0298 

0.0524 

0.0428 

0.0595 

0.0326 

0.0285 

1.0 

0.0157 

0.0419 

0.0471 

0.0793 

0.0774 

0.1063 

0.0929 

0.1044 


Tableau II. 4 Valeurs numeriques des angles de perturbation pour k- 17 


m 

Pi 

P2 

Ps 

P4 

Ps 

Pe 

P7 

CO 

Ps 

P10 

0.1 

0.0064 

- 0.0112 

0.0194 

- 0.0227 

0.0326 

- 0.0347 

0.0461 

- 0.0473 

0.0601 

- 0.0605 

0.2 

0.0056 

- 0.0086 

0.0168 

- 0.0178 

0.0283 

- 0.0281 

0.0403 

- 0.0397 

0.0532 

- 0.0529 

0.3 

0.0049 

- 0.0057 

0.0147 

- 0.0123 

0.0247 

- 0.0208 

0.0351 

- 0.0317 

0.0465 

- 0.0451 

0.4 

0.0044 

- 0.0024 

0.0131 

- 0.0063 

0.0217 

- 0.0129 

0.0305 

- 0.0232 

0.0400 

- 0.0371 

0.5 

0.0041 

0.0011 

0.0122 

0.0003 

0.0197 

- 0.0042 

0.0266 

- 0.0138 

0.0339 

- 0.0289 

0.6 

0.0041 

0.0050 

0.0120 

0.0077 

0.0188 

0.0056 

0.0240 

- 0.0034 

0.0282 

- 0.0201 

0.7 

0.0045 

0.0094 

0.0128 

0.0160 

0.0195 

0.0169 

0.0231 

0.0087 

0.0234 

- 0.0104 

0.8 

0.0052 

0.0144 

0.0151 

0.0257 

0.0227 

0.0305 

0.0255 

0.0240 

0.0205 

0.0011 

0.9 

0.0068 

0.0204 

0.0198 

0.0378 

0.0307 

0.0486 

0.0353 

0.0466 

0.0233 

0.0185 

1.0 

0.0093 

0.0273 

0.0279 

0.0521 

0.0457 

0.0715 

0.0597 

0.0804 

0.0513 

0.0578 


Tableau II. 5 Valeurs numeriques des angles de perturbation pour k- 21 


m ,- 

Pi 

Pi 

Ps 

P4 

0.1 

0.3352 

0.0730 

0.5858 

0.2298 

0.2 

0.3085 

0.0940 

0.5671 

0.2502 

0.3 

0.2862 

0.1093 

0.5487 

0.2707 

0.4 

0.2662 

0.1277 

0.5305 

0.2913 

0.5 

0.2476 

0.1465 

0.5125 

0.3124 

0.6 

0.2304 

0.1655 

0.4952 

0.3337 

0.7 

0.2144 

0.1850 

0.4786 

0.3559 

0.8 

0.1999 

0.2051 

0.4634 

0.3791 

0.9 

0.1870 

0.2266 

0.4550 

0.4051 

1.0 

0.1768 

0.2503 

0.4438 

0.4369 


Tableau II. 6 Valeurs numeriques des angles de perturbation pour k = 9 
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m 

Pi 

P 2 

Ps 

P4 

Ps 

Pe 

0.1 

0.2249 

0.0664 

0.3801 

0.1885 

0.5731 

0.3132 

0.2 

0.2075 

0.0777 

0.3685 

0.2020 

0.5599 

0.3273 

0.3 

0.1929 

0.0891 

0.3571 

0.2156 

0.5467 

0.3414 

0.4 

0.1797 

0.1005 

0.3460 

0.2295 

0.5338 

0.3556 

0.5 

0.1673 

0.1120 

0.3353 

0.2437 

0.5210 

0.3700 

0.6 

0.1558 

0.1234 

0.3253 

0.2580 

0.5087 

0.3843 

0.7 

0.1445 

0.1356 

0.3155 

0.2736 

0.4963 

0.3996 

0.8 

0.1341 

0.1477 

0.3071 

0.2899 

0.4847 

0.4153 

0.9 

0.1245 

0.1604 

0.3005 

0.3079 

0.4744 

0.4322 

1.0 

0.1160 

0.1740 

0.2976 

0.3292 

0.4673 

0.4525 


Tableau II. 7 Valeurs numeriques des angles de perturbation pour k- 13 


mi 

Pi 

P2 

Ps 

P 4 

Ps 

Pe 

P 7 

Ps 

0.1 

0.1861 

0.1406 

0.1620 

0.0800 

0.1372 

0.0092 

0.1068 

- 0.0636 

0.2 

0.1758 

0.1411 

0.1540 

0.0904 

0.1304 

0.0186 

0.0955 

- 0.0579 

0.3 

0.1652 

0.1411 

0.1462 

0.1010 

0.1237 

0.0283 

0.0837 

- 0.0527 

0.4 

0.1545 

0.1405 

0.1387 

0.1120 

0.1172 

0.0383 

0.0714 

- 0.0483 

0.5 

0.1435 

0.1394 

0.1317 

0.1236 

0.1109 

0.0487 

0.0583 

- 0.0450 

0.6 

0.1323 

0.1376 

0.1253 

0.1362 

0.1049 

0.0596 

0.0438 

- 0.0434 

0.7 

0.1207 

0.1353 

0.1202 

0.1507 

0.0995 

0.0714 

0.0274 

- 0.0444 

0.8 

0.1087 

0.1325 

0.1182 

0.1698 

0.0948 

0.0843 

0.0077 

- 0.0497 

0.9 

0.0964 

0.1300 

0.1288 

0.2028 

0.0921 

0.0989 

- 0.0170 

- 0.0615 

1.0 

0.0860 

0.1321 

0.1747 

0.2604 

0.0874 

0.1150 

- 0.0417 

- 0.0709 


Tableau II. 8 Valeurs numeriques des angles de perturbation pour k- 17 


m 

Pi 

P2 

Ps 

Ps 

Ps 

Pe 

P7 

Ps 

Ps 

P10 

0.1 

0.1485 

0.1109 

0.1147 

0.0557 

0.2570 

0.1689 

0.2217 

0.1029 

0.1821 

0.0356 

0.2 

0.1398 

0.1103 

0.1074 

0.0640 

0.2530 

0.1797 

0.2165 

0.1108 

0.1720 

0.0391 

0.3 

0.1309 

0.1093 

0.1001 

0.0728 

0.2500 

0.1911 

0.2113 

0.1189 

0.1613 

0.0418 

0.4 

0.1218 

0.1078 

0.0930 

0.0822 

0.2485 

0.2037 

0.2062 

0.1274 

0.1501 

0.0437 

0.5 

0.1126 

0.1060 

0.0866 

0.0931 

0.2496 

0.2176 

0.2012 

0.1363 

0.1379 

0.0445 

0.6 

0.1032 

0.1039 

0.0820 

0.1071 

0.2548 

0.2345 

0.1963 

0.1456 

0.1247 

0.0440 

0.7 

0.0936 

0.1020 

0.0831 

0.1272 

0.2648 

0.2535 

0.1912 

0.1553 

0.1107 

0.0431 

0.8 

0.0844 

0.1013 

0.0953 

0.1549 

0.2755 

0.2725 

0.1856 

0.1653 

0.0968 

0.0432 

0.9 

0.0758 

0.1024 

0.1176 

0.1852 

0.2815 

0.2894 

0.1794 

0.1757 

0.0831 

0.0447 

1.0 

0.0683 

0.1059 

0.1439 

0.2132 

0.2842 

0.3055 

0.1739 

0.1871 

0.0680 

0.0456 


Tableau II. 9 Valeurs numeriques des angles de perturbation pour k- 21 
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Maintenant nous presentons en detail le calcul du parametres A,- pour un taux de modulation 
egal a 0.70. Les tableaux ci-dessous regroupent les differentes valeurs numeriques pour 
creer le signal de la porteuse modifiee. 


k~- 

= 9 

k = 

13 

k = 

: 17 

k = 

21 

9i (rad) 

mi 

0i (rad) 

m) 

9i (rad) 

fm 

9[ (rad) 

fm 

0.3666 

-0.3584 

0.2497 

-0.2471 

0.1894 

-0.1883 

0.1526 

-0.1520 

0.5919 

-0.5579 

0.4275 

-0.4146 

0.3355 

-0.3292 

0.2762 

-0.2727 

1.1199 

-0.9000 

0.7533 

-0.6840 

0.5696 

-0.5393 

0.4584 

-0.4425 

1.2509 

-0.9493 

0.8717 

-0.7654 

0.6770 

-0.6264 

0.5552 

-0.5271 

* 

* 

1.2698 

-0.9550 

0.9541 

-0.8158 

0.7659 

-0.6932 

* 

* 

1.3527 

-0.9763 

1.0317 

-0.8582 

0.8399 

-0.7445 

* 

* 

* 

* 

1.3450 

-0.9746 

1.0765 

-0.8803 

* 

* 

* 

* 

1.4057 

-0.9864 

1.1337 

-0.9060 

* 

* 

* 

* 

* 

* 

1.3903 

-0.9837 

* 

* 

* 

* 

* 

* 

1.4384 

-0.9912 


Tableau II. 10 Nouvelles coordonnees des points d’intersection (Onduleur monophase) 


k~- 

= 9 

k = 

13 

k = 

17 

k = 

21 

9i (rad) 

fm 

9-, (rad) 

fm 

9-, (rad) 

fm 

9i (rad) 

fm 

0.1782 

-0.1772 

0.1171 

-0.1168 

0.0757 

-0.0756 

0.0634 

-0.0633 

0.4434 

-0.4290 

0.3134 

-0.3083 

0.2137 

-0.2120 

0.1836 

-0.1825 

0.6995 

-0.6438 

0.4698 

-0.4527 

0.4688 

-0.4518 

0.3881 

-0.3784 

0.9008 

-0.7838 

0.6242 

-0.5844 

0.5474 

-0.5204 

0.4440 

-0.4295 

* 

* 

0.8126 

-0.7260 

0.8823 

-0.7722 

0.5206 

-0.4974 

* 

* 

0.9470 

-0.8116 

0.9758 

-0.8281 

0.6033 

-0.5673 

* 

* 

* 

* 

1.3470 

-0.9750 

0.9083 

-0.7884 

* 

* 

* 

* 

1.4406 

-0.9915 

0.9871 

-0.8334 

* 

* 

* 

* 

* 

* 

1.3030 

-0.9643 

* 

* 

* 

* 

* 

* 

1.3849 

-0.9827 


Tableau II. 11 Nouvelles coordonnees des points d’intersection (Onduleur triphase) 

On remplace chaque valeur de la paire (#, , f (#,.)) dans 1’ equation (II. 11), par consequence 

on abouti a des systemes d’ equations lineaires de la variable A„ la taille de ces systemes 
varie selon la valeur de k. 

Les valeurs numeriques de q ; pour k = 9, 13, ..., n, n+4 sont definies comme suit : 


k = 9 : qi= 3, q 2 =5, q 3 = 7, q 4 = 9 ; 

k= 13 : qi= 3, q 2 = 5, q 3 = 7, q 4 = 9, q 5 =ll, qe = 13 ; 

k = n: q 3 = 3, q 2 = 5, q 3 = 7, , q(„-q/ 2 = n (n est un nombre entier). 
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Onduleur monophase 

Onduleur triphase 

CD 

II 

k = 13 

k = 17 

k = 21 

CD 

II 

k = 13 

k= 17 

k = 21 

Ai 

-0.16609 

-0.14734 

-0.14663 

-0.14645 

-0.80268 

-0.85243 

-0.14554 

-0.71558 

a 2 

0.11394 

0.02864 

0.02421 

0.02408 

0.50319 

0.61334 

-0.29449 

0.41612 

A3 

-0.34049 

-0.02499 

-0.00535 

-0.00402 

-0.22814 

-0.36382 

0.33996 

-0.27987 

a 4 

0.86026 

0.09098 

0.00540 

0.00112 

0.05671 

0.16943 

-0.38037 

0.29008 

A5 

0 

-0.33841 

-0.02101 

-0.00136 

0 

-0.05593 

0.20722 

-0.39762 

Ae 

0 

0.86367 

0.09061 

0.00427 

0 

0.00992 

0.34049 

0.55270 

a 7 

0 

0 

-0.33864 

-0.02066 

0 

0 

-0.69763 

-0.61356 

Aq 

0 

0 

0.86456 

0.09039 

0 

0 

0.34314 

0.46994 

A9 

0 

0 

0 

-0.33830 

0 

0 

0 

-0.22315 

Ajo 

0 

0 

0 

0.86412 

0 

0 

0 

0.05063 


Tableau 11.12 Valeurs numeriques de A, 

II ne reste que de noter, que pour un systeme triphase on a un dephasage de 2.jt/3 entre les 
trois phases, done les trois phases ont un dephasage de 0, -2.n/3, 2.n/3 respectivement. 
Ainsi, a chaque terme de la porteuse modifiee liee a la premiere phase, la deuxieme phase et 
la troisieme phase on ajoute des termes de correction qui sont : n, 7t-q ; .2.7t/3, 7t+q ; .2.7t/3. 


Un exemple de l’expression du signal de la porteuse modifiee pour le cas d’un onduleur 
monophase et triphase respectivement avec k = 9 et m ; = 0.7 : 


f (d) = -0.16609.sin(3.£ + k) + O.11394.sin(5.0 + n)~ 0.34049.sin(7.£ + n) 
+ 0.86026.sin(9 .0 + x) 


' f[6) = -O.8O268.sin(3.0 + n) + O.5O319.sin(5.0 + n)- 0.22814.sin(7.^ + n) 

+ 0. 05671. sin(9.0 + /r) 

f(0) = -0.80268.sin(3.(^ - — ) + n) + 0.50319.sin(5.(^ - — ) + n) 

3 3 

< -0.22814.sin(7.(^-— ) + n) + 0.05671. sin(9.(#-— ) + ;r) ( IL14 ) 

3 3 

rj r-* 

f{0) = -0.80268.sin(3.(^ + — ) + k) + 0.50319.sin(5.(^ + — ) + n) 

3 3 

- 0.22814.sin(7.(^ + — ) + it) + 0.05671. sin(9.(^ + — ) + n) 

, 3 3 


Avec : 6 = 2.n.\ X 
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II.3.1- Resultats de simulation pour un onduleur monophase 

Pour demontrer les performances de 1’ onduleur monophase commande par notre technique 
proposee pour differentes valeurs de k (k = 9, 13, 17, 21), le signal de la porteuse modifiee, 
la tension V a b et son spectre d’harmoniques sont montres dans les figures II.2 et II. 3. 




t(s) 

Figure II. 2c k = 17 



t(s) 

Figure II. 2d k = 21 


Figure II. 2 Signal de la porteuse modifiee et le signal de reference 
pour un onduleur monophase 


Les resultats de simulation obtenus montrent, que la methode proposee est tres satisfaisante 
pour differentes valeurs de k. Les allures de la porteuse modifiee montrent clairement une 
forme oscillatoire, avec une frequence de k fois la frequence de la reference et ayant une 
symetrie par rapport au quart et a la demi-periode, comme illustre dans la figure II. 2. 
A partir des spectres d’harmoniques des tensions de sortie de 1’ onduleur dans la figure II. 3, 
on constate que tous les harmoniques d’ordre non desire sont nuls pour differentes valeurs 
de k. Ce qui confirme done la faisabilite de notre methode proposee. Les premiers rangs 
d’harmoniques non controlables sont le 9 dme , 13™ e , 17" me et le 21 toe . 


Tension V ab (V) Tension V ab (V) Tension V ab (V) Tension V ab (V) 
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Figure II. 3a k = 9 
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Figure II. 3b k= 13 
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Figure II. 3c k = 17 
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Rang des harmoniques 


Figure II. 3d k = 21 


Figure II. 3 Differentes tensions de sortie et leurs spectres d’harmoniques 
pour un onduleur monophase 
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II.3.2- Resultats de simulation pour un onduleur triphase 

Pour demontrer les performances de 1’ onduleur triphase commande par notre technique 
proposee pour differentes valeurs de k (k = 9, 13, 17, 21), le signal de la porteuse modifiee, 
la tension V am , V as et son spectre d’harmoniques sont montres dans les figures II.4, 
II. 5 et II. 6. 




t(s) t ( s ) 

Figure II. 4a k = 9 Figure II. 4b Ac = 13 



Us) 

Figure II. 4c Ac = 17 



0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 


t(s) 

Figure II. 4d Ac = 21 


Figure II. 4 Signal de la porteuse modifiee et le signal de reference 
pour un onduleur triphase 


Les resultats obtenus s’averent satisfaisants pour differentes valeurs de k. Les allures 
de la porteuse modifiee montrent clairement une forme quasi oscillatoire et ayant une 
symetrie par rapport au quart et a la demi-periode, comme illustre dans la figure II.4. 
A partir des spectres d’harmoniques des tensions de sortie de l’onduleur dans les figures 
II. 5a, 5b, 6a et 6b, on constate que tous les harmoniques de rang has d’ordres non desires 
seront elimines, c’est-a-dire les premiers rangs d’harmoniques non controlables sont le 
13™ 16 , 19 pme , 25 pmp et le 31 dme . 


Tension V am (V) Tension V am (V ) 
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Rang des harmoniques 


Figure II. 5b k= 13 


Figure II. 5 Differentes tensions de sortie et leurs spectres d’harmoniques 
pour un onduleur triphase 






Tension V am (V) Tension V am (V ) 
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II.4- Resultats experimentaux 

Dans le but de verifier experimental ement la strategie d’ elimination selective 
d’harmoniques en temps reel, la figure II. 7a presente la photographie du banc d’essais 
utilise, dont le schema synoptique est donne par la figure II. 7b. Les algorithmes de 
commande sont implementes sur la carte DS1103 TMS320F240 branchee a un PC mene du 
logiciel MATLAB / Simulink. Pour visualiser les signaux de commande et l’acquisition des 
differentes grandeurs mesurees (tension, courant) en temps reel, le logiciel control-desk est 
employe, la fenetre de ce dernier est illustree dans la figure II. 7c [dSP-06], [Mak-14.1], 

La partie puissance est constitute d’un l’onduleur qui est compose de six transistors IGBT 
(1000V/25A) avec une frequence de commutation maximale de 10 kHz, la tension continue 
a 1’ entree de l’onduleur est assuree a travers un redresseur triphase, enfin notre charge est 
une charge inductive-resistive (RL). On adopte les memes parametres que la simulation. 
Une breve description du materiel utilise est donnee en annexe A. 



Figure II. 7a Photographie du banc d’essai 



Figure II. 7b Schema synoptique du banc d’essai 


*•* U £• O <r c»V <t • ♦ St 8 /' * aj 

9 y 



~ — I 




Figure II. 7c Fenetre du logiciel control-desk 


Figure II. 7 Description generate de la plate-forme experimental 




T ension V ab ( V ) T ension V ab ( V ) 
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15 lO 15 20 


Rang des harmoniques 


Figure II. 8a Tension de sortie et son spectre d’harmoniques 



t(s) 

Figure II. 8b Courant de la charge 


Figure II. 8 Resultats experimentaux pour un onduleur monophase (k = 17) 



t(s) 


Rang des harmoniques 


Figure II. 9a Tension de sortie et son spectre d’harmoniques 



t(s) 

Figure II. 9b Courant de la charge 


Figure II. 9 Resultats experimentaux pour un onduleur monophase ( k = 21) 





Tension V as (V) Tension V as (V) 
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Figure II. 10a Tension de sortie et son spectre d’harmoniques 



O O.Ol 0.02 0.03 0.04 0.05 0.06 


t(s) 

Figure II. 10b Courant de la charge 

Figure II. 10 Resultats experimentaux pour un onduleur triphase (k = 17) 




t ( s ) Rang des harmoniques 

Figure II. 11a Tension de sortie et son spectre d’harmoniques 
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Figure II. lib Courant de la charge 


Figure II. 11 Resultats experimentaux pour un onduleur triphase (R — 21) 
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Les tensions de sortie et leurs spectres d’harmoniques sont presentees dans les figures II. 8a, 
II. 9a et II. 10a, II. 11a pour le cas du monophase et du triphase respectivement. 
D’apres l’analyse harmonique des tensions, nous constatons que les resultats experimentaux 
sont en conformite avec la simulation. Elle apparait une legere difference a cause de la 
tension d’entree de l’onduleur qui n’est pas parfaitement constante. On peut voir que la 
technique proposee elimine avec succes les harmoniques non desires tout en maintenant le 
fondamental a sa valeur desiree. Ces resultats montrent clairement la validite de la methode 
proposee. 

En plus, ces tensions sont appliquees sur une charge RL, les figures II. 8b, II. 9b, II. 10b 
et II. lib montrent le courant traversant la charge qui a une nature sinusoidal e avec des 
bruits a haute frequence. 

II.5- Extension de la strategic proposee aux resultats de la strategic MLI-ESH 

On va demontrer par la suite comment notre strategic proposee peut utiliser les valeurs des 
angles de commutation calculees pour la MLI-ESH. 

On prend comme exemple les valeurs des angles de commutation pour les quatre cas 
indiques dans le tableau 11.13 (pour un onduleur triphase avec m, = 0.70) : 



La symetrie par rapport au quart 
et a la demi-periode 

La symetrie par rapport 
a la demi-periode 


MLI-ESH 

Bipolaire 

MLI-ESH 

Unipolaire 

MLI-ESH 

Bipolaire 

MLI-ESH 

Unipolaire 

a i 

0.2269 

0.7503 

0.3057 

0.5221 

0,2 

0.2734 

0.8362 

0.4909 

1.6209 

«3 

1.1527 

0.9840 

0.6868 

1.8770 

a.4 

1.2937 

1.1596 

0.8405 

1.9441 

as 

1.4896 

1.2303 

1.0617 

2.1955 

a 6 

* 

* 

1.2085 

2.2389 

a 7 

* 

* 

2.4689 

2.6961 

«8 

* 

* 

2.6353 

2.8536 

«9 

* 

* 

2.8109 

3.0315 

aio 

* 

* 

3.0218 

3.1067 


Tableau 11.13 Valeurs numeriques des angles de commutation 


On utilise la meme expression du signal de la porteuse modifiee dans l’equation II. 11 pour 
la symetrie par rapport au quart et a la demi-periode. Pour la symetrie par rapport a la demi- 
periode, on va effectuer quelques modifications : 
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1/ Le parametre q, prend les valeurs suivantes : 

qr=2, q 2 = 3, q 3 =4, q 4 = 5, q 5 = 6, q 6 = 7, q 7 =8, q 8 = 9, q 9 =10, q 10 = 11 


2/ Le signal de la porteuse modifiee a une symetrie par rapport au quart et a la demi- 
periode, done pour etablir une forme qui a une symetrie par rapport a la demi-periode tout 
d’abord, on prend une demi periode de la porteuse et on inverse cette derniere pour creer la 
forme de la deuxieme demi periode de la porteuse. 

Pour les autres etapes de calcul du parametre A, elles restent inchangeables et les valeurs 
numeriques sont determinees et reportees dans le tableau 11.14. 



La symetrie par rapport au quart 
et a la demi-periode 

La symetrie par rapport 
a la demi-periode 


MLI-ESH 

Bipolaire 

MLI-ESH 

Unipolaire 

MLI-ESH 

Bipolaire 

MLI-ESH 

Unipolaire 

Az 

1.12341 

2.89997 

0.49301 

-0.96097 

a 2 

-0.51640 

3.04776 

0.60454 

0.21314 

a 3 

-0.75507 

2.45159 

-0.64883 

1.66698 

A 4 

0.14036 

1.15862 

-0.14113 

2.43758 

A 5 

0.27579 

0.42115 

0.45943 

2.23210 

Ag 

0 

0 

-0.09660 

1.45786 

A 7 

0 

0 

-0.18674 

0.69212 

A 8 

0 

0 

0.09733 

0.23082 

A 9 

0 

0 

0.03223 

0.04890 

A 10 

0 

0 

-0.03085 

0.00498 


Tableau 11.14 Valeurs numeriques de A ; 

Les figures II. 12a et 13a montrent bien la symetrie par rapport au quart et a la demi-periode 
de la porteuse modifiee d’une part, d’autre part les figures II. 14a et 15a montrent aussi la 
symetrie par rapport a la demi periode de la porteuse modifiee. En outre, si on compare les 
allures des tensions de sortie et leurs spectres d’harmoniques presentes dans les figures 
II. 12b, 13b, 14b et 15b avec les figures 1. 10b, lib, 12b et 13b respectivement, on ne 
constate aucune difference enregistree au niveau des amplitudes que ce soit pour le 
fondamental ou pour les harmoniques non elimines e’est-a-dire on obtient les memes 
performances que la strategic MLI-ESH. Nous constatons aussi que notre strategic proposee 
a un avantage considerable dans le but d’exploiter touts les resultats des travaux effectues 
sur la technique MLI-ESH. 




Tension V as ( V ) Tension V as (V) 
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t ( s ) 

Figure II. 12a Signal de la porteuse modifiee et le signal de reference 
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Figure II. 12b Tension de sortie et son spectre d’harmoniques 

Figure 11.12 Resultats obtenus par la MLI-ESH-TR avec des angles de la MLI-ESH-B 
pour une symetrie par rapport au quart et a la demi-periode 
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Figure II. 13a Signal de la porteuse modifiee et le signal de reference 
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Figure II. 13b Tension de sortie et son spectre d’harmoniques 


Figure 11.13 Resultats obtenus par la MLI-ESH-TR avec des angles de la MLI-ESH-U 
pour une symetrie par rapport au quart et a la demi-periode 


Tension V as (V) Tension V as (V) 
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t ( s ) 

Figure II. 14a Signal de la porteuse modifiee et le signal de reference 
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Figure II. 14b Tension de sortie et son spectre d’harmoniques 


Figure 11.14 Resultats obtenus par la MLI-ESH-TR avec des angles de la MLI-ESH-B 
pour une symetrie par rapport a la demi-periode 
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Figure II. 15a Signal de la porteuse modifiee et le signal de reference 
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Figure II. 15b Tension de sortie et son spectre d’harmoniques 


Figure 11.15 Resultats obtenus par la MLI-ESH-TR avec des angles de la MLI-ESH-U 
pour une symetrie par rapport a la demi-periode 
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Conclusion 

Le principe de la MLI-ESH presentee dans le chapitre precedent est de calculer les angles 
de commutation par la resolution d’un systeme d’ equations non lineaires pour eliminer 
certains harmoniques du spectre de la tension desiree. Les resultats de ces calculs sont 
stockes dans des memoires puis utilises au cours du fonctionnement selon les besoins. Un 
recent avancement dans le domaine des techniques de commande des systemes 
d'alimentation donne la naissance d’une strategie qui permet la generation des angles de 
commutation a partir d’une comparaison entre deux signaux, le premier c’est la reference et 
le deuxieme c’est la porteuse sinuso'rdale modifiee, Par consequent on n’a pas besoin 
d’utiliser les memoires. Cette nouvelle technique est appelee l’elimination selective 
d’harmoniques en temps reel (MLI-ESH-TR). 

Les resultats presenters dans ce chapitre pour un onduleur monophase et triphase demontrent 
clairement la generalisation de notre strategie du point de vue le nombre d’harmoniques a 
eliminer, cela est due a la formule de la porteuse sinuso'rdale modifiee adoptee, sous forme 
lineaire ce qui facilite le calcul des differents parametres. De plus, notre strategie a 
l’avantage d’etre applicable a tous les resultats de la strategie MLI-ESH c’est-a-dire, on 
peut exploiter les angles de commutation calcules pour la strategie MLI-ESH et de creer un 
signal de la porteuse sinusoidal modifiee, tous qa en gardant le meme spectre 
d’harmoniques de la tension desiree. Les resultats experimentaux sont tres satisfaisants et 
prouvent clairement la validite de la technique proposee pour 1’ onduleur monophase et 
triphase. 




Chapitre III 

Simplification 
et implementation de ia 
MLI par elimination 
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Introduction 

Plusieurs travaux de recherche ont ete consacres a la simplification et 
l’implementation de la MLI a elimination selective d’harmoniques (MLI-ESH). 
La commande en temps reel par cette technique exige le stockage de tous les angles de 
commutation obtenus dans des memoires, la capacite de ces dernieres est relative au 
nombre d’harmoniques a eliminer et aussi au pas d’echantillonnage de l’intervalle de 
variation de l’indice de modulation. Cette methode est connue sous le nom la methode des 
tableaux (Look up table : LUT) [Jia-96], [Alb-12]. 

D’autres methodes sont basees sur la recherche d’une relation mathematique de la 
variation de la trajectoire des angles en fonction de l’indice de modulation. La strategie 
MLI-ESH a echantillonnage regulier (MLI-ESH-ER) est definie par des equations 
mathematiques sur des intervalles bien determines suivant les allures des trajectoires 
des angles de commutation obtenus par la MLI-ESH qui varient entre 0° et 60° ou entre 
0° et 90°, c’est-a-dire cette technique est applicable pour une symetrie par rapport au quart 
et a la demi-periode bipolaire avec un nombre pair d’harmoniques a eliminer pour un 
onduleur monophase et triphase [Bow-95], [Bow-99], [Bow-01], [Bow-06]. Par contre la 
technique MLI-ESH linearisant (MLI-ESH-L) repose sur la linearisation des trajectoires des 
angles de commutation obtenus pour la MLI-ESH par des equations du premier ordre, avec 
le cas d’une symetrie par rapport au quart et a la demi-periode unipolaire et un nombre 
impair d’harmoniques a eliminer, pour un onduleur monophase [Azl-06]. Une MLI-SHE 
simplified (MLI-ESH-S) basee sur l’algorithme des surfaces egales (Equal Area Algorithm) 
[Che-97] dans le cas d’une symetrie par rapport au quart et a la demi-periode unipolaire 
pour un onduleur monophase. [Vil-0.5] propose une amelioration de cet algorithme basee 
sur le centre de gravite (Centroid PWM : CPWM). Une autre approche est basee sur l’une 
des techniques de l’intelligence artificielle, il s’agit des reseaux de neurones. Son 
implementation pratique demande l’utilisation de la carte dite field programmble gate array 
(FPGA) et le langage VHDL [And-94], [Woe-08], [Ima-11], 

Dans les techniques decrites ci-dessus les variables de sortie sont des angles, par 
contre dans la MLI par elimination selective d’harmoniques en temps reel (MLI-ESH-TR) 
developpee dans le chapitre precedent les angles sont transformes au parametre A, a travers 
le signal de la porteuse modified, done on cherche comment generer le parametre A, le long 
de l’intervalle de variation de l’indice de modulation. Pour cela nous proposons d’utiliser 
l’interpolation polynomiale de Lagrange [Mak-14.4], 

Dans ce chapitre, nous allons demontrer mathematiquement et d’une maniere claire 
les differentes simplifications (approximations) et implementations des trajectoires des 
angles de commutation calcules par la MLI-ESH. Des resultats de simulation 
et experimentales de ces methodes sont presented. 
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III.1- La strategic MLI-ESH avec la methode des tableaux (MLI-ESH-LUT) 


On prend comme exemple les resultats des angles de commutation obtenu dans le premier 
chapitre dans le cas de la symetrie par rapport au quart et a la demi-periode pour un 
onduleur monophase et triphase. L’intervalle de variation de l’indice de modulation ([0, 1]) 
est discretise avec un pas de 0.1. Les resultats sont presentes dans les 
tableaux III. la et III. lb. 


m i 

O.i 

02 

0-3 

04 

Os 

0.0 

0.2856 

0.5712 

0.8568 

1.1424 

1.4280 

0.1 

0.2894 

0.5631 

0.8671 

1.1291 

1.4420 

0.2 

0.2926 

0.5541 

0.8765 

1.1149 

1.4558 

0.3 

0.2952 

0.5444 

0.8848 

1.0999 

1.4695 

0.4 

0.2971 

0.5338 

0.8918 

1.0838 

1.4832 

0.5 

0.2982 

0.5224 

0.8971 

1.0662 

1.4967 

0.6 

0.2982 

0.5099 

0.9000 

1.0464 

1.5103 

0.7 

0.2970 

0.4961 

0.8995 

1.0234 

1.5239 

0.8 

0.2941 

0.4805 

0.8931 

0.9950 

1.5377 

0.9 

0.2885 

0.4620 

0.8755 

0.9562 

1.5517 

1.0 

0.2767 

0.4364 

0.8312 

0.8930 

1.5661 


Tableau III. la : Onduleur monophase 


rrii 

O.i 

02 

o 3 

04 

Os 

0.0 

0.3338 

0.3338 

1.0472 

1.3941 

1.3941 

0.1 

0.3246 

0.3299 

1.0623 

1.3808 

1.4091 

0.2 

0.3040 

0.3149 

1.0773 

1.3658 

1.4223 

0.3 

0.2860 

0.3029 

1.0922 

1.3511 

1.4357 

0.4 

0.2696 

0.2931 

1.1072 

1.3365 

1.4491 

0.5 

0.2544 

0.2851 

1.1223 

1.3221 

1.4626 

0.6 

0.2403 

0.2785 

1.1374 

1.3078 

1.4761 

0.7 

0.2269 

0.2734 

1.1527 

1.2937 

1.4896 

0.8 

0.2142 

0.2694 

1.1682 

1.2799 

1.5031 

0.9 

0.2022 

0.2666 

1.1843 

1.2666 

1.5165 

1.0 

0.1907 

0.2647 

1.2020 

1.2547 

1.5299 


Tableau III. lb : Onduleur triphase 


Tableau III. 1 Valeurs numeriques des angles de commutation 
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L’ implementation en temps reel de cette methode dans le logiciel Matlab/simulink exige 
1’ utilisation du bloc « Look up Table » pour Stocker les differents angles. Dans le cas ou la 
valeur de l’indice de modulation ne correspond pas aux valeurs dans le tableau, on utilise 
une methode d’interpolation lineaire pour les deux points qui limitent le point a chercher. 




Figure III. la m, = 0.15 Figure III. lb m, = 0.45 



Figure III. lc m, = 0.75 


Figure III. 1 Spectres d’harmoniques de la tension pour un onduleur monophase 



Rang des harmoniques Rang des harmoniques 

Figure III.2a m, = 0.15 Figure III. 2b m, = 0.45 
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Figure III.2c m, = 0.75 


Figure III. 2 Spectres d’harmoniques de la tension pour un onduleur triphase 

La figure III.l et III. 2 montre l’analyse spectrale de la tension de sortie (La tension continue 
a 1’ entree de 1’ onduleur est egale a 100 V), on remarque que les harmoniques 3, 5, 7 et 9 
(5, 7, 11 et 13) sont nuls, et le fondamental prend la valeur souhaitee. 

La figure III.3 montre un schema simplifie de cette strategie. 



Look up Table 
tableau (III. lb) 


CXi 


OCs 


Generation 
du signal 
MLI 


f 


Gi G 6 | I 

^ L 



Figure III.3a Onduleur monophase 


Figure III. 3b Onduleur triphase 


Figure III. 3 Schema simplifie de l’implementation de la MLI-ESH-LUT 


III.2- La strategie MLI-ESH a echantillonnage regulier (MLI-ESH-ER) 


Le developpement de cette technique s’appuit sur des remarques sur la variation des 
trajectoires des angles de commutation obtenus par la MLI-ESH. II est applicable pour les 
onduleurs monophase et triphase, de plus il est valable pour toute valeur de N impaire 
( N : nombre d’angles). Les figures ci-apres montrent un exemple des angles de 
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commutation presente dans [Bow-99], [Bow-01] obtenus par la MLI-ESH en utilisant 
l’equation (1.3) avec la contrainte (1.4) : 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure III. 4a N-7 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure III. 4b N = 9 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

m,- m ; 

Figure III.5a N = 5 Figure III.5b N = 7 

Figure III. 5 Trajectoire des angles de commutation en fonction de m, pour un 

onduleur triphase 


A partir des figures III.4 et III. 5 on peut extraire les caracteristiques suivantes [Bow-99], 
[Bow-01] : 


- Les trajectoires des angles de commutation sont approximativement lineaires dans 
l’intervalle de variation de l’indice de modulation [0, 0.95] et [0, 0.80] au cas d’un onduleur 
monophase et triphase respectivement ; 
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- Les trajectoires des angles de commutation qui ont un indice impair sont paralleles et ont 
des pontes negatives. Au contraire ceux qui ont un indice pair sont aussi paralleles mais 
avec des pentes positives ; 

- Pour la valeur de l’indice de modulation nulle, la difference entre deux angles successifs 
pour un onduleur monophase est donne par la relation suivante : 

T = 2,71 (III.l) 

2.1V +1 

- Pour la valeur de l’indice de modulation nulle, la difference entre deux paires des angles 
successives pour un onduleur triphase est donne par la relation suivante 

2 71 

T = — — — (III. 2) 

3. (IV +1) 

- II faut noter que ces caracteristiques ne sont pas applicables aux allures illustrees a la 
figure 1.6. 

La modelisation mathematique de cette strategie est definie par les equations suivantes 

[Bow-99], [Bow-01] : 

1/ Pour un onduleur monophase : 

** Avec un indice des angles impair c’est-a-dire z = 1, 3, 5, , N 

a, = i | ~ M, |sin(z | + q>,) - a| (IH.3) 

M i = <P,=^, A = 0.2.(l-cos(^)) (HI-4) 

** Avec un indice des angles pair done i = 2, 4, 6, , IV- 1 

= * | - M 2 |sin(z | + cp 2 ) + B j (HI- 5) 


n 


4 ’ 


B = 0.145.sin(^ ? ^)) 
1.4 


M 2 =0.343.m ; , cp y = 


(III-6) 
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2/ Pour un onduleur triphase : 

** Avec un indice des angles impair (z = 1, 3, 5, ,N) 

T M.T T 

a, = O' + 1) — sin((z' + 1) - 


+ <Pi) 


** Avec un indice des angles pair (z = 2, 4, 6, , N- 1) 


. T M.T . .. T 
«, = «- — — sm(z- + p 2 ) 


Les valeurs des constantes T, cp x , cp 2 sont presentees au tableau III.2. 


(HI- 7) 


(III.8) 


zrti 

T 

<Pi 

<Pi 

0 < Z7i ; < 0.8 

2 .71 

3. (IV +1) 

N.T 

2 

T 

4 

0.8 < zrij <1.27 

n .. ( rrii -0.8^| n 

(1 + cos — .— 

3. (IV +1) ^ 0.4 ) 2) 

N.T 

2 

T 

4 


Tableau III. 2 Valeurs numeriques des constantes T, (p x , (p 2 

A partir des equations III.l jusqu'a III.8, les resultats correspondants a la mise en oeuvre de 
la technique MLI-ESH-ER sont montres dans les figures III. 6 et III. 7. 



mi 

Figure III. 6a N = 7 


90 

^ 80 

O 

CD 

3 70 
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d 

d 50 

CD 
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d 

■q- 

s 30 

CO 

d 

a 1 20 

H 

^ 10 
0 











































































































































0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


mt 

Figure III. 6b N = 9 


Figure III. 6 Trajectoire des angles de commutation en fonction de m, pour un 

onduleur monophase 
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0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

m ( m,- 

Figure III. 7a N = 5 Figure III. 7b N = 7 

Figure III. 7 Trajectoire des angles de commutation en fonction de m, pour un 

onduleur triphase 

Les valeurs numeriques des angles de commutation montrees dans les figures III.4-III.7 
pour la MLI-ESH et la MLI-ESH-ER sont presentees sous forme des tableaux dans 
l'annexe B. 


Notons que Act,- c’est l’erreur entre les valeurs de 1’ angle calculees par la technique 
MLI-ESH et MLI-ESH-ER. 



Figure III.8a N -7 


m i 

Figure III. 8b N = 9 


Figure III. 8 Evolution de l’erreur des angles de commutation en fonction de m,- 
pour un onduleur monophase 
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0 0.2 0.4 0.6 0.8 1 

m i 

Figure III.9a N = 5 



Figure III.9b N = 7 


Aaj 

Aa 3 

Aa 4 

Aa 5 

— Aa 6 
Aa 7 


Figure III. 9 Evolution de l’erreur des angles de commutation en fonction de m ; 
pour un onduleur triphase 


A partir des figures III. 8a et III. 8b on remarque que pour m, e [0, 0.9] et m , e [0, 1] l’erreur 
varie dans l’intervalle [-3, 3] degres, par contre pour m, e [0.9, 1] l’erreur atteint sa valeur 
maximale de -6 degres, done on deduit que les approximations appliquees ne sont pas 
vraiment fiables surtout dans l’intervalle [0.9, 1]. On remarque aussi que l’erreur diminue 
avec P augmentation de N. 

Pour les resultats illustres dans les figures III.9a et III. 9b on remarque que l’erreur est douze 
fois plus faible que celle du cas d’un onduleur monophase pour m, compris entre 0 et 0.9. 

Pour voir l’effet de l’erreur des angles de commutation sur la forme du spectre 
d’harmoniques de la tension de sortie, on considered deux valeurs de l’indice de 
modulation. Les spectres sont presentes aux figures III. 10 et III. 11 (La tension continue a 
l’entree de l’onduleur est egale a 100 V) 

Avec N = 9 et m ; = 0.5 (m, = 1), le premier harmonique non nul est le 19 eme et l’amplitude 
du fondamental prend la valeur 50 V (100V). A partir des figures III. 10a et 10b, on constate 
que les harmoniques ciblees (3 eme jusqu'a 17 eme ) restent non nuls et presented une 
amplitude maximale de 30% par rapport a P amplitude du fondamental. 

La figure III. 10c montre que l'amplitude du fondamental n’est pas proportionnelle a l’indice 
de modulation et presente une difference maximale de 10% (Figure Ill.lOd). 
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Rang des harmoniques 

Figure III. 10a Spectre d’harmoniques 
(m, = 0.5) 




Rang des harmoniques 

Figure III. 10b Spectre d’harmoniques 

(m, - 1) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

m i 

Figure III. 10c Evolution du fondamental 
en fonction de m. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 .! 

m i 

Figure Ill.lOd Evolution de l’erreur 
en fonction de m. 


Figure III. 10 Performances de la strategie MLI-ESH-ER 
pour un onduleur monophase (IV = 9) 


Rang des harmoniques 

Rang des harmoniques 

Figure III. 11a Spectre d’harmoniques 

Figure III. lib Spectre d’harmoniques 

(m,= 0.5) 

(m, = 1) 
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0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

m i 

Figure III. 11c Evolution du fondamental 
en fonction de m, 



m i 

Figure III. lid Evolutions de l’erreur 
en fonction de m t 


Figure III. 11 Performances de la strategie MLI-ESH-ER 
pour un onduleur triphase ( N = 7) 


Avec N = 7 et m, = 0.5 (m, = 1), le premier harmonique non nul est le 23 eme et l’amplitude 
du fondamental prend la valeur 25 V (50V). Au contraire du cas precedant, pour l’onduleur 
triphase nous avons obtenu des meilleurs resultats du point de vue qualite du spectre 
d’harmoniques de la tension, tout qa est montre a la figure III. 11. 


Les equations III.9 et III. 10 et le schema de la figure III. 12 resument 1’ implementation en 
temps reel de la strategie MLI-ESH-ER pour un onduleur monophase ( N = 9) et triphase 
(N= 7): 


cq = 0.0331. cos(m r n/2) - 0.0203.171; + 0.1323 
a 2 = 0.0240.sin(m i .^'/2) + 0.0249.771; + 0.3307 
« 3 = 0.0331.cos(771;.^-/2) - 0.0452.771; + 0.4630 
a 4 = 0.0240.sin(T7i;.^-/2) - 0.0070.TT1; + 0.6614 
< a 5 = 0.0331.cos(tt1;.^/2) - 0.0652.777; + 0.7937 (III.9) 

a 6 = 0.0240.sin(T7i;.^-/2) - 0.0116.771; + 0.9921 
a 7 = 0.0331.cos(77i;.^-/2) - 0.0782.771,. +1.1244 
a 8 = 0.0240. sin(77i;.^-/2) - 0.0290.771; +1.3228 
« 9 = 0.0331. cos(77i;.^-/2) - 0.0827.171; +1.4551 
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0 < m i < 0.8 

a 1 = -0.1209.rn,. + 0.2618 
a 2 = -0.0421.rn,- +0.2618 
« 3 = -0.1298.m ; + 0.5236 
< a A = -0.0727.rn,. + 0.5236 (III. 10a) 
a 5 = -0.1298.m,. + 0.7854 
« 6 = -0.0984.m, + 0.7854 
a 7 = -0.1209.m, +1.0472 


0.8 < m . 

<1.27 


— (1 + cos 
24 

2 m i -0.8^ 

2 'N 

K 

l 0.4 y 

'2J 


m T 9 X 

a,=T ^sin(— ) 

2 2 

m t J . 5 J 

a 7 =T '■ — sm( ) 

2 2 4 

mj . 11J 

«3 = 21 


a 4 = 2.T - 


2 2 

m..T . .9 J. /"ttt i 

— — sin( ) (HI. 10b) 

2 4 

m,J . .13 J\ 


= 3 - r - -^- sin (— ) 

m l . 13.X 
«e = 3 - T 

m.J . 15 J 
a 7 = 4.T ! — sin( ) 




Figure III. 12a Onduleur monophase 


■ 

m, <= 0.8 

Systeme d’equations (III. 10a) 

■1 

m, > 0.8 

Systeme d’equations (III. 10b) 

OCi 

• 

• 

• 

• 

• 

• 

• 

► 

0.7 


Generation 
du signal 
MLI 

f 




G x G 6 | I 

^ L 



Figure III. 12b Onduleur triphase 


Figure III. 12 Schema simplifie de l’implementation de la MLI-ESH-ER 
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III.3- La strategic MLI-ESH linearisante (MLI-ESH-L) 

Le travail presente dans [Azl-06] est base sur le developpement des equations simples 
du premier ordre, derivees de la linearisation des trajectoires des angles de commutation 
pour un onduleur monophase. La methode de linearisation proposee est etudiee pour le cas 
d’une symetrie par rapport au quart et a la demi-periode unipolaire, avec un nombre 
pair des angles. Les figures III. 13 montrent les allures des angles de commutation presentes 
dans [Azl-06], ces angles sont calcules par l’equation (1.13) qui satisfait 
la contrainte (1.14) : 



Figure III. 13a IV = 10 



Figure III. 13b IV = 12 




Figure III. 13c N - 14 Figure III. 13d N = 16 


Figure III. 13 Trajectoire des angles de commutation en fonction de m, 
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Le developpement de cette methode consiste a exploiter les observations sur les allures de 
la figure III. 13 comme premiere etape, on peut citer : 

- Chaque trajectoire des angles de commutation est divisee en deux intervalles adequats, 
cette segmentation permet une representation lineaire de ces trajectoires. La premiere 
segmentation pour m, e [0, 0.85] et la deuxieme pour m, e [0.85, 1], 

- Les trajectoires des angles de commutation qui ont un indice impair sont paralleles et ont 
des pentes negatives. Au contraire que ceux qui ont un indice pair qui sont aussi paralleles 
mais avec des pentes positives ; 

- Pour la valeur de l’indice de modulation nulle, la difference entre deux paires des angles 
successives est donnee par la relation suivante : 

T = — — — (III.11) 

N +1 

La deuxieme etape est de presenter chaque trajectoire par une equation du premier ordre de 
la forme : 


a K =a.m i +b,K= [1 : N] (III. 12) 

Les tableaux III.3, III.4, III.5 et III.6 donnent les valeurs numeriques des coefficients a et b 
de l’equation (III. 12). Le calcul de ces coefficients pour chaque trajectoire se fait a partir de 
la fonction Curve Fitting dans le logiciel MATLAB. 


l’indice de l’angle 
K 

m i <= 0.85 

rrii >= 0.85 

a 

b 

a 

b 

1 

-2.9765 

16.4439 

-5.27 

18.33 

2 

1.1958 

16.5241 

-3.46 

20.33 

3 

-5.7250 

32.8992 

-10.82 

37.09 

4 

2.4780 

33.0196 

-7.53 

41.26 

5 

-7.9698 

49.3489 

-17.14 

56.92 

6 

3.9649 

49.4348 

-13.08 

63.63 

7 

-9.3564 

65.7756 

-25.02 

78.81 

8 

5.8213 

65.7183 

-21.10 

88.56 

9 

-9.4538 

82.0476 

-35.60 

104.2 

10 

7.9985 

81.8184 

8.57 

81.34 


Tableau III. 3 Valeurs numeriques des coefficients a et b pour N = 10 
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l’indice de Tangle 
K 

rrii <= 0.85 

mi >= 0.85 

a 

b 

a 

b 

1 

-2.1910 

13.9100 

-4.0340 

15.4300 

2 

0.8459 

13.9600 

-2.5370 

16.7300 

3 

-4.2530 

27.8300 

-8.2260 

31.1000 

4 

1.7310 

27.9100 

-5.3740 

33.7600 

5 

-6.0410 

41.7400 

-12.7900 

47.3100 

6 

2.7090 

41.8200 

-8.9100 

51.4500 

7 

-7.3860 

55.6400 

-18.0700 

64.5000 

8 

3.8570 

55.6700 

-13.7400 

70.4300 

9 

-8.0820 

69.4900 

-24.6300 

83.3800 

10 

5.2530 

69.4200 

-20.1900 

91.1200 

11 

-7.8990 

83.2600 

-32.9000 

104.6000 

12 

6.7920 

83.0800 

7.2610 

82.6800 


Tableau III.4 Valeurs numeriques des coefficients a et b pour N - 12 


l’indice de Tangle 
K 

m, <= 0.85 

m,- >= 0.85 

a 

b 

a 

b 

1 

-1.6830 

12.06 

-3.194 

13.3 

2 

0.6249 

12.09 

-1.933 

14.19 

3 

-3.2880 

24.11 

-6.476 

26.74 

4 

1.2720 

24.17 

-4.033 

28.54 

5 

-4.7280 

36.17 

-9.962 

40.49 

6 

1.9690 

36.24 

-6.509 

43.25 

7 

-5.9100 

48.21 

-13.82 

54.77 

8 

2.7540 

48.27 

-9.661 

58.61 

9 

-6.7270 

60.24 

-18.32 

69.91 

10 

3.6720 

60.24 

-13.88 

75.04 

11 

-7.0500 

72.23 

-23.84 

86.38 

12 

4.7550 

72.15 

-19.04 

92.52 

13 

-6.7670 

84.14 

-30.34 

104.3 

14 

5.9000 

84.00 

6.277 

83.69 


Tableau III. 5 Valeurs numeriques des coefficients a et b pour N = 14 
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l’indice de l’angle 
K 

m, <= 0.85 

rrii >= 0.85 

a 

b 

a 

b 

1 

-1.3250 

10.6300 

-2.5920 

11.6700 

2 

0.4891 

10.6500 

-1.5230 

12.3100 

3 

-2.6040 

21.2700 

-5.2370 

23.4300 

4 

0.9865 

21.3100 

-3.1470 

24.7100 

5 

-3.7800 

31.9000 

-8.0020 

35.3800 

6 

1.5120 

31.9500 

-4.9920 

37.3200 

7 

-4.7960 

42.5300 

-10.9800 

47.6400 

8 

2.0860 

42.5700 

-7.2200 

50.3000 

9 

-5.5880 

53.1400 

-14.3100 

60.3900 

10 

2.7340 

53.1700 

-10.0600 

63.8800 

11 

-6.0850 

63.7400 

-18.2000 

73.8900 

12 

3.4830 

63.7300 

-13.7500 

78.3100 

13 

-6.2100 

74.3000 

-22.8900 

88.4200 

14 

4.3400 

74.2400 

-17.2200 

93.3300 

15 

-5.9060 

84.8200 

-28.0600 

103.800 

16 

5.2200 

84.7100 

-5.5270 

84.4500 


Tableau III. 6 Valeurs numeriques des coefficients a et b pour IV = 16 

Dans le but de reduire le nombre d’ equations (car chaque angle ( a K ) est represente par une 
equation), on va presenter les constantes a et b chacune par une seule equation en fonction 
des indices des angles ( K ). II est preferable de presenter les coefficients a et b par des 
equations du deuxieme ordre et du premier ordre respectivement pour un indice de 
modulation inferieur a 0.85 et par des equations du deuxieme ordre pour un indice de 
modulation superieur a 0.85 respectivement [Azl-06] : 

a = d.K 2 + e.K + 1 (III. 13) 

b = h.K + q ou b = h.K 2 + q.K +v (III. 14) 

Les valeurs numeriques des coefficients d, e, l, h, q et v pour les differentes valeurs de N 
sont resumees dans les tableaux III.7 et III. 8 : 
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Indice des angles impair 

N 

d 

e 

l 

h 

q 

10 

0.1100464 

-1.9296291 

-1.0793990 

8.2042333 

8.2806310 

12 

0.0696106 

-1.4263959 

-0.7350828 

6.9376946 

7.0201072 

14 

0.0468942 

-1.0989697 

-0.5203273 

6.0099450 

6.0968365 

16 

0.0328576 

-0.8698174 

-0.3709166 

5.3015053 

5.3791336 

Indice des angles pair 

N 

d 

e 

/ 

h 

q 

10 

0.0385333 

0.3848488 

0.2867169 

8.1638842 

0.3170128 

12 

0.0220435 

0.2834859 

0.2097356 

6.9137648 

0.2492214 

14 

0.0139540 

0.2141497 

0.1628913 

5.9946794 

0.2086254 

16 

0.0092283 

0.1698528 

0.1363938 

5.2911726 

0.1703916 


Tableau III. 7 Valeurs numeriques des coefficients des equations III. 13 et III. 14 ( tit; <= 0.85) 


Indice des angles impair 

N 

d 

e 

/ 

h 

q 

V 

10 

-0.2074048 

-1.668593 

-3.5827694 

0.2747347 

7.9291278 

10.3635012 

12 

-0.1266383 

-1.3206454 

-2.8147927 

0.1707640 

6.8042432 

8.7479127 

14 

-0.0820399 

-1.0754506 

-2.2760433 

0.1124617 

5.9565599 

7.5497686 

16 

-0.05601594 

-0.89242431 

-1.8840180 

0.0778055 

5.2922308 

6.62994072 

Indice des angles pair 

N 

d 

e 

/ 

h 

q 

V 

10 

1.1608241 

-13.4052245 

22.0343197 

-0.9594558 

19.9788804 

-18.6347156 

12 

0.6240102 

-8.7400530 

16.0760857 

-0.5147764 

14.6469215 

-13.6030285 

14 

0.3616802 

-6.0147772 

2.2150204 

-0.2974016 

11.3341693 

-10.3344734 

16 

0.2225960 

-4.3258588 

9.5927857 

-0.1823041 

9.1424123 

-8.1110660 


Tableau III. 8 Valeurs numeriques des coefficients des equations III. 13 et III. 14 (m t >= 0.85) 


A partir des resultats indiques dans les tableaux III. 7 et III.8 et pour donner plus d’efficacite 
a cette methode, on peut presenter les coefficients d, e, l, h, q et v en fonction de la 
constante N par des equations du deuxieme ordre [Azl-06]. On obtient une equation 
generale qui englobe toutes les variables : 


** angles avec un indice impair : 


* rrti e[0, 0.85] 
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a K = [(0.0016499JV 2 - 0.0556129 .JV + 0.5007285 ).K 2 + (-0.0171300 .JV 2 + 
0.6207245 .JV - 6.4199924 ).K + (-0.0121815 .JV 2 + 0.4337315 .JV 
-4.1953445)].™,. +(0.0348811 .JV 2 -1.3887076.JV + 18.5972165).K + 
(0.0339263 .JV 2 -1.3634722.JV +18.5161382) 


(III. 15) 


*m,- e[0.85, 1.0] 


a K = [(-0.0034214 .JV 2 + 0.1138947.JV -1.0033324).K 2 + (-0.0103076 .JV 2 
+ 0.3966827 .JV - 4.6026316).K + (-0.0234969. JV 2 + 0.8926711.JV - 
10.1556696)].™, + (0.043321.JV 2 - 0.14509050. JV +1.2913232).K 2 + (III. 16) 

(0.028784 7.JV 2 -1.1863212.JV +16.9091770).K + (0.04348503.JV 2 - 
1.7505522.JV + 23.5135636) 

** angles avec un indice pair : 

* m; e[0, 0.85] 

a K = [(0.0007352.JV 2 - 0.0239169 .JV + 0.2039250 ).K 2 + (0.0035666 .JV 2 
- 0.1284484.JV +1.3123212)JC + (0.0031552 .JV 2 - 0.1069268 .JV + 

(III. 17) 

1.0399720)].™,. + (0. 0341632. JV 2 -1.3651064.JV +18.3928471).K + 

(0.0018473.JV 2 - 0.0720541.JV + 0.8515772) 

*mi e[0.85, 1.0] 

a K = [(0.0248581 .JV 2 - 0.8001616.JV + 6.6690673 ).K 2 + (-0.1860158.JV 2 + 

6.3345800 .JV - 58.1042658).K + (0.2084999. JV 2 - 7.4802823 .JV + 

75.94442303)].™,. + (-0.02059886 .JV 2 + 0.663011937.JV - 5.5234375).K 2 + (III. 18) 
(0.19626261 .JV 2 -6.8939358 .JV + 69.2470669).K + (-0.1755174.JV 2 + 

6.30542963 .JV - 64.1013648) 

On remarque une difference au niveau de quelques coefficients des equations III. 15-III. 18 
et ceux presentes dans [Azl-06], cette difference est due a 1 ’u t ilisation du logiciel Curve 
Expert par [Azl-06]. Afin de verifier les performances de cette strategic, un test est effectue 
pour un indice de modulation egal a 0.7, les figures III.14 et III. 15 montrent les spectres 
d’harmoniques pour differentes valeurs de m, et JV (La tension continue a l’entree de 
l’onduleur est egale a 100 V). 




Amplitude des harmoniques (V) Amplitude des harmoniques (V) Amplitude des harmoniques (V) 
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0 I ■ Lb ■— UJ 

1 5 10 15 20 25 30 35 

Rang des harmoniques 

Figure III.15d N~ 16 
Figure III. 15 Spectres d’harmoniques de la tension 

A partir des spectres d’harmoniques de la figure III. 14 et III. 15, on peut observer que la 
plupart des harmoniques cibles sont elimines. Le 27 eme harmonique pour IV = 14 apparait 
avec une amplitude de 5% par rapport au fondamental (Figure III. 15c). L’amplitude du 
fondamental a une bonne concordance avec la valeur de l’indice de modulation. Cette 
technique presente des avantages considerables du point de vue simplicity au niveau de la 
generation des angles de commutation et reduction du nombre d’equations. II faut noter que 
les tableaux III.3-III.6 presentent 52 equations dans chaque intervalle de l’indice de 
modulation. Avec cette technique ce nombre se reduit a 4 equations globales sans aucun 
risque sur la precision 

Pour une implementation en temps reel, la generation des angles de commutation se fait de 
la maniere indique sur le schema de la figure III. 16 avec le systeme d’equations III. 19. 


a. sot 


6 


.? 50 \- 


.9" 5 


10 15 20 2E 

Rang des harmoniques 

Figure III. 15c IV = 14 


E 

CO 4 


E 2 
< 


cq = -2.8923374.rn,. +16.4722966 
a 2 =1. 2083224. m ; +16.6319724 
« 3 = -5. 8671162. m ; + 32.8687976 
a, = 2.4366287.m,. +32.9481786 
cq = -7.9651789.m + 49.2652986 

,m,€[0, .85] (III.19a) 

a 6 = 3.9711432.m, + 49.2643848 
a 7 =- 9.1865257.m ; +65.6617996 
a 8 = 5.8118655.m, + 65.5805910 
a 9 = -9.5311566.m,. + 82.0583006 
a w = 7. 9587959. m ; +81.8967972 
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a l = -5.4517285.nl; +18.5546077 
a 2 = -0.1156491.m,. +17.4563147 
a 3 = -10.4370609.rn,. +36.5925033 
a 4 = -12.9966051.rn,. + 45.88580501 
a. = -17.0745966 .m. + 56.8194246 

, m, e[0.85, 1] (III. 19b) 

a 6 = -16.6514706.m, +66.6896622 
« 7 = -25.3643353.m, + 79.2353713 
a 8 = -ll.0802458.rn,. + 79.8678863 
« 9 = -35.3062773.m,. +103.8403438 
a 10 = 3.7170693.m, +85.4204774 




Figure III. 16 Schema simplifie de l’implementation de la MLI-ESH-L 

III.4- La strategic MLI-ESH simplifiee (MLI-ESH-S) 

Nous avons deja vu aux deux chapitres precedents que le calcul des angles de commutation 
ou les angles de perturbation se fait par la resolution des systemes d’ equations non lineaires. 
Plusieurs methodes mathematiques on ete proposees chacune par son niveau de complexite, 
dont les solutions sont exactes. La methode proposee par [Vil-05] est une amelioration de la 
methode proposee par [Che-97], cette methode permet la determination des angles de 
commutation, sans passer par la formulation et la resolution mathematique des systemes 
d’equations non lineaires. 

On doit noter que cette methode est valable pour une symetrie par rapport au quart et a la 
demi-periode unipolaire avec un onduleur monophase. 
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Le principe de cette strategie est decrit comme suit : on commence par la methode proposee 
par [Che-97], elle est reconnue par le nom : algorithme des aires egaux (Equal Area 
Algorithm). On prend l’exemple de la figure III. 17 pour une demi periode d’un signal 
sinusoidal (v ref = m i .sin(2.^. f .t), m, = 1 et f = 50Hz), la demi periode est divisee en cinq 

secteurs de largeurs egaux limites par deux angles ^,-et fr+i. 



Figure III. 17 Positionnement des differents angles 


L’aire limitee par la sinusorde et l’axe des abscisses en chaque secteur est donne par la 
formule suivante : 


fa+i 

aire i = m r J sin(o.t)d(Z»t = m i .(cos^ i - cos </> i+1 ),i =[1:3] (III. 20) 


A l’interieur de chaque secteur on a deux angles de commutation a, et a i+1 qui presentent 
l’extremite de l’impulsion de commande, done le principe de cette methode est de rendre 
l’aire de l’impulsion qui a une forme rectangulaire (avec une longueur unitaire) egale 
a l’aire calculee par l’equation (III. 20). 

Les angles de commutation dans chaque secteur pour un quart de periode sont calcules 
comme suit: 


_ A + aire i 

2 2 , z =[1:3] (HI-21) 

A +0M , aire i 


2 


2 


Amplitude des harmoniques (V) Amplitude des harmoniques (V) 
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La figure III. 18, presente les spectres d’harmoniques de la tension delivree par l’onduleur 
monophase controle par La strategic MLI-ESH-S. Avec cinq angles de commutation on 
elimine quatre harmoniques et on impose le fondamental a une valeur desiree, pour notre 
cas les rangs d’harmoniques concernes sont le 3 eme , 5 eme , 7 eme et le 9 eme . 

On remarque 1’ apparition des harmoniques de rang qui sont normalement nuls surtout 
l’harmonique du rang 9 qui a une amplitude considerable, qui diminue avec 1’ augmentation 
de l’indice de modulation. Aussi, pour l’amplitude du fondamental, on peut conclure que 
cette derniere est presque egale a la valeur desiree. D’autre part, la methode developpee par 
[Vil-05] est basee sur le centre de gravite de l’aire limitee par la sinuso'rde et l’axe des 
abscisses pour chaque secteur, elle est donnee par les relations suivantes : 


- _ sinfr/f ,) - tj) i+1 . cos(^. +1 ) - sin(^) + ^.cos(^) 
cos($)-cos(^ i+1 ) 


— cure. 

= x t L 

2 

, f =[1:3] 

— aire , 

X; + L 


(III. 22) 


(III. 23) 


Le tableau III. 10 donne les valeurs numeriques des angles pour quatre valeurs de l’indice de 
modulation : 


" ^ «/ (rad) 

a.i 

a.2 

a 3 

a.4 

a 5 







0.1 

0.4065 

0.4256 

0.9415 

0.9915 

1.5399 

0.4 

0.3779 

0.4543 

0.8665 

1.0665 

1.4472 

0.7 

0.3492 

0.4829 

0.7915 

1.1415 

1.3545 

1.0 

0.3206 

0.5165 

0.7165 

1.2165 

1.2618 


Tableau III. 10 Exemple des angles de commutation pour differentes valeurs de m, 

On constate une amelioration notable aux niveaux des spectres d’harmoniques presentes 
dans la figure III. 19 par rapport aux spectres de la figure III.18. Mais l’amplitude du 9 eme 
harmonique reste non nulle. 





Amplitude des harmoniques (V) Amplitude des harmoniques (V) 
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15 10 

Rang des harmoniques 


Figure III. 19a m,= 0.1 



15 10 

Rang des harmoniques 


Figure III. 19b 0.4 


70 
60 
50 
40 
30 
20 
10 
0 

15 10 15 10 

Rang des harmoniques Rang des harmoniques 

Figure III. 19c m, = 0.7 Figure III.19d m, - 1 



Figure III. 19 Spectres d’harmoniques de la tension 


Pour une implementation en temps reel, la generation des angles de commutation pour 
un quart de periode se fait de la maniere suivante avec le schema presente dans la 
figure III. 20 (pour cinq angles de commutation) : 

- calcul les aires dans chaque secteur en prenant m, = 1 : 

(a /re, = cos <j> r - cos tf> 2 ;aire 2 = cos <f) 2 - cos tf> 3 , aire 3 = cos tf> 3 - cos <j) A (III. 24) 


- calcul du centre de gravite de chaque aire : 

- _ sin(^ 2 ) - &.cos(#,) - sin(^) + ^.cos(^) — _ sin(^) - ^.cos(^) - sin(^ 2 ) + ^.cos(^) 

cos(^,)-cos(^ 2 ) ’ 2 cos(^ 2 )-cos((/S 3 ) 

< 

- _ sin(^ 4 ) - ^.cos(^) - sin(^) + &.cos(#,) 

cos(#,)-cos(^ 4 ) 


(III. 25) 
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— aire, 

<*i= x i 

— aire , 
a 2 =x 1 + m i .^±; 


a 3 = x 2 ~ m r 

a 4 = x 2 + m r 


aire 2 
2 

aire 2 
2 


— aire , 

«s = X 3~ m r^ A 


(III. 26) 



Figure III. 20 Schema simplifie de l’implementation de la MLI-ESH-S 

III.5- La strategic MLI-ESH temps reel (MLI-ESH-TR) 

L’inconvenient des methodes expliquees precedemment reside dans la specification 
de l’application, c’est-a-dire chaque technique a son propre modele mathematique qui 
ne peut pas etre applique a une autre technique, done une pour l’onduleur triphase et 1’ autre 
pour l’onduleur monophase soit bipolaire ou unipolaire et plus le nombre d’ angles a utiliser 
soit pair ou impair. Dans la MLI-ESH-TR developpee dans le chapitre precedent on n’a pas 
besoin de generer des angles de commutation mais de generer le signal de la porteuse 
modifiee afin de le comparer avec un signal de reference. Cette strategie est presque 
similaire a celle de la MLI triangulo-sinuso'idale. La generation du signal de la porteuse 
modifiee necessite la connaissance du parametre A,- le long de l’intervalle de variation de 
l’indice de modulation, tous qa pour un ajustement continu de l’amplitude de tension de 
sortie de l’onduleur [Mak-14.4], 

Les tableaux III. 11 et III. 12 presentent les resultats du parametre A, apres la resolution de 
l’equation 11.11, avec k = 17, pour un onduleur monophase et triphase respectivement : 
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rrii 

A] 

a 2 

a 3 

A4 

A5 

A.6 

A7 

Ag 

0.1 

-0.00346 

-0.00407 

-0.00239 

-0.00122 

0.00279 

-0.00183 

-0.04326 

6.39670 

0.2 

-0.00895 

0.00010 

-0.00109 

0.00004 

0.00146 

0.00163 

-0.08208 

3.18401 

0.3 

-0.01974 

0.00027 

-0.00066 

-0.00074 

-0.00041 

0.00486 

-0.12381 

2.12279 

0.4 

-0.03694 

0.00140 

-0.00013 

-0.00011 

-0.00074 

0.01213 

-0.16878 

1.58610 

0.5 

-0.06126 

0.00358 

-0.00054 

0.00025 

-0.00252 

0.02608 

-0.21882 

1.26103 

0.6 

-0.09598 

0.01024 

-0.00170 

0.00158 

-0.00782 

0.04933 

-0.27473 

1.03719 

0.7 

-0.14663 

0.02421 

-0.00534 

0.00540 

-0.02100 

0.09061 

-0.33864 

0.86456 

0.8 

-.022581 

0.05750 

-0.01963 

0.02142 

-0.05673 

0.16324 

-0.41055 

0.71197 

0.9 

-0.36446 

0.15204 

-0.08640 

0.09329 

-0.15859 

0.29082 

-0.47657 

0.54714 

1.0 

-0.76620 

0.62580 

-0.53956 

0.48833 

-0.45397 

0.41568 

-0.34644 

0.21114 


Tableau III. 11 Valeurs numeriques du parametre A, pour un onduleur monophase (k - 17) 


m, 

A; 

a 2 

a 3 

A4 

A5 

Ag 

A7 

Ag 

0.1 

-0.17836 

-0.11155 

0.11256 

-0.19441 

0.09543 

2.32702 

-4.64902 

2.32374 

0.2 

-0.17260 

-0.13418 

0.13610 

-0.21020 

0.10541 

1.16714 

-2.33312 

1.16600 

0.3 

-0.16846 

-0.15780 

0.16370 

-0.23013 

0.11690 

0.78168 

-1.56472 

0.78171 

0.4 

-0.16332 

-0.18477 

0.19604 

-0.25532 

0.13112 

0.59181 

-1.18667 

0.59221 

0.5 

-0.15745 

-0.21603 

0.23473 

-0.28701 

0.14974 

0.47794 

-0.96167 

0.47882 

0.6 

-0.15108 

-0.25271 

0.28253 

-0.32858 

0.17523 

0.39949 

-0.81034 

0.40197 

0.7 

-0.14563 

-0.29482 

0.34077 

-0.38113 

0.20737 

0.34070 

-0.69777 

0.34318 

0.8 

-0.14181 

-0.34023 

0.40574 

-0.44496 

0.24952 

0.29070 

-0.60268 

0.29070 

0.9 

-0.13748 

-0.37944 

0.45617 

-0.51028 

0.29895 

0.25028 

-0.51051 

0.23211 

1.0 

-0.11089 

-0.38998 

0.44263 

-0.54257 

0.32409 

0.24186 

-0.43184 

0.17088 


Tableau III. 12 Valeurs numeriques du parametre A, pour un onduleur triphase (k = 17) 


Pour les applications en temps reel, la valeur de l’indice de modulation est variable par la 
consequence de la variation de l’amplitude de la tension de reference. Pour atteindre cet 
objectif il faut chercher la relation mathematique qui lie A, et m,. Ceci peut etre realise en 
utilisant l'interpolation polynomial [Mak-14.4]. 

La methode d’interpolation la plus frequente est appelee l’interpolation polynomiale 
de Lagrange. Le principe fondamental de l’interpolation est de trouver une formule 
mathematique, souvent un polynome d’ordre n, dont le graphe passe par n+1 points 
(X(h (xu g(*i)), (x r , g(x r )) [Aln-08]: 

n+1 

G(x) = Y j g(x r )xH r (x),r=[l:n+l] 

r=l 


(III.27) 
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H r (x ) est donnee par : 


H r (x) 


pr ( x x s) 

7*7 ( X r ~ x s) 


(III. 28) 


Dans notre cas (x r , g(x,)) et C(x) correspondent a ( m, , A,) et P 7 y(m,). Suivant les tableaux 
III. 11 et III. 12 on a dix valeurs de l’indice de modulation. Par consequence l’ordre du 
polynome est n = 9 et z = 1, 2, 3, 4,..., 8. Apres quelques simplifications 

l’equation III. 27 devient : 


p A ( m i) = \ m f + b i, m i + + \ m i + b w, 


(III. 29) 


b t ,b 2 , ,b g ,b l0 : represente les coefficients du polynome qui sont reportes dans les 

tableaux III. 13 et III. 14. 



Pai 

PA3 

P A3 

Pa 4 

bn 

0.0282953 xl0 ! 

-0.0411466 xlO" 

0.0488640 xlO" 

-0.0279907 xlO" 

b 2 i 

-0.2778946 xl0 ! 

0.4016817x10" 

-0.4815674 xlO" 

0.2710865x10" 

b3i 

0.8218264 xlO a 

-1.1723942 xlO" 

1.4357534x10" 

-0.7777410 xlO" 

b 4i 

-1.1890344 xlO 3 

1.6692411x10" 

-2.1030589 xlO" 

1.0857222 xlO" 

bsz 

0.9801155 xl0 a 

-1.3506326 xlO" 

1.7638085 xlO" 

-0.8589452 xlO" 

b6i 

-0.4871285 xl0 ! 

0.6562796x10" 

-0.8958304 xlO" 

0.4063417x10" 

b7i 

0.1468623 xl0 a 

-0.1921582 xlO" 

0.2771129 xlO" 

-0.1152617x10" 

bsi 

-0.0261450 xlO 3 

0.0324862 xlO" 

-0.0502519 xlO" 

0.0188119x10" 

b<>i 

0.0024284 xl0 a 

-0.0028230 xlO" 

0.0048148 xlO" 

-0.0015868 xlO" 

blOi 

-0.0000917 xl0 ! 

0.0000918 xlO" 

-0.0001845 xlO" 

0.0000515x10" 


PA5 

P A 6 

PA7 

Pa8 

bn 

-0.0246250 xlO 3 

0.0190765x10" 

-0.0291519 xlO" 

-0.1690322 xlO" 

b 2 i 

0.2598678 xl0 ! 

-0.1881788 xlO" 

0.2857186x10" 

1.8633988 xlO" 

b3i 

-0.8819162 xl0 ! 

0.5612831x10" 

-0.8411482 xlO" 

-6.9074125 xlO" 

b 4i 

1.4637179 xl0 a 

-0.8184353 xlO" 

1.2102832 xlO" 

13.0206969 xlO" 

bsi 

-1.3710799 xl0 ! 

0.6825037x10" 

-0.9925312 xlO" 

-14.4355350 xlO" 

bsi 

0.7632368 xl0 ! 

-0.3452993 xlO" 

0.4907677x10" 

10.0012388 xlO" 

b7i 

-0.2533237 xlO 3 

0.1073355x10" 

-0.1474017 xlO" 

-4.4106076 xlO" 

bsi 

0.0482531 xl0 ! 

-0.0197494 xlO" 

0.0257902 xlO" 

1.2221492 xlO" 

b<ii 

-0.0047714 xlO ! 

0.0019588 xlO" 

-0.0027546 xlO" 

-0.2026699 xlO" 

blOi 

0.0001867x10" 

-0.0000790 xlO" 

0.0000814x10" 

0.0179846x10" 


Tableau III. 13 Coefficients des polynomes pour un onduleur monophase 
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Pai 

PA2 

PA3 

Pa 4 

bu 

0.7605852x10" 

-0.4434125x10" 

0.7590223x10" 

-0.1414783x10" 

k>2i 

-3.6784158x10" 

2.1890473x10" 

-3.6477323x10" 

0.7089485x10" 

b,3i 

7.6848172x10" 

-4.6199764x10" 

7.4728619x10" 

-1.5157721x10" 

b 4i 

-9.0544375x10" 

5.4707527x10" 

-8.5733643x10" 

1.8097155x10" 

bsi 

6.5945232x10" 

-3.9929737x10" 

6.0473550x10" 

-1.3230711x10" 

bei 

-3.0634078x10" 

1.8542665x10" 

-2.6981220x10" 

0.6098890x10" 

b7i 

0.9031156x10" 

-0.5466442x10" 

0.7555135x10" 

-0.1759058x10" 

bsi 

-0.1616215x10" 

0.0970453x10" 

-0.1247975x10" 

0.0300974x10" 

bgi 

0.0161713x10" 

-0.0114609x10" 

0.0131715x10" 

-0.0028854x10" 

blOi 

-0.0024389x10" 

-0.0005437x10" 

0.0005181x10" 

-0.0000803x10" 


PA5 

Pa6 

PA7 

Pa8 

bu 

0.1850208x10" 

-0.0825501xl0 4 

0.1228362xl0 4 

-0.0523748xl0 4 

b 2 i 

-0.9295588x10" 

0.4419015xl0 4 

-0.6783987xl0 4 

0.2960142xl0 4 

b3i 

1.9857895x10" 

-1.0295907xl0 4 

1.6331842xl0 4 

-0.7283171xl0 4 

b 4i 

-2.3590411x10" 

1.3719411xl0 4 

-2.2510420xl0 4 

1.0245357xl0 4 

bsi 

1.7094231x10" 

-1.1548834xl0 4 

1.9610538xl0 4 

-0.9096461xl0 4 

bsi 

-0.7789399x10" 

0.6397511xl0 4 

-1.1239734xl0 4 

0.5305325xl0 4 

b7i 

0.2218221x10" 

-0.2352903xl0 4 

0.4271566xl0 4 

-0.2048217xl0 4 

bsi 

-0.0376677x10" 

0.0564386xl0 4 

-0.1056156xl0 4 

0.0513457xl0 4 

b9i 

0.0035142x10" 

-0.0083816xl0 4 

0.0161073xl0 4 

-0.0079222xl0 4 

blOi 

-0.0000381x10" 

0.0006880xl0 4 

-0.0013517xl0 4 

0.0006710xl0 4 


Tableau III. 14 Coefficients des polynomes pour un onduleur triphase 


Nous connaissons bien que si on teste un polynome par les memes valeurs (l’indice de 
modulation) que nous utilisons dans l’interpolation on aura les meilleurs resultats, pour 
cela il est necessaire d’evaluer notre strategic par d’autre valeurs qui sont : m, = 0.15, 
0.45, 0.75. 


Les figures III.21 (III. 22), presentent les spectres d’harmoniques de la tension delivree 
par l’onduleur monophase (triphase) avec une tension continue a l’entree de ce dernier 
egale a 100V. 

Avec des valeurs de l’indice de modulation egales 0.15, 0.45 et 0.75 le fondamental prend 
les amplitudes suivantes 15V, 45V et 75V (7.5V, 22.5V et 37.5V). D’autre part avec 
k = 17 le premier harmonique non nul est le I7 eme (25 eme ). Du spectre d’harmoniques, 
on observe clairement, que notre strategic donne des resultats vraiment tres satisfaisants 
et demontrent la faisabilite de notre implementation en temps reel. 
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Q. 

E 

< 



1 5 10 15 20 25 

Rang des harmoniques 


Figure III.22c m, = 0.75 


Figure III. 22 Spectres d’harmoniques de la tension pour un onduleur triphase 
Un schema d’ implementation en temps reel est presente dans la figure III. 23a et III. 23b 


Equation (III. 29) 

•PaI P*A8 

Tableau (III. 13) 
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Figure III. 23a Onduleur monophase 


Figure III. 23b Onduleur triphase 


Figure III. 23 Schema simplifie de la strategie 


III.6- Resultats experimentaux 


Dans cette partie, nous presentons les resultats experimentaux des differentes 
implementations en temps reel de la strategie MLI-ESH. Les releves experimentaux choisis 
sont : le courant d’une phase et la tension et son spectre d’harmoniques. Ces essais ont ete 
effectues dans les memes conditions que les simulations. La figure III. 24 montre les 
differents dispositifs necessaires pour realiser cette experience. 
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Figure III. 24 Photo graphie du banc d’essai 

III.6.1- Resultats experimentaux de la MLI-ESH-LUT 




t(s) t(s) 

Figure III.25a Profil de l’indice de modulation et la tension de sortie 
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Figure III.25b Tension de sortie et son spectre d’harmoniques (m, = 0.15) 
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t ( s ) Rang des harmoniques 

Figure III.25c Tension de sortie et son spectre d’harmoniques ( m, = 0.45) 




t ( s ) Rang des harmoniques 

Figure III.25d Tension de sortie et son spectre d’harmoniques (m, = 0.75) 
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Figure III.25e Courant de la charge 

Figure III. 25 Resultats experimentaux de la MLI-ESF1-LUT monophase 
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T ension V ab ( V ) T ension V ab ( V ) 
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t(s) 

Figure III.27d Courant de la charge 


Figure III.27 Resultats experimentaux de la MLI-ESF1-ER monophase 





Tension V as (V) Tension V as ( V ) 
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Figure III. 28a Profil de l’indice de modulation et la tension de sortie 
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Figure III.28b Tension de sortie et son spectre d’harmoniques (m, = 0.5) 
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Figure III.28c Tension de sortie et son spectre d’harmoniques (m, = 1) 
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- 0 . 4 1 1 1 1 1 1 1 1 1 

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

t(s) 

Figure III.28d Courant de la charge 


Figure III.28 Resultats experimentaux de la MLI-ESF1-ER triphase 


III.6.3- Resultats experimentaux de la MLI-ESH-L 




t(s) t(s) 

Figure III. 29a Profil de l’indice de modulation et la tension de sortie 



Figure III.29b Tension de sortie et son spectre d’harmoniques (m, = 0.1) 
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0 


Figure III.28f Courant de la charge 
Figure III.29 Resultats experimentaux de la MLI-ESF1-L 


III.6.4- Resultats experimentaux de la MLI-ESH-S 
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Figure III.30b Tension de sortie et son spectre d’harmoniques (m, = 0.1) 
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Figure III.30c Tension de sortie et son spectre d’harmoniques (m, = 0.4) 
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Figure III.30d Tension de sortie et son spectre d’harmoniques (m, = 0.7) 
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Figure III.30e Tension de sortie et son spectre d’harmoniques (m, = 1) 
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0 0.1 


Figure III.30f Courant de la charge 
Figure III.30 Resultats experimentaux de la MLI-ESF1-S 


III.6.5- Resultats experimentaux de la MLI-ESH-TR 



Figure III. 31a Profil de l’indice de modulation et la tension de sortie 




t ( s ) Rang des harmoniques 

Figure III.31b Tension de sortie et son spectre d’harmoniques (m, = 0.15) 
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t(s) 


Rang des harmoniques 


Figure III.32c Tension de sortie et son spectre d’harmoniques (m, = 0.45) 
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Figure III.32d Tension de sortie et son spectre d’harmoniques (m, = 0.75) 
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Figure III.32e Courant de la charge 
Figure III. 32 Resultats experimentaux de la MLI-ESH-TR triphase 

Les resultats experimentaux sont presentes dans les figures III. 25 a III.32, ces resultats 
mettent en evidence la similarite existant entre les resultats simules et experimentaux. Une 
legere difference remarquable entre les resultats de simulation et experimentaux due a la 
forme de la tension d’entree de l’onduleur, car elle possede une forme ondulee. 

Conclusion 

Cinq methodes de l’implementation en temps reel de la MLI-ESH sont presentees dans 
ce chapitre. On peut distinguer entre ces methodes par deux choses : le processus 
mathematique pour laquelle on developpe la methode (simple ou complexe) et aussi par la 
qualite du spectre d’harmoniques de la tension delivree. 

Pour notre cas on utilise comme partie logiciel le Matlab/Simulink et materiel 
le dSPACE 1103. Dans ce cas on ne peut pas introduire le critere est ce que telle methode 
est simple a implementer ou non ? Car toute implementation avec ce materielle est 
simplifiee beaucoup si on compare avec 1 ’u t ilisation d’un micro-processeur. 
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On commence par la MLI-ESH-LUT qui est la plus simple, car il suffit d’utiliser 
le bloc « Look up Table » sous Matlab pour Stocker les differents angles de commutation 
calcules par la MLI-ESH. Les allures des spectres d’harmoniques sont satisfaisantes. 
Le developpement mathematique de la MLI-ESH-ER est tres complique et pour le cas d’un 
onduleur monophase on n’aura pas une elimination complete des harmoniques ciblees mais 
une minimisation. Par contre dans le cas d’un onduleur triphase les spectres montrent 
l’elimination desiree. La modelisation mathematique correspondante a la MLI-ESH-L est 
donnee par des equations simples (premier et deuxieme ordre) mais le processus du calcul 
et de la simplification est tres long. L’ analyse harmonique montre 1’ absence des 
harmoniques cibles, ce qui valide l’efficacite de cette technique. Les simplifications ne 
donnent pas toujours des bons resultats, comme dans le cas de la MLI-ESH-S. 
La formulation mathematique de ce dernier est basee sur quelques equations simples, ce qui 
repercute sur la qualite de la tension delivree car on trouve le 9 eme (le cas etudie) 
harmonique avec une amplitude non nulle. Les resultats obtenus par ces technique MLI- 
ESH-TR sont vraiment tres satisfaisants, point de vue l’elimination desiree d’harmoniques 
tout en assurant un bon control du fondamental. La formulation mathematique necessite le 
passage par l’interpolation polynomiale de Lagrange pour generer le parametre A,-. La 
complexity de cette methode est liee au pas de discretisation de l’intervalle de l’indice de 
modulation c’est-a-dire au degre du polynome, done pour notre cas avec un pas 0.1 on 
obtient des bon resultats. 

Les resultats des travaux dans [Jia-96], [Alb-12], [Bow-95], [Bow-99], [Bow-01], 
[Bow-06], [Azl-06], [Che-97] et [Vil-0.5] presentent un manque d’aspect temps reel, 
c’est-a-dire les allures des tensions ou des courants pour une seule valeur de l’indice de 
modulation. Par contre les resultats presentes dans ce chapitre bases sur une variation 
continue de l’indice de modulation. Les resultats experimentaux des cinq strategies, 
montrent 1’ applicability de nos algorithmes developpes pour une implementation 
temps reel. 
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Introduction 

Dans le chapitre precedent, nous avons presente seulement les algorithmes de 
^implementation en temps reel de la MLI-ESH en boucle ouverte avec une charge 
inductive. On dit en boucle ouverte car ils operent en temps reel mais ne necessitent en 
entree que la reference de la tension desiree qui est en relation avec l’indice de modulation. 
Par opposition aux techniques dites en boucle fermee qui exigent la connaissance d’autres 
informations. Maintenant, notre charge c’est une machine asynchrone triphasee (MAS), 
done on a besoin de mesurer les tensions et les courants de phase et aussi la vitesse. 

Au debut des annees 1970, Blaschke a presente le principe de la commande 
vectorielle permettant un decouplage entre le flux et le couple de la MAS. Ce principe 
consiste a orienter le flux suivant une direction determinee, d’oii l’appellation commande a 
flux oriente, avec une phase et une amplitude imposees a partir des references souhaitees, 
d’oii l’autre appellation commande vectorielle. Cette commande a fortement contribue au 
remplacement de la machine a courant continu, qui dominait dans les processus industriels 
travaillant a vitesse et a couple variables par des machines a courant alternatif. Le controle 
vectoriel transforme le modele de la MAS d’un systeme multivariable couple, en deux sous- 
systemes completement decouples. Ainsi, en controlant de maniere separee les deux 
variables principals qui sont le flux magnetique et le couple electromagnetique [Nai-99], 
[Bez-03], [Zha-08], [Rut-11], 

La commande vectorielle de la MAS necessite une alimentation alternative capable 
de delivrer une tension la plus sinuso'rdale possible a amplitude et a frequence variables. 
Les onduleurs de tension sont de nos jours tres largement utilises dans les systemes 
d’entrainements electriques, l’onduleur presente l’interface physique entre le dispositif de 
commande et la machine. Ainsi, les techniques de controle par MLI ont fait l’objet de 
recherches intensives pendant ces deux dernieres decennies. Un nombre important des 
methodes differentes de par leur concept et leur performance ont ete developpees. 
La triangulo-sinuso'rdale et la modulation vectorielle qui sont largement connues et utilisees 
mais l’inconvenient de ces techniques est les pertes par commutation. L’elimination 
selective d’harmoniques avec l’avantage de faible frequence de commutation permet la 
diminution des pertes par commutation, l’augmentation de la duree de fonctionnement des 
composants du convertisseur et aussi de reduire le dimensionnement du filtre utilise dans le 
systeme et par consequent son cout [Kha-04], [Bou-08], [Ima-11]. La commande 
vectorielle avec la technique de la MLI-ESH connue dans la litterature utilise la 
MLI-ESH-LUT, mais aucun detail sur la methode de generation des impulsions de 
commande des interrupteurs de puissance done les tensions de sortie de l’onduleur, de plus 
les allures de la vitesse et du couple ne sont pas presentees sauf l’allure de la tension et du 
courant dans le regime permanent [Xun-10], [Yon-11]. 
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Dans ce chapitre, nous presentons le modele mathematique triphase de la machine 
asynchrone et sa transformation dans le systeme biphase. Ensuite, nous passerons a la 
commande vectorielle avec des regulateurs PI flous, 1’ onduleur est pilote par la technique 
de la MLI-ESH-TR qui est caracterisee par sa constante du nombre de commutations sur la 
periode de la tension de reference dans tous les regimes de fonctionnement. Des resultats de 
simulation pour differents tests sont presentes. 

IV.l- Moderation de la MAS 

La synthese des lois de commande d’une machine electrique ou la simulation d’un 
ensemble convertisseur-machine-commande est intimement liee au modele mathematique 
decrivant le comportement dynamique et statique de la machine. La MAS triphasee 
comporte un stator fixe de trois enroulements identiques dephases de 2;r/3 et un rotor 
mobile autour de l’axe de symetrie de la machine. On distingue deux types des rotors : soit 
des rotors bobines soit des rotors a cage d’ecureuil. 

Pour notre etude, nous considerons les hypotheses suivantes [Bou-08], [Lou-10] : 

- Le circuit magnetique n’est pas sature, ce qui permet d’ exprimer les flux comme 
fonctions lineaires des courants. 

- Les pertes par hysteresis et courant de Loucault sont negligees. 

- Les forces magnetomotrices creees par chacune des phases des deux armatures sont 
a repartition sinuso'rdale d’ou resulte du fait que l’entrefer est constant, que les inductances 
mutuelles entre deux enroulements varient sinuso'rdalement en fonction de 1’ angle entre 
leurs axes magnetiques. 

La machine est representee a la figure IV.l par ces six enroulements; l’angle 0 repere l’axe 
de la phase rotorique R a par rapport a l’axe fixe de la phase statorique S a . 



Figure IV. 1 Representation des enroulements de la MAS 
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Les equations des tensions des phases statoriques et rotoriques s’ecrivent comme suit : 
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Les equations magnetiques sous forme matricielle sont donnees par les expressions 
suivantes : 
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La matrice des inductances mutuelles (matrice de couplage rotor-stator) s’ecrit : 
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(IV. 3) 
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L’ equation mecanique qui determine la dynamique de la machine est donnee par : 
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IV.2- Modele biphase de la MAS 

La transformation de la machine triphasee en une machine biphasee consiste a remplacer 
tout enroulement triphase [(S a ,S b ,S c ) ou (R a ,R b ,R C )] par deux enroulements, l’un sur 

l’axe direct (d) et l’autre sur l’axe en quadrature (q) (Figure IV. 2). La transformation qui 
traduit ce passage du systeme triphase (a,b,c) au systeme biphase (d,q) est dite de Park. 



Figure IV.2 Representation de la MAS dans le repere de Park 


Les matrices de passage directe \p{0)] et inverse [p(< 9)] 1 sont donnees par : 
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(IV. 5) 


et on a : [ ; X q; Xo] 1 = [P(0)U X a ; X b ; X c ]\ [ ; X b ; X c ] f = [P(d)\' .[ X d ; Aq ; X 0 f 

Le vecteur X peut etre le vecteur de tension, de courant, ou de flux. 
Pour les grandeurs statoriques 9 = 6 S et pour les grandeurs rotoriques 9 = 9 r [Bag-99] : 


0=6 + 0 , 


co „ = 


(O, =- 


d6 s _ 

dt 

d0, 

dt 

d0 


: la pulsation statorique. 


: la pulsation rotorique. 


co = — = pd : la pulsation mecanique. 


(IV.6) 
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A partir de ce qui est precede, on obtient les equations suivantes : 
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Avec : 

L s =l s -M s : Inductance cyclique statorique 
L r = l r -M r : Inductance cyclique rotorique 
3 

L m = —M rs : Inductance mutuelle cyclique entre le rotor et le stator 


(IV. 7) 


(IV.8) 


(IV.9) 


IV.3- Principe de la commande a flux rotorique oriente ou la commande vectorielle 

Nous avons choisi le cadre du controle vectoriel, qui est tres largement applique 
industriellement et est devenu un classique dans le domaine de 1’ electrotechnique. 
La commande vectorielle a fait preuve par ses nombreuses applications industrielles, de ces 
qualites de performances et de precision. La commande vectorielle n’est pas une structure 
figee et definitive mais plutot un cadre general qui offre d’autres voies d’ investigation dans 
lesquelles de nombreuses strategies de commande peuvent etre implantees : observation 
du flux, application de 1’ intelligence artificielle dans la regulation comme la logique floue 
[Nas-02], [Raj-07], 

Le controle vectoriel (champ oriente) consiste a choisir un systeme d’axes et un type de 
commande (directe ou indirecte) permettant de decoupler le couple et le flux analogue a une 
machine a courant continu (figure IV. 3) [Bag-99], [Nai-99] : 

- Dans la commande indirecte, l'angle de Park 0 S est calcule a partir de la pulsation 
statorique, elle-meme reconstitute a l'aide de la vitesse de la machine et de la pulsation 
rotorique co r . 

- Pour la commande directe, 0 S est calculee directement a l'aide des grandeurs mesurees ou 
estimees. 
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Dans les cas de commande avec un capteur de vitesse (ou via un codeur incremental) on 
choisi une commande indirecte par orientation du flux rotorique car elle est plus facile a 
implanter. 



Figure IV.3 Analogie entre la MCC et la MAS commandee par flux oriente 

On oriente le systeme d’axes (d-q) de maniere a ce que l’axe d soit en phase avec le flux 
rotorique, c'est-a-dire : 


(/>dr = <t>r 


< 


rtr=0 



(IV. 10) 


En remplaqant (j) dr et <j) qr , apres quelques simplifications le systeme d’equations (IV. 7) 
devient 


V* = R.4 +a.L.^+^^-a.L,.a..L 


dt L, dt 


s s qs 


V = R s iqs + V-k d ^k +( °s-J L< / , r+ V-k M s k 


(IV.ll) 


Les expressions du systeme d’equations (IV.ll) peuvent etre exploitees telles quelles pour 
realiser la commande vectorielle mais v ds et v qs influent a la fois sur z ds et i qs done sur le flux 
et le couple, il est done necessaire de realiser un decouplage afin de rendre les axes d et q 
completement independants. Ce couplage est supprime generalement par une methode 
classique, dite de compensation. On va definir deux nouvelles variables de commande 
L,. S| et V, telles que v dS[ n’agit que sur i ds et n’agit que surz qs . 


v dSl =R s .i ds +cr.L t 


^ qSl =Rs-k +a - L s 


di 


ds 


dt 

di 


qs 


dt 


(IV. 12) 
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Avec : a =1 



LL 



On prend en consideration la supposition que le module du flux rotorique ne varie que 
lentement par rapport aux courants. Le principe consiste a faire la regulation des courants 
en negligeant les termes de couplage, ces derniers seront rajoutes a la sortie des regulateurs 
selon le schema de la figure IV.4 [Bag-99], [Bou-08] : 


e c = ™ s LJ q5 


e c=^s L s i ds + ^ nL< />, 


(IV. 13) 


1 Vj- Vj. -l- 

* 

v ds 


t 

ids 

rr® tr 

< lqs 5 O) .[FdEi qSl tn 

* 

Vqs 

► 

e<if+ 

R s +oL s p 

1 

l qs 

TV \e q 

\ Regulation 


i 

R s +oL s p 

vlachine 



Figure IV.4 Decouplage par compensation 

Le schema de la figure IV. 5 montre les differentes boucles de regulation de la vitesse et des 
courants pour une commande vectorielle classique (Reg : regulateur). 




Figure IV.5 Schema de regulation de la commande vectorielle de la MAS 
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IV.4- Application d’un controleur flou dans une commande vectorielle 

Afin de pouvoir appliquer la technique de la logique floue a la commande vectorielle, nous 
presentons dans cette partie les concepts de la logique floue et nous exposons tous les 
aspects methodologiques necessaries a la comprehension de cette methode. Nous etudions 
ensuite la conception d'un controleur flou pour la commande vectorielle. 

IV.4.1- Principe et definition de la logique floue 

La logique floue a ete introduite par le professeur Lotfi Zadeh en 1965 comme une 
generalisation de la logique binaire. L’interet de la logique floue reside dans sa capacite a 
traiter l’imprecis, l’incertain et le vague. Elle est issue de la capacite de l’homme a decider 
d’agir de faqon pertinente malgre le flou des connaissances disponibles. En effet, la logique 
floue a ete introduite pour approcher le raisonnement humain a l’aide d’une representation 
adequate des connaissances. Ces dernieres annees, la logique floue a ete appliquee dans 
divers domaines. Cependant, la commande floue est le domaine qui suscite le plus de 
curiosite et s’avere etre un champ d’ application tres actif [Bag-99], [Boub-09]. 

IV.4.1. 1- Ensembles flous 

Un ensemble classique A de U (par exemple U = R N ) est defini par une fonction 
caracteristique notee p,\(x) telle que : 


Soit U une collection d’objets appelee univers de discours, un ensemble flou A dans U est 
caracterise par une fonction d’appartenance, notee Pa(x) qui appliquee a un element x de U, 
retourne un degre d’appartenance Pa(x) de x a A. La fonction d’appartenance d’un ensemble 
classique peut prendre seulement deux valeurs {0,1}. Un ensemble flou peut etre represente 
comme un ensemble de paires ordonnees {(x, Pa(*)}, comme le montre la figure IV. 6. 



(IV. 14) 






1 


o 


> 



X 

Figure IV.6a Logique classique 


X 


Figure IV.6b Logique floue 


Figure IV.6 Comparaison entre la logique classique et la logique floue 
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IV.4.1.2- Variable linguistique 

C’est une variable dont les valeurs ne sont pas des nombres, mais des mots ou des phrases 
exprimes en langage naturel. La raison pour laquelle on utilise cette representation, est que 
le caractere linguistique est moins specifique que le caractere numerique. Line variable 
linguistique est generalement representee par un triplet (x, T(x), U) dans lequel [Boub-09] : 

x est le nom de la variable linguistique (vitesse, erreur, position, chaleur, ) ; 

T(x) est 1’ ensemble des valeurs linguistiques qui sont utilisees pour caracteriser x ; 

U est l’univers de discours de la variable linguistique x ; 

Par exemple, si l’erreur est considered comme une variable linguistique definie sur l’univers 
de discours U=[- 1, 1], ses valeurs linguistiques peuvent etre definies comme suit : 
T(Erreur)= {Negatif Grand (NG), Negatif Petit (NP), Environ Zero (EZ), Positif Petit (PP), 
Positif Grand (PG)}. Ces symboles linguistiques peuvent etre considered comme des 
ensembles flous dont les fonctions d’appartenance sont representees sur la figure IV. 7. 



Figure IV. 7 Representation de la variable linguistique 

IV.4.1.3- Fonctions d’appartenance 

Afin de pouvoir traiter numeriquement les variables linguistiques, il faut les soumettre a une 
definition mathematique a base des fonctions d’appartenance qui montrent le degre de 
verification de ces variables. 

Dans ce contexte, on associe a chaque valeur de la variable linguistique une fonction 
d’appartenance p,\{x). Les fonctions d’appartenance peuvent theoriquement prendre 
n’importe quelle forme. Les plus souvent sont [Jin-02]: 
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La fonction triangulaire : definie par trois parametres a,b etc : 




0 x < a 
x-a 


b-a 

c-x 


a<x<b 


f 


ou par //(x) = max 


b<x<c 


min 

v V 


f x-a c-x^ ^ 


b-a c-b 


0 

J ) 


c-b 
0 x < c 


(IV. 15) 


La fonction gaussienne : definie par deux parametres oetm: 


ju(x) = exp 


r 7 \ 

(■ x-m ) 2 
2xr 2 


(IV. 16) 


La fonction trapezoidale : definie par quatre parametres a, b, c et d : 


M (x) = 


0 x < a 

x-a 


a<x<b 

b-a 

1 b < x < c ou par //(x) = max 
d -x 


r 


min 

V V 


x-a d -x 


b-a d-c 


,0 

J J 


b<x<c 


d-c 
0 d < x 


(IV. 17) 


La fonction sigmoide : definie par deux parametres a et c : 


m=- — — — - (iv. is) 

l + exp(-a(x-c) 




Figure IV.8a Fonction triangulaire Figure IV.8b Fonction trapezoidale 
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Figure IV.8c Fonction gaussienne Figure IV.8d Fonction sigmoi'de 

Figure IV.8 Differentes fonctions d’appartenance 

On notera que la forme rectangulaire correspond exactement a la logique classique, et que 
le triangle peut etre considere comme un cas particulier du trapeze, si les deux sommets 
correspondent aux valeurs a et c se coincident. 

IV.4.1.4- Implication floue 


L’implication floue est un operateur qui permet d’evaluer le degre de verite d’une regie de 
la forme : « Si x est A Alors y est B », a partir des valeurs de la premisse d’une part et de 
celle de la conclusion d’autre part. Ce degre de verite est evalue a partir des degres 
d’appartenance de x a A et de y a B comme suit : 


ju R (x,y) = imp(u A (x), y. B (y)) 
Les operateurs les plus utilises en commande floue sont : 
Implication de Mamdani : 


jU R (x, y) = min (ju A (x), ju B (y)) 


Implication de Larsen : 


JU R (x, y) = p A (x).p n (y) 


(IV. 19) 


(IV. 20) 


(IV. 21) 


IV.4.1.5- Raisonnement flou 

En commande, le raisonnement flou se base sur deux regies : a savoir, le modus ponens 
et le modus ponens generalise. En logique classique le modus ponens permet a partir de la 
regie « Si x est A alors y est B » et du fait « x est A » de conclure le fait « y est B ». 
Zadeh a etendu ce principe au cas flou, sous le nom modus ponens generalise pour 
permettre de raisonner lorsque les regies ou les faites sont connus d’une faqon imparfaite. 
Les deux modus sont resumes comme suit [Boub-09]: 
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Modus ponens Modus ponens generalise 
Premisse (fait) xestA x est A ’ 

Implication (regie) Si x est A alors y est B Six est A alors y est B 

Consequence (deduction) y est B y est B’ 

IV.4.2- Conception d’un controleur flou 

Le controleur flou est compose des quatre elements suivants : regies floues, fuzzification, 
moteur d’inference, et defuzzification. La figure IV. 9 montre le bloc de ce controleur flou 

[Jin-02], [Boub-09], [Piu-10] : 


Base de regies 
et de definitions 


Variables d’entree 
reelles 


Variables d’entree 
floues 


Interface de 



Interface de 

fuzzification 



defuzzification 


Variables de sorties 
reelles 


Moteur 

d’inference 


Variables de sortiel 
floues 


Figure IV.9 Schema general d’un controleur flou 


La base de regies floues : ou base de connaissances, contient les regies floues decrivant le 
comportement du systeme, elle permet de definir les fonctions d’appartenance. 
D’une maniere generale, la base de regies d’un systeme flou doit respecter les conditions de 
completude et de consistance afin d’ assurer le bon fonctionnement de ce dernier. Une base 
de regies d’un systeme flou est dite complete, si pour chaque vecteur d’entree il existe au 
moins une regie floue activee. Afin d’assurer cette propriete, les fonctions d’appartenance 
doivent couvrir tout l’espace des variables d’entree. Une base de regies d’un systeme flou 
est dite inconsistante, s’il existe deux regies floues ayant la meme premisse mais des 
conclusions differentes. La propriete de consistance permet d’eviter les contradictions dans 
une base de regies. 


La fuzzification : permet de transformer la grandeur physique d’entree en une grandeur 
floue (Numerique/Symbolique). On 1 ’u t ilise pour pouvoir activer les regies qui sont 
exprimees a l’aide de variables linguistiques associees a des ensembles flous. 
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Le moteur d’inference floue : une action sous forme symbolique, est decidee a l’aide des 
techniques de raisonnement flou en fonction des variables floues precedemment definies. 
Les trois moteurs d’inference floue les plus utilises sont : le moteur de Mamdani, celui de 
Sugeno et celui de Tsukumoto. 

La defuzzification : ce module traduit l’action floue issue de l’inference en une grandeur 
physique (Symbolique/Numerique) directement applicable au processus a commander. 
Parmi les methodes de defuzzification on a : la methode du maximum, la methode des 
maxima, la methode du centre de gravite, la methode des hauteurs ponderees et la methode 
des hauteurs ponderees modifiees. 

Le schema de base d’un controleur PI flou repose sur la structure d’un regulateur PI 
classique a la difference que nous allons retenir la forme incrementale. Cette derniere donne 
en sortie, non pas la grandeur de commande a appliquer au processus, mais plutot 
l’increment de cette grandeur, ainsi qu’ a ses entrees on a "e" et "de" qui representent 
respectivement, l’erreur et sa derivee. Les gains K e , K de , K A u sont des gains d’adaptation, ils 
permettent de changer la sensibilite du regulateur flou sans changer sa structure 
(Figure IV. 10) [Jin-02], [Lar-13] : 



e : est la difference entre la grandeur de reference (consigne) et la grandeur mesuree. 
La forme generale de la loi de commande adoptee est la suivante : 


U(t) = U(t-l) + G AU .AU(t) 


(IV.22) 


Avec : 

U(t), U(t-l): sont respectivement les commandes appliquees au systeme aux instants t et t-1. 
AU(t) : est la variation de la commande (l’increment a rajouter a la commande U(t)). 

Le regulateur de vitesse et de courant ont la meme structure, la seule difference est au 
niveau des valeurs des gains d’adaptation. 
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L’univers de discours des variables floues e et de est limite a l’intervalle [-1, 1], on a sept 
variables linguistiques qui sont (la terminologie est en anglais) : Negative Big (NB), 
Negative Medium (NM), Negative Small (NS), Zero (ZE), Positive Small (PS), Positive 
Medium (PM), Positive Big (PB). La forme triangulaire est utilisee pour presenter tous les 
variables linguistiques sauf la variable NB et PB sont presentes par la fonction trapezo'idale. 
Apres plusieurs essais en simulation nous pouvons etablir les plages suivantes des 
grandeurs d'entrees et de sortie (Figure IV. 11). 




-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 -1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 

Figure IV. 11a Erreur Figure IV. lib Derivee de l’erreur 

Figure IV. 11 Fonctions d’appartenance 

Le nombre total des regies floues est de 49, elles sont representees dans le tableau IV. 1, 
avec : Negative Very Big (NVB), Negative Big (NB), Negative Medium (NM), 
Negative Small (NS), Zero (ZR), Positive Small (PS), Positive Medium (PM), 
Positive Big (PB), Positive Very Big (PVB). 



Erreur (e) 

NB NM NS ZE PS PM PB 

Derivee de l’erreur (de) 

NB 

NVB ( Cl ) j NVB (c 2 ) | NVB (c 3 ) j NB (c 4 ) j NM (c 5 ) j ZE (c 6 ) j PVB (c 7 ) 

NM 

NVB (c 8 ) ! NVB ( C g) | NB (c 10 ) | NB (c u ) | NM (c 12 ) | ZE (c 13 ) | PVB (c 14 ) 

NS 

NVB (c 15 ) | NB (c 16 ) j NB (c 17 ) j NM (c 18 ) j PS (c 19 ) j PB (c 20 ) j PVB (c 21 ) 

ZE 

NVB (c 22 ) | NB (c 23 ) | NM (c 24 ) j ZE (c 25 ) | PM (c 26 ) | PB (c 27 ) | PVB (c 28 ) 

PS 

NVB (c 29 ) ! NB (c 30 ) ! NS (c 31 ) j PM (c 32 ) \ PB (c 33 ) j PB (c 34 ) j PVB (c 35 ) 

PM 

NVB (c 36 ) | ZE (c 37 ) | PM (c 38 ) | PVB (c 39 ) | PB (c 40 ) | PVB (c 41 ) ! PVB (c 42 ) 

PB 

— 

NVB (c 43 ) | ZE (c 44 ) | PM (c 45 ) ! PVB (c 46 ) \ PVB (c 47 ) \ PVB (c 48 ) \ PVB (c 49 ) 


Tableau IV. 1 Table des regies floues 

Le moteur d’inference utilise est de Takagi-Sugeno. Dans ce cas, des regies floues de type 
Sugeno sont utilisees. En effet, les conclusions des regies floues sont directement des 
constantes reelles (c,) (Tableau IV. 2). 
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Cl, C2, 
C3, Cg, 
C9, C15, 
C22, C29, 
C36, C43 

Cl , c 10 , 
Cib Ci6, 
Cl7, C23, 
C30 

c 5 , C 12 , 

Ci8, C24 

C 31 

C6, C13, 
C 25 , C37, 
C44 

C 19 

C26, C32, 
C38, C45 

C20, C27, 
C33, C34, 
C40 

C7, C14, 

C21, C28, 
C35, C39, 
C 41 , C42, 

C46, C47, 
C48, C49 

-1 

-0.3 

-0.1 

-0.02 

0 

0.02 

0.1 

0.3 

1 


Tableau IV.2 Valeurs numeriques des conclusions des regies floues 


L’implication floue est realisee par Poperateur min ou par le produit algebrique. 
La sortie finale est egale a la moyenne ponderee des conclusions des regies. Le produit 
algebrique pour Poperation d’implication floue et la moyenne ponderee pour la 
defuzzification [Boub-09], [Ara-12] : 


V, = Me, ( e ( f ))-/L e , ( de(t )) (IV.23) 


49 



AU(t) = ^ r - (IV.24) 



1=1 


La figure IV. 12 presente les tensions de reference v ds et v qs pour une commande vectorielle 
classique d’une MAS pour une vitesse de rotation donnee. L’onduleur est commande par 
MLI triangulo-sinuso'idale avec une frequence de la porteuse egale a 500 Hz. 

Notre choix de tel type de regulation soit par des regulateurs proportionnel integrale (PI) 
classique ou par des regulateurs PI flous (FLC : Fuzzy Logic Controller) est base sur 
Pobservation montree dans la figure IV. 12 : 



-50 >- 
0 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

t(s) 



t(s) 


Figure IV. 12a Regulateur PI classique 
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t(s) t(s) 

Figure IV. 12b Regulateur PI flou 


Figure IV. 12 Tensions de reference v ds et v qs 


IV.5- Commande vectorielle avec la MLI-ESH-TR 


La commande vectorielle avec un controle de l’onduleur par l’une des techniques 
de la MLI-ESH, exige le passage par l’indice de modulation (qui est l’image de la tension 
de reference) pour generer la tension desiree a la sorite de l’onduleur. On utilise aussi la 
valeur moyenne des tensions de reference v* et v qs dans notre commande vectorielle 
developpee de la MAS avec la technique MLI-ESH-TR. On remarque que, avec des 
regulateurs PI classiques, les tensions de reference v ds et v qs presentent des oscillations a 
grande amplitude, par contre avec des regulateurs PI flous on obtient des tensions avec le 
minimum d’ oscillations, pour ces raisons on choisi des regulateurs PI flous. 

On va modifier le schema de la figure IV. 5 pour developper notre schema de la commande 
vectorielle de la MAS avec la technique MLI-ESH-TR. II faut noter que le schema reste le 
meme jusqu'a la generation des tensions de reference v* et v qs . Pour elaborer ce schema, 
quatre etapes sont necessaries. 

Avant de presenter ces etapes, on passe par un petit rappel : 


V ds- 


•qs' 


^ sc =&s+Qc’ O c =artg 


V( v ds 2+v i 2 J-sin(6» sc ) 

(IV. 2 5) 

* o- * 
> 

V 

(IV.26) 


V V ds j 


Dans la premiere etape, on commence par le calcul de la valeur de l’indice de modulation 
(mi) a partir de la valeur moyenne des tensions de reference v* et v qs par l’equation 
suivante : 



(IV. 2 7) 
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Dans la seconde etape, on insere la valeur de m ; dans l’equation (III. 29) et avec le 
tableau (III. 14) on determine les valeurs du parametre A, ; 

La troisieme etape s’occupe de la construction du signal de la porteuse modifiee, on utilise 
les valeurs du parametre A, et l’angle 0 SC . L’expression du signal de la porteuse modifiee est 
donnee par l’equation (II. 11). 

Finalement, les signaux de commande de l’onduleur sont obtenus a partir d’une 
comparaison entre le signal de la porteuse modifiee et le signal de reference. 

Le schema de la figure IV. 13 montre d’une maniere simple et detaillee notre commande 
vectorielle de la MAS avec la technique MLI-ESH-TR. 



Figure IV. 13 Schema propose de la commande vectorielle de la MAS 


IV.6- Resultats de simulation 

Afin de verifier l'efficacite de la commande proposee, plusieurs simulations ont ete 
effectuees. La MAS utilisee dans la simulation est un moteur a deux poles 
avec une puissance de 3kW, la tension continue a l’entree de l’onduleur est egale a 514 V. 
Les parametres electriques et mecaniques sont : R s = 1.31 F2; R, = 2.013 Q; 
L s = L r = 0.184 H; L m = 0.167 H; J = 0.0097 kg.m 2 . Une methode d’identification de ces 
parametres est presentee dans l’annexe C. 








Tension v, s (V) T em (N.m) 


Chapitre IV: Commande vectorielle de la MAS alimentee par un onduleur pilote par la MLI 
a elimination selective d’harmoniques en temps reel 


112 



t (S) 

Figure IV. 14a Couple electro magnetique 
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Figure IV. 14b Vitesse de rotation 
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Figure IV.14g Tension de reference v ds 
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Figure IV. 14f Angle 6 C 
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Figure IV.14h Tension de reference v qs 


Figure IV. 14 Evolution des differentes grandeurs mecaniques et electriques de la MAS 


Amplitude 
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Figure IV. 15 Signal de la porteuse modifiee 
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t (s) 

Figure IV. 16 Signal de commande du transistor T n 



t (s) 

Figure IV. 17a Q r * = 1300 Trs/min 




t (s) 

Figure IV. 17c O, = 1200 Trs/min 


Figure IV. 17 Signal de la porteuse modifiee pour chaque vitesse de reference 


Tension v, s (V) Tension v, s (V) Amplitude 
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Figure IV. 19 Tensions de sortie de l’onduleur et leurs spectres d’harmoniques 




Figure IV.20a Couple electromagnetique Figure IV.20b Vitesse de rotation 
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Figure IV.20c Tension de sortie de l’onduleur 
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Figure IV.20d Courant statorique 
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Figure IV.20e Indice de modulation 
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Figure IV.20f Angle Q c 
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Figure IV.20g Tension de reference v ds 



t (s) 

Figure IV.20h Tension de reference v qs 


Figure IV.20 Evolution des differentes grandeurs mecaniques et electriques de la MAS 




Figure IV.21 Signal de la porteuse modifiee 




Figure IV.22 Tension de sortie de l’onduleur et son spectre d’harmoniques 



Figure IV.23a Couple electromagnetique 


Figure IV.23b Vitesse de rotation 
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t(s) 

Figure IV.23c Tension de sortie de 1’ onduleur 



t(s) 

Figure IV.23d Courant statorique 




Figure IV.23e Indice de modulation Figure IV.23f Angle 6 C 



Figure IV.23g Tension de reference v ds 
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Figure IV.23h Tension de reference v qs 


Figure IV.23 Evolution des differentes grandeurs mecaniques et electriques de la MAS 
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Figure IV.24a Q r * = 1200 Trs/min (C r = 0 N.m) 
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Figure IV.24b Q, = 1200 Trs/min (C r = 10 N.m) 


Figure IV.24 Tensions de sortie de l’onduleur et leurs spectres d’harmoniques 


Notre objectif de ces differents tests, est d’assurer une bonne regulation de la vitesse 
et de confirmer que Pelimination desiree des harmoniques dans la tension de sortie de 
l’onduleur est garantie. 


La vitesse de la machine suit parfaitement sa reference, soit par un changement de consigne 
(1300, 800, 1200 Trs/min), inversion du sens de rotation (1000, -1000 Trs/min) ou par une 
application d’un couple resistant (1200 Trs/min et C r = 10 N.m) (Figure IV. 14b, 
20b, 23b) 


Le couple electromagnetique apres chaque changement de vitesse oscille autour de zero en 
tant qu’il n’y a pas de couple resistant. Le couple electromagnetique oscille autour de la 
valeur 10 N.m qui est la valeur du couple resistant applique a la machine (Figure IV. 14a, 
20a, 23a). 

Les tensions a la sortie de l’onduleur et les courants statoriques sont presentes dans les 
figures IV. 14c, 20c, 23c et figures IV.14d, 20d, 23d respectivement. 

L’indice de modulation prend presque la meme forme que la tension de reference v qs , 
d’autre part l’influence de la tension de reference v ds est negligeable (Figure IV.14e, 14g, 
14h, 20e, 20g, 20h, 23e, 23g, 23h). L’angle 0 C presente des valeurs importantes au 
demarrage et lors de l’inversion du sens de rotation, par contre en regime permanent, sa 
valeur est faible de quelques degres mais ils sont necessaires comme les montrent 
les figures IV.14f, 20f, 23f. 

Le signal de la porteuse modifiee possede une symetrie par rapport au quart et a la 
demi-periode et aussi, l’apparition des fluctuations lors de l’inversion du sens de rotation 
(Figure IV.15, 17, 21). La comparaison entre le signal de la porteuse modifiee et le signal 
de reference produit le signal de commande des interrupteurs de l’onduleur. On remarque 
bien que dans une periode de la tension de reference, on obtient le meme nombre 
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d’impulsions de commande qui est 17 done 33 changements d’etat (commutations) de 
l’interrupteur de puissance, pour tous les tests effectues (Figure IV. 16, 18). 

L’ analyse harmonique des tensions de sortie montre des tres bons resultats du point de vue 
F elimination desiree e’est-a-dire tous les harmoniques cibles sont nuls. L’amplitude du 
fondamental prend la valeur souhaitee qui est en relation avec l’indice de modulation 
(Figure IV. 19, 22, 24). 

Le signal de commande du transistor Tn de la figure IV.25 est obtenu pour une commande 
vectorielle avec la MLI triangulo-sinusoidale (MLI-TS) a partir d’une comparaison entre le 
signal de reference et une porteuse triangulaire. La frequence de la porteuse triangulaire est 
fixe et egale a 737 Hz, cette valeur permet de generer 17 impulsions de commande du 
transistor Tn pendant une periode de la tension de reference avec une vitesse de rotation 
egale a 1300Trs/min (Figure IV. 25a). 

En comparant les resultats de la figure IV. 18 et IV.25, on constate que pour une vitesse de 
rotation egale a 800 Trs/min (1200 Trs/min) le nombre d’impulsions augmente de 17 a 28 
(17 a 19) ce qui donne une augmentation du nombre de commutations de 66% (12%). 
La commande vectorielle avec la MLI-ESH-TR donne des performances comparables 
a la MLI-TS, tout en gardant le meme nombre de commutations lors des changements des 
consignes, e’est-a-dire vitesse ou couple de reference. 
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Figure IV.25a Q, = 1300 Trs/min 
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Figure IV.25b II, = 800 Trs/min 
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Figure IV.25 Signal de commande du transistor Tn pour chaque vitesse de reference 
pour une commande vectorielle avec la MLI triangulo-sinusoidale 
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Conclusion 

Dans ce chapitre, nous avons presente le modele mathematique de la MAS a cage en tenant 
compte des hypotheses simplificatrices, 1 'u t ilisation de la transformation de Park a simplifie 
le modele de la MAS. Ensuite, nous avons presente notre schema de la commande 
vectorielle de la MAS alimentee par l'onduleur commande par la technique de la 
MLI-ESH-TR. 

Pour pouvoir juger notre schema propose, les trois tests essentiels sont effectues comme 
dans la commande vectorielle classique. On peut conclure, dans tous ces resultats que 
le decouplage est maintenu dans tout le regime de fonctionnement de la MAS, de plus 
la generation avec succes du parametre A,- par la methode d’ interpolation polynomiale avec 
la variation de la frequence et de l’indice de modulation en meme temps. 

La comparaison entre la MLI-ESH-TR et la MLI-TS dans une commande vectorielle de la 
MAS, montre des performances notables de la MLI-ESH-TR en terme de sa constante du 
nombre de commutations dans les differents tests appliques. 

Finalement, nous avons montre d’une maniere claire notre methode developpee qui est 
differente des travaux presentes dans la litterature. 
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Cette these a pour objectif d’apporter une nouvelle contribution dans le cadre des 
etudes menees sur la technique d’elimination selective d’harmoniques pour des onduleurs 
monophase et triphase. Ceci est accompli par le developpement d’une nouvelle technique 
optimisee d’elimination selective d’harmoniques similaire a la MLI triangulo-sinuso'idale, 
afin de simplifier la generation des signaux de commande des onduleurs. 

Nous avons commence par l’etude des differentes strategies MLI-ESH qui sont 
caracterisees par la forme et par la symetrie. Nous avons vu que la determination des angles 
de commutation (instants de commutation) revient a la resolution d’un systeme d’equations 
non lineaires a fonction trigonometrique. Le degre de complexity augmente avec le nombre 
d’harmoniques a eliminer et avec la symetrie appliquee. On remarque que pour une 
symetrie par rapport a la demi-periode le nombre d’angles de commutation a chercher 
devient le double par rapport au quart et a la demi-periode pour le meme nombre 
d’harmoniques a eliminer. L’algorithme des essaims particulaires est exploite pour resoudre 
les systemes d’equations obtenus. La forme unipolaire donne une ties bonne valeur du taux 
de distorsion harmonique car cette forme est proche de la forme sinuso'rdale. 
Son implementation est un peut compliquee qa est due a l’independance des signaux de 
commande des interrupteurs du meme bras c’est-a-dire qu’ils ne sont pas complementaires. 

L’implementation en temps reel de cette MLI-ESH est generalement mise en oeuvre 
par la methode des tableaux, qui exige pratiquement 1 ’u t ilisation des memoires 
electroniques pour Stocker les angles de commutation. La generation des instants de 
commutation par la MLI triangulo-sinuso'idale est en ligne (temps reel) c’est-a-dire par la 
consequence de la comparaison entre le signal de reference (sinuso'rdale) et le signal de la 
porteuse (triangulaire), pour cette raison le developpement d’une MLI-ESH temps reel est 
necessaire. Dans notre MLI-ESH-TR, au lieu d’une porteuse triangulaire, on utilise une 
porteuse sinuso'rdale modifiee. Cette derniere est une somme des constants multiplies par 
des fonctions sinusordales. En effet, pour chaque valeur du parametre k, un ensemble des 
angles de perturbation a ete calcule pour chaque pas de l’indice de modulation. 
On a aussi demontre 1’ applicability de la strategic proposee aux resultats calcules par 
la MLI-ESH, c’est-a-dire les angles de commutation calcules au chapitre I. 
Les resultats obtenus soit de simulation ou experimentaux pour un onduleur monophase et 
triphase sont en accord et satisfaisants. Ainsi ils permettent de prouver la validite de notre 
technique. 

Le vrai concept de l’implementation en temps reel de la MLI-ESH est d’obtenir une 
tension souhaitee a la sortie de 1’ onduleur avec 1’ elimination des harmoniques desires, mais 
d’une maniere continue. Si on suppose que notre parametre d’entree est l’indice de 
modulation done on obtient une tension souhaitee sur toute la plage de variation de l’indice 
de modulation. Pour cet objectif et pour completer les travaux presentes dans la litterature, 
nous avons demontre par simulation et experimentalement la procedure a suivre pour 
assurer une implementation en temps reel de la MLI-ESH. 
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En premier temps, notre charge connectee a l’onduleur est une charge inductive, 
malgre les bons resultats obtenus mais on cherche d’ exploiter la technique de la 
MLI-ESH-TR avec une charge dynamique en vue de sa commande. Done la charge utilisee 
e’est une machine asynchrone triphasee. L’etape de la modelisation de la machine utilisee 
est indispensable, nous avons presente les equations electriques, magnetiques et mecaniques 
dans le repere triphase et grace a la transformation biphasee de Park on diminue la 
complexite du modele. Nous avons elabore une commande vectorielle a flux rotorique 
oriente de la machine asynchrone alimentee par un onduleur pilote par la MLI-ESH-TR. 
A cet effet, nous avons presente une nouvelle structure de la commande vectorielle detaillee 
avec des controleurs PI flous qui est different du schema de la commande vectorielle 
classique. La partie principale dans ce nouveau schema est la methode de generation des 
signaux de commande de l’onduleur, done la methode appliquee consiste a calculer 
la valeur de l’indice de modulation a partir des tensions de reference directe et en 
quadrature. Les resultats obtenus pour cette nouvelle structure sont satisfaisants et montrent 
la faisabilite de la methode. 

Finalement, nous recommandons la poursuite des recherches sur la technique de la 
MLI-ESH-TR avec les entrainements electriques, ou plusieurs ameliorations peuvent etre 
apportees a ce travail, a savoir : 

• La simplification de la generation du parametre A,- par une autre methode qui permet 
de reduire les degres des polynomes ; 

• L'extrapolation de la strategic de la MLI-ESH-TR dans les onduleurs multi- 
niveaux ; 

• L’implementation experimentale de notre commande vectorielle proposee. 
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Description generate du banc d’essai 

Ce banc d’essai a ete conqu dans le but de pouvoir tester les differentes strategies de 
MLI-ESH et aussi les lois de commande. L’avantage du banc est la possibility de compiler 
les schemas de simulation Matlab/Simulink pour une simulation temps-reel, c’est a dire de 
permettre la generation des signaux de commande de l’onduleur a partir d’un schema 
Simulink sous Matlab realisee sur l’ordinateur avec la carte dSPACE 1103. Le banc d’essai 
est compose d’un onduleur, d’une carte de commande numerique dSPACE 1103, 
d’une carte d’interface, d’une carte des capteurs de courant et de tension et enfin une charge 
RL. La Figure A.l illustre les differents dispositifs du banc d’essai. 



Alimentation de 15 V pour les circuits 
de I’ onduleur 


Carte d’interface 


Figure A.l Photographie d’une partie du banc d’essai 
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A.l Convertisseur 

II est constitue d’un redresseur de tension a base de diodes de puissance (1000V/25A), d’un 
onduleur de tension compose de six transistors IGBT (1000V/25A) qui sont commutes a 
une frequence maximale de 10 kHz et aussi d’un hacheur abaisseur. Chaque IGBT est 
equipe d’une diode de roue libre (1000V/20A) mise en antiparallele (Figure A.2). 



Figure A.2 Photographie du convertisseur 

A.2 Carte d’interface dSPACE-convertisseur 

Cette carte permet d’amplifier les signaux MLI generes par la carte dSPACE pour attaquer 
les IGBT du convertisseur. L’ amplitude des signaux generes par la carte dSPACE varie 
entre 0-10V, mais les IGBT du convertisseur exigent une amplitude de 0-15V. 
Pour ces raisons on a realise une carte d’amplification et d’inversion a base des circuits 
integres SN7407 (Hex buffers/drivers with open-collector high-voltage outputs) et SN7406 
(Hex inverter buffers/drivers with open-collector high-voltage outputs) (Figure A.3). 
On genere seulement les signaux de commande des IGBT des bras du haut pour le cas d’un 
onduleur triphase, done on genere trois signaux par la carte dSPACE. L’ amplification 
et l’inversement (complementaire) des signaux est assure par notre carte. L’avantage de 
cette carte permet de reduire le nombre des DAC (digital to analog chanal) utilises ce qui 
permet de diminuer le pas de simulation. 
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Figure A.3a Vue externe 
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Figure A.3b Vue interne 
Figure A.3 Photographie de la carte d’interface 


A.3 Capteur de tension 

La mesure des tensions est faite a l’aide d’un capteur de tension a effet Hall, de type 
LV25-P. II delivre a sa sortie de mesure un courant proportionnel au courant dans son 
entree avec un rapport de (2500/1000). Une resistance R e a l’entree du capteur permet de 
convertir la tension a mesurer en un courant. Le courant mesure a sa sortie passe a travers 
une resistance shunt de mesure R m . Dans notre application, nous avons fixe la valeur 
de R e a 47KX2 et R m a 470Q (Figure A.4). 
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Figure A.4 Photographie de la carte des capteurs de tension 

A.4 Capteur de courant 

La mesure des courants est faite a l’aide d’un capteur de courant a effet Hall, de type 
LA25-P. II delivre a sa sortie de mesure un courant proportionnel au courant reel avec un 
rapport de (1/1000). Le courant mesure a sa sortie passe a travers une resistance shunt de 
mesure R m . Dans notre application, nous avons fixe la valeur de R m a 270F2 (Figure A.5). 



Figure A.5 Photographie de la carte des capteurs de courant 
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Valeurs numeriques des angles de commutation pour la MLI-ESH et MLI-ESH-ER 


m t 

oq 

a 2 

a 3 

a 4 

a 5 

«6 

a 7 

0.0 

12.0000 

24.0000 

36.0000 

48.0000 

60.0000 

72.0000 

84.0000 

0.1 

11.8671 

24.2285 

35.6281 

48.4254 

59.4625 

72.5589 

83.3986 

0.2 

11.7185 

24.4244 

35.2179 

48.8079 

58.8875 

73.0944 

82.7854 

0.3 

11.5545 

24.5842 

34.7685 

49.1428 

58.2702 

73. 6055 

82.1561 

0.4 

11.3752 

24.7029 

34.2769 

49.4210 

57.6019 

74.0879 

81.5038 

0.5 

11.1801 

24.7733 

33.7377 

49.6277 

56.8684 

74.5330 

80.8175 

0.6 

10.9674 

24.7844 

33.1418 

49.7373 

56.0453 

74.9229 

80.0768 

0.7 

10.7338 

24.7190 

32.4735 

49.7040 

55.0894 

75.2176 

79.2395 

0.8 

10.4718 

24.5455 

31.7023 

49.4367 

53.9154 

75.3104 

78.1970 

0.9 

10.1605 

24.1884 

30.7518 

48.7115 

52.3134 

74.7850 

76.5340 

1.0 

9.5963 

23.0082 

28.9343 

45.8733 

48.7076 

68.4929 

69.2369 


Tableau B.l Valeurs numeriques des angles de commutation pour N -7 (MLI-ESH monophase) 




a 2 

a 3 

a 4 

«5 

a 6 

a 7 

a 8 

a g 

0.0 

9.4737 

18.9474 

28.4211 

37.8947 

47.3684 

56.8421 

66.3158 

75.7895 

85.2632 

0.1 

9.3904 

19.0908 

28.1818 

38.1693 

47.0033 

57.2234 

65.8695 

76.2409 

84.7880 

0.2 

9.2967 

19.2124 

27.9152 

38.4098 

46.6045 

57.5730 

65.3968 

76.6772 

84.3054 

0.3 

9.1928 

19.3103 

27.6209 

38.6126 

46.1693 

57.8878 

64.8936 

77.0976 

83.8124 

0.4 

9.0785 

19.3817 

27.2972 

38.7720 

45.6931 

58.1609 

64.3528 

77.4997 

83.3047 

0.5 

8.9533 

19.4224 

26.9412 

38.8785 

45.1676 

58.3805 

63.7621 

77.8783 

82.7752 

0.6 

8.8160 

19.4261 

26.5474 

38.9172 

44.5790 

58.5251 

63.1003 

78.2223 

82.2112 

0.7 

8.6641 

19.3830 

26.1063 

38.8627 

43.9042 

58.5535 

62.3272 

78.5067 

81.5858 

0.8 

8.4924 

19.2744 

25.5990 

38.6677 

43.0985 

58.3760 

61.3564 

78.6639 

80.8291 

0.9 

8.2873 

19.0564 

24.9790 

38.2209 

42.0572 

57.7513 

59.9580 

78.4195 

79.6667 

1.0 

7.9994 

18.5593 

24.0838 

37.1191 

40.4379 

55.7600 

57.3270 

75.2765 

75.7258 


Tableau B.2 Valeurs numeriques des angles de commutation pour N = 9 (MLI-ESH monophase) 


m t 

a 4 

a 2 

a 3 

a 4 

a 5 

0.0 

20.0000 

20.0000 

40.0000 

40.0000 

60.0000 

0.1 

19.1215 

20.4537 

39.0881 

40.7230 

59.1299 

0.2 

18.2316 

20.9053 

38.1603 

41.4458 

58.2504 

0.3 

17.3289 

21.3507 

37.2133 

42.1671 

57.3592 

0.4 

16.4118 

21.7843 

3 6.2426 

42.8846 

56.4533 

0.5 

15.4779 

22.1986 

35.2418 

43.5950 

55.5281 

0.6 

14.5242 

22.5826 

34.2010 

44.2928 

54.5766 

0.7 

13.5462 

22.9191 

33.1049 

44.9674 

53.5871 

0.8 

12.5371 

23.1789 

31.9273 

45.5983 

52.5370 

0.9 

11.4855 

23.3086 

30.6199 

46.1367 

51.3753 

1.0 

10.3669 

23.1920 

29.0769 

46.4319 

49.9495 


Tableau B.3 Valeurs numeriques des angles de commutation pour N - 5 (MLI-ESH triphase) 
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mi 

a.! 

a 2 

a 3 

a 4 

«5 

a 6 

a 7 

0.0 

15.0000 

15.0000 

30.0000 

30.0000 

45.0000 

45.0000 

60.0000 

0.1 

14.3500 

15.2713 

29.3239 

30.4525 

44.3176 

45.5767 

59.3482 

0.2 

13.6932 

15.5402 

28.6360 

30.9003 

43.6244 

46.1523 

58.6910 

0.3 

13.0286 

15.8036 

27.9342 

31.3409 

42.9181 

46.7261 

58.0272 

0.4 

12.3550 

16.0577 

27.2158 

31.7708 

42.1955 

47.2971 

57.3550 

0.5 

11.6710 

16.2975 

26. 4764 

32.1851 

41.4515 

47.8637 

56.6716 

0.6 

10.9747 

16.5158 

25.7103 

32.5761 

40.6789 

48.4227 

55.9735 

0.7 

10.2628 

16.7020 

24.9078 

32.9311 

39.8651 

48.9689 

55.2535 

0.8 

9.5311 

16.8390 

24.0541 

33.2270 

38.9879 

49.4909 

54.4993 

0.9 

8.7718 

16.8970 

23.1216 

33.4169 

38.0024 

49.9623 

53.6817 

1.0 

7.9691 

16.8151 

22.0541 

33.3850 

36.7988 

50.2985 

52.7135 


Tableau B.4 Valeurs numeriques des angles de commutation pour N -7 (MLI-ESH triphase) 


m, 

oq 

a 2 

a 3 

a 4 

«5 

«6 

a 7 

0.0 

12.0000 

24.0000 

36.0000 

48.0000 

60.0000 

72.0000 

84.0000 

0.1 

11.7850 

24.5347 

35.5690 

48.3656 

59.4223 

72.2003 

83.3705 

0.2 

11.5117 

25.0500 

35.0796 

48.7119 

58.7863 

72.3812 

82.6825 

0.3 

11.1822 

25.5274 

34.5340 

49.0202 

58.0940 

72.5243 

81.9384 

0.4 

10.8000 

25.9504 

33.9357 

49.2742 

57.3491 

72.6129 

81.1416 

0.5 

10.3700 

26.3052 

33.2897 

49.4600 

56.5564 

72.6334 

80.2971 

0.6 

9.8982 

26.5814 

32.6018 

49.5671 

55.7219 

72.5752 

79.4107 

0.7 

9.3917 

26.7725 

31.8792 

49.5892 

54.8527 

72.4320 

78.4896 

0.8 

8.8584 

26.8764 

31.1298 

49.5240 

53.9566 

72.2015 

77.5416 

0.9 

8.3068 

26.8952 

30.3621 

49.3738 

53.0423 

71.8859 

76.5754 

1.0 

7.7459 

26.8354 

29.5852 

49.1450 

52.1187 

71.4918 

75.6000 


Tableau B.5 Valeurs numeriques des angles de commutation pour N -7 (MLI-ESH-ER monophase) 


m, 

a 4 

a 2 

a 3 

a 4 

«5 

a 6 

a 7 

«s 

«9 

0.0 

9.4737 

18.9474 

28.4211 

37.8947 

47.3684 

56.8421 

66.3158 

75.7895 

85.2632 

0.1 

9.3341 

19.3958 

28.1386 

38.2406 

46.9713 

57.0811 

65.8444 

75.9288 

84.7661 

0.2 

9.1484 

19. 8288 

27.8102 

38.5711 

46.5281 

57.3048 

65.3270 

76.0528 

84.2231 

0. 3 

8.9183 

20.2320 

27.4373 

38.8718 

46.0405 

57.4985 

64.7652 

76.1469 

83.6356 

0.4 

8.6467 

20.5922 

27.0229 

39.1295 

45.5113 

57.6494 

64.1619 

76.1981 

83.0066 

0.5 

8.3373 

20.8986 

26.5707 

39.3334 

44.9444 

57.7463 

63.5207 

76.1954 

82.3398 

0.6 

7.9950 

21.1429 

26.0855 

39.4751 

44.3446 

57.7812 

62.8466 

76.1307 

81.6400 

0.7 

7.6252 

21.3201 

25.5729 

39.5498 

43.7173 

57.7490 

62.1451 

75.9988 

80.9128 

0.8 

7.2342 

21.4283 

25.0392 

39.5555 

43.0688 

57.6479 

61.4224 

75.7980 

80.1645 

0.9 

6.8288 

21.4695 

24.4910 

39.4941 

42.4059 

57.4796 

60.6853 

75.5301 

79.4017 

1.0 

6.4161 

21.4485 

23.9355 

39.3707 

41.7357 

57.2493 

59.9409 

75.2001 

78.6316 


Tableau B.6 Valeurs numeriques des angles de commutation pour N = 9 (MLI-ESH-ER monophase) 
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m, 

a.: 

a 2 

«3 

a 4 

a 5 

0.0 

20.0000 

20.0000 

40.0000 

40.0000 

60.0000 

0.1 

19.0603 

20.42 26 

39.0000 

40.7071 

59.0603 

0.2 

18.1206 

20.8452 

38.0000 

41.4142 

58.1206 

0.3 

17.1809 

21.2679 

37.0000 

42.1213 

57.1809 

0.4 

16.2412 

21.6905 

36.0000 

42.8284 

56.2412 

0.5 

15.3015 

22.1131 

35.0000 

43.5355 

55.3015 

0.6 

14.3618 

22.5357 

34.0000 

44.2426 

54.3618 

0.7 

13.4222 

22.9583 

33.0000 

44.9497 

53.4222 

0.8 

12.4825 

23.3809 

32.0000 

45.6569 

52.4825 

0.9 

11.5240 

23.7377 

30.9254 

46.2403 

51.4046 

1.0 

10.3622 

23.1588 

28.8754 

45.0732 

48.4610 


Tableau B.7 Valeurs numeriques des angles de commutation pour N = 5 (MLI-ESH-ER triphase) 


m, 

oq 

a 2 

a 3 

a 4 

a 5 

a 6 

a 7 

0.0 

15.0000 

15.0000 

30.0000 

30.0000 

45.0000 

45.0000 

60.0000 

0.1 

14.3071 

15.2411 

29.2564 

30.4167 

44.2564 

45.5639 

59.3071 

0.2 

13.6142 

15.4822 

28.5128 

30.8334 

43.5128 

46.1278 

58.6142 

0. 3 

12.9213 

15.7232 

27.7692 

31.2500 

42.7692 

46.6916 

57.9213 

0.4 

12.228 4 

15.9643 

27.0257 

31.6667 

42.0257 

47.2555 

57.2284 

0.5 

11.5355 

16.2054 

26.2821 

32.0834 

41.2821 

47.8194 

56.5355 

0.6 

10.8425 

16.4465 

25.5385 

32.5001 

40.5385 

48.3833 

55.8425 

0.7 

10.1496 

16.6876 

24.7949 

32.9167 

39.7949 

48.9472 

55.1496 

0.8 

9.4567 

16.9286 

24.0513 

33.3334 

39.0513 

49.5110 

54.4567 

0.9 

8.7500 

17.1234 

23.2547 

33.6609 

38.2120 

49.9453 

53.6223 

1.0 

7.8942 

16.6629 

21.7677 

32.7309 

36.0885 

48.5546 

50.8635 


Tableau B.8 Valeurs numeriques des angles de commutation pour N = 7 (MLI-ESH-ER triphase) 
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Identification des parametres de la machine asynchrone [Bou-08] 

II existe de nombreuses methodes d’ identification des parametres de la MAS, parmi ces 
methodes, la methode classique qui est basee sur le schema equivalent en regime permanent 
avec les trois essais suivants : a courant continu, a vide et a rotor bloque. 

Les parametres necessaries pour la simulation du fonctionnement de la MAS se divisent en 
parametres electriques (R r , R s , L r , L s , L m ) et mecaniques (J, Kf). Les essais sont effectues sur 
une machine asynchrone de 3 kW (Figure C.l). 



Figure C.l Photographie de la MAS, frein a poudre et generatrice tachymetrique 

Comme les essais experimentaux seront effectues en regime sinusoidal, nous faisons appel 
au modele de la figure C.2 representant le schema monophase equivalent ramene au stator. 



Rf 

3 


Figure C.2 Schema monophase de la MAS 
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Dans ce schema, N r et R r represented: respectivement l’inductance de fuite totale et la 
resistance rotorique ramenees au stator. L’ expression de N r est donnee par : 



C.l Mesure de la resistance statorique 


(C.l) 


La mesure de la resistance statorique est effectuee en alimentant la machine en continu, 
en realisant le montage de la figure C.3. 



La valeur de la resistance est obtenue simplement par la relation : 



(C.2) 


E (V) 

2.62 

5.30 

7.90 

10.60 

Is (A) 

1 

2 

3 

4 


Tableau C.l Differentes valeurs de la tension et du courant 
La valeur moyenne de la resistance R s egale a 1.3 ID 

C.2 Essai a vide 

Dans le fonctionnement a vide le glissement est faible done on peut negliger les pertes 
Joules et ferromagnetiques rotorique. Le moteur asynchrone absorbe une puissance P so qui 
correspond a la somme des pertes mecaniques, ferromagnetiques et Joules au stator : 

Pso=Pjs+Pf+P m eca (C-3) 


Pjs =3 r s i; 


(C.4) 
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Ou les pertes mecaniques sont la somme des pertes par frottements (proportionnelles 
a la vitesse) et des pertes par ventilation (proportionnelles au carre de la vitesse). Les pertes 
ferromagnetiques sont les pertes par courants de Foucault et par Hysteresis 
(proportionnelles au carre de la tension V s , de l’induction B et de la frequence). 

Par la variation de la tension d’ alimentation V s on mesure le courant I s la puissance P so 
on gardera la vitesse constante. La puissance reactive Q so est deduite par la relation 
suivante : 


L’inductance cyclique L s (L r ) est donnee par la relation : 

r _ 3V * 2 

Qso®s 

Le tableau C.2 regroupe les differentes mesures effectuees. 


(C.5) 


(C.6) 


Sv s (V) 

Iso (A) 

Pso( W) 

Qso (VAR) 

Pjs (W) 

Pf + P meco (W) 

3V s 2 (V) 

400 

4 

240 

2760 

63.36 

176.64 

160000 

370 

3.2 

170 

2044 

40.55 

129.45 

136900 

340 

2.7 

150 

1583 

28.86 

121.14 

115600 

300 

2.2 

125 

1137 

19.16 

105.84 

90000 

270 

1.9 

100 

883 

14.29 

85.71 

72900 


Tableau C.2 Mesures effectuees pour differentes valeurs de la tension 


L s = L r = 0.184 H 


C.3 Essai en court-circuit 


Le rotor etant bloque (g = 1), la tension statorique doit etre ajustee de telle sorte que le 
courant statorique soit proche du courant nominal. Dans ces conditions, on peut negliger le 
courant de la branche magnetisante ce qui conduit aux formules suivantes : 


R = 


31 


— R 


(C.7) 


L 


m 


-N r +^N?+4L 2 S 


2 


(C.8) 
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V SC (V) I sc (A) P SC (W) Q sc (VAR) 

56.4 5 250 808 

R r = 2.013 n 
L m = 0.167 H 

C.4 Determination des parametres mecaniques : J, Kf 

La determination du coefficient de frottement visqueux et du moment d’inertie est basee sur 
la mesure des pertes mecaniques, lorsque la machine tourne a une vitesse donnee et sur le 
releve de la courbe de ralentissement, en tragant la courbe de la puissance Pf s + P mecQ en 
fonction du carre de la tension d’ alimentation, on obtient en principe une droite dont le 
prolongement jusqu’a la tension nulle donne les pertes mecaniques. 



Tension (V ) x io 4 

Figure C.4 Determination des pertes mecanique 
A partir de la figure C.4 on trouve P mec = 40W. 



Figure C.5 Vitesse de rotation pour Fessai de ralentissement 
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Le moment d’inertie J peut etre calcule par : 


J = 


40 


a 


' dCl 

v dt 




7n=n„ 


Q 


AD 

At 


154 


188 


0.0097 kg.rn 


Le moment d’inertie calcule presente le moment d’inertie totale c’est-a-dire de la machine, 
du frein a poudre et de la generatrice tachymetrique. 

Les parametres electriques et mecaniques sont regroupes dans le tableau ci-dessous. 


Constructeur 

LANGLOIS / SIEMENS MAS 32 

Puissance nominale P n 

3 kW 

Nombre de paires de poles p 

2 

Inductance cyclique statorique L s 

0.184 H 

Inductance cyclique rotorique L r 

0.184 H 

Inductance cyclique mutuelle L m 

0.167 H 

Resistance statorique R s 

1.310 Q 

Resistance rotorique R r 

2.013 D 

Moment d’inertie J 

0.0097 kg.rn 2 


Tableau C.3 Parametres electriques et mecaniques de la MAS 





